Preprint typeset in JHBP style - HYPER VERSION 



Five Easy Pieces: The Dynamics of Quarks 
in Strongly Coupled Plasmas 



^ I Mohammed Mia, Keshav Dasgupta, Charles Gale, Sangyong Jeon 

, Ernest Rutherford Physics Building, McGill University, 



=3 

1—5 

o 



o 



3600 University Street, Montreal QC, Canada H3A 2T8 
miam, keshav, gale, jeonShep. physics. mcgill.ca 



^ _ Abstract: We revisit the analysis of the drag a massive quark experiences and the 

D \ wake it creates at a temperature T while moving through a plasma using a grav- 



ity dual that captures the renormalisation group runnings in the dual gauge theory. 
Our gravity dual has a black hole and seven branes embedded via Ouyang embed- 
\ ding, but the geometry is a deformation of the usual conifold metric. In particular 

\ the gravity dual has squashed two spheres, and a small resolution at the IR. Using 

' this background we show that the drag of a massive quark receives corrections that 

^ _ are proportional to powers of log T when compared with the drag computed using 

O \ AdS/QCD correspondence. The massive quarks map to fundamental strings in the 

dual gravity theory. We use the perturbation produced by these strings to compute 
the wake and compare with the results obtained using AdS/QCD correspondence. 
We also study the shear viscosity in the theory with running couplings, analyze the 



X 

^ ' viscosity to entropy ratio and compare the result with the bound derived from AdS 

backgrounds. In the presence of higher order curvature square corrections from the 
back-reactions of the embedded D7 branes, we argue the possibility of the entropy to 
viscosity bound being violated. Finally, we show that our set-up could in-principle 
allow us to study a family of gauge theories at the boundary by cutting off the 
dual geometry respectively at various points in the radial direction. All these gauge 
theories can have well defined UV completions, and more interestingly, we demon- 
strate that any thermodynamical quantities derived from these theories would be 
completely independent of the cut-off scale and only depend on the temperature at 
which we define these theories. Such a result would justify the holographic renor- 
malisabilities of these theories which we, in turn, also demonstrate. We give physical 
interpretations of these results and compare them with more realistic scenarios. 
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1. Introduction 

There is no question that understanding the behavior of many-body QCD in the 
strong coupling regime is a hard problem to solve, but this needs to be done. The 
wealth of intriguing experimental data obtained at the Relativistic Heavy Ion Col- 
lider (RHIC) has made this situation abundantly clear. One of the main goals of 
the RHIC program is the creation and the analysis of the quark gluon plasma, a 
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new phase of matter predicted by lattice QCD. The goal of this paper is to discuss 
several quantities that can be related to observables measured at RHIC, or to be 
measured at the LHC Calculations are done in a regime where the gauge sector is 
nonperturbative, using techniques borrowed from string theory. 

From the early days of the RHIC experiments, the appearance of a strong ellip- 
tic hydrodynamic flow was taken as a consequence of early thermalization (before 
1 fm/c), and indicative of QGP formation [jl|. The elliptic flow, defined as the sec- 
ond harmonic component of the momentum distribution, develops when the system 
undergoing the hydrodynamic expansion has an elliptical shape with different short 
and long axes. This difference in the spatial shape causes difference in the pressure- 
gradient, which in turn causes the particles to accelerate more in the short axis 
direction. The anisotropy in the acceleration then causes the final momentum distri- 
bution to be anisotropic. The efficiency of this process of generating the momentum 
space anisotropy from the spatial anisotropy, however, depends on the size of the 
shear viscosity 77: a quantity that will be discussed here in some detail. Another 
experimental observable that has been linked to the formation of a plasma of quarks 
and gluons is the amount of energy lost by a fast parton travelling through this hot 
and dense medium. This phenomenon has also been dubbed jet quenching. Interest- 
ingly, there is a theoretical link between the concept of a small shear viscosity and 
that of a large jet quenching f^. The hard partons that travel through the strongly 
interacting plasma may also leave a wake behind, owing to the medium's response to 
the source which is the hard jet. We shall also discuss this in this paper. Somewhat 
related, the amount of energy lost by a heavy quark has been of great interest as well: 
the drag force can be related to the properties of the strongly interacting medium. 
Interestingly, the kinematics of bound states with a heavy quark can be observed 
through semileptonic decay channels. We compute the drag force experienced by a 
quark, as it looses energy to the surrounding medium. A brief introduction to some 
of these topics follows, before the general organization of the paper is outlined. 

1.1 Shear viscosity and the viscosity to entropy bound 

The shear viscosity represents the strength of the collective interaction between the 
two laminally flowing layers. Roughly speaking, large shear viscosity means faster 
mixing of the particles in two neighboring laminas. Somewhat counter-intuitively, 
the strength of this collective interaction is actually smaller when the microscopic 
interaction is stronger. This is because the rate of mixing is controlled by the mean 
free path. When the mean free path is small compared to the typical size of the 
flow velocity variation, two laminas with different flow velocities cannot easily mix 
since the exchange of particles is limited to the small volume near the interface of 
the two laminas: Most particles in the fluid just flows along as if there is no other 
layers nearby. On the other hand, if the mean free path is comparable to the typical 
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size of the flow velocity variation, then mixing between different layers can proceed 
relatively quickly. 

When the elliptic flow develops, the fluid has anisotropic fluid velocity distribu- 
tion. Since the shear viscosity controls the mixing, a large shear viscosity can quickly 
wash out these difference in the local fluid variables. Therefore, one can say that the 
smaller the shear viscosity, the stronger the elliptic flow. 

In the weakly coupled QCD, the ratio t]/ s is parametrically large since it is of the 
order l/a^ln(as) 0. If one is to believe that the QGP created in relativistic heavy 
ion collisions is in the weak coupling regime, one would then expect the elliptic flow to 
be small. One of the big surprises from the RHIC experiments is that this expectation 
is almost maximally violated. It turned out that the ideal hydrodynamics where the 
value of ?7 is just set to consistently describes the elliptic flow at RHIC pointing to 
a strongly interacting plasma well above the phase transition temperature. In fact, 
the QGP created at RHIC behaves like the most perfect fluid ever observed. 

The analytic tools available to a theorist to tackle this issue, however, have been 
rather limited. Perturbation theory is obviously not valid in the strong coupling 
regime and the lattice QCD study has been so far limited to the Euclidean space 
where extraction of the dynamic quantities such as the perfect-fluidity can be difficult. 

This situation changed dramatically when Policastro, Son and Starinets [^] dis- 
covered that Maldacena conjecture p can be used to calculate the shear viscosity 
of a certain strongly coupled quantum field exactly. More intriguingly, it turned out 
that in their solution the shear viscosity and the entropy density ratio, 77/5, had a 
minimum value at = ^ IS where fc^ is the Boltzmann constant. The authors 
P then made a conjecture that this is indeed the lower bound for any quantum field 
theory which has a gravity dual. Whether l/47r is a true bound or not have been 
debated many times in the literature J^, |], ^. But in all similar calculations, rj/s 
remains to be (9(l/47r) |^. The question is thus: What can we learn from this about 
strongly coupled QCD? 

In this paper, we use the non- AdS/QCD theory to investigate this question 



using a modified Klebanov-Strassler construction [0, ^ In the IR limit, 
this theory eventually leads to a confining SU{M) theory with a large M and Nf 
flavors^. This is, of course not a full QCD, but it shares many features with the 
real world QCD including the appearance of the renormalisation scale^. There is no 
gravity-dual construction of QCD yet but our study may shed some light on how the 
strongly coupled QCD should behave. 



^An alternative possibility is to get an approximately conformal theory with Nf flavors under 
cascade. This can also be realized in our set-up but we will ignore this possibility. 

^Note however that by QCD we will always mean a theory that resembles at IR large N QCD 
(but UV strongly coupled and almost conformal) throughout the text because it has a large number 
of colors. For finite N there is no known gravity dual. We will however need to keep finite number 
of fundamental flavors because of certain constraints that will become clearer as we go along. 
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Before we go into the details of our calculations, it is instructional to consider 
the physical meaning of the shear viscosity and its behavior and also why it makes 
sense to talk about the ratio rj/s as a measure of its strength. 

The definitions of the shear viscosity 77 and the bulk viscosity C, are given by the 
following constitutive equation 

= (v.(i;°) + VAT?) - l^i^'^iiT'')) - j^S,,Vi{T^') (1.1) 

where ST^^ is the deviation from the ideal fluid stress tensor in the fluid rest frame 
and e and P are the local energy density and the pressure, respectively. Upon using 
the thermodynamic identity Ts — e + P where s is the entropy density, the two 
coefficients can be also written as rj/Ts and C/Ts. Since the temperature is the 
only relevant energy scale in the highly relativistic fluid, one can easily see that the 
importance of the viscous terms depends on the size of the dimensionless ratios rj/s 
and C/s- 

To see what the shear viscosity signifies, consider situation where the fiuid is 
fiowing in the z direction but the speed of the fiow varies in the x direction (c.f. 
figure 1). That is, we have (T°^) = (T°^) = everywhere, but {T^^{x)) does not 
have to vanish everywhere. Recall that the stress part, T*-^, of the stress-energy 
tensor has the interpretation of the i-th component of the current for the conserved 
momentum density T'^-' . 

If the value of the conserved density {T^^{x)) at two different points are different, 
then there must be a net current between these two points. Now consider a z-y plane 
at a fixed x and think about the amount of microscopic current across this plane (c.f. 
figure 1). A particle crossing this plane from above in figure 1 has the average 
according to the thermal distribution and average according to the fiow velocity 
above the plane. A particle crossing the plane from the below share the same (but 
opposite sign) Vx but the average Pz in this case is proportional to the flow velocity 
below the plane. The net flow of through this z-y plane is then 

Txz{x) ^ T^^x - l){\vx\)t^ - To,{x + 0(k.|)th (1.2) 

Now for this expression to hold, the particles must not suffer a coUision while at- 
tempting to cross the plane. Hence the point of origin for the border-crossing particles 
must not exceed the size of the mean free path. Hence, Taylor expanding the above 
expression to the first order and setting / to be the mean free path A yields 

)r^-X{e + P){\vx\)tAuz{x) (1.3) 

Here we used the fact that in case of local equilibrium T^^ = {e + P)u'^u\ We also 
chose the frame where Uz{x) — 0. Comparing with the constitutive equation, one 
sees that 

r}^ X{e + P){v)t^ (1.4) 
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Figure 1: Difference of flows for a particle moving above and below the x-plane. 



where we have substituted with the average speed {v) for this rough estimate 

and discarded all 0(1) constants. Upon using the thermodynamic identity Ts = 
e + P, this becomes 

- ~ XT{v)tt (1.5) 
s 

From the discussion above, it is clear that the shear viscosity controls the rate of 
the momentum diffusion in the transverse direction to the flow. In fact, using the 
constitutive equation, it can be easily shown that the momentum densities satisfy a 
diffusion equation with the diffusion constant give by 

D = = ^ ~ A(t;)th (1.6) 

e + P Ts ^ ' ^ ' 



From Eq.( |1.5| ), we can also argue the existence of the lower bound following [|1( 
The mean free path A is defined as the average distance between two collisions. 
This length scale cannot become arbitrarily small due to the uncertainty principle 
AxAp > 1/2. The distance between two collisions must be at least longer than 
the Compton wavelength l/m and the de Broglie wavelength 1/p. Since the factor 
T{v)th in Eg. ( p.. 51) can be though of as the typical momentum scale in the medium, 
the Tj/s ratio must satisfy 

- > Aow (1.7) 

s 

where i?iow is a non-zero 0{1) constant. Therefore, the fact that rj/s is bounded from 
below by an C(l) constant is a rather robust consequence of quantum mechanics as 
long as the entropy is not completely dominated by large chemical potential; cf. Ref. 



2J]. What the value of is and whether the strongly coupled QCD has the r]/s 
ratio close to this bound, of course, is the issue that we are discussing here. 

In the theory we are using in this paper it will turn out that rj/s is smaller than 
the 1/Att bound of pj. Of course it was argued earlier in M that if one incorporates 
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corrections (where R is the curvature tensor) then the bound is automatically 
violated in the AdS space. What we find here is that if one goes beyond the AdS 
case by incorporating RG fiows in the gauge theory the bound is violated from the 
curvature square corrections but the non-trivial RG fiow do contribute to this also. 
We will discuss this in full details in sections 3.4 and 3.5 where we will argue that 
the required curvature squared corrections can be achieved by adding appropriate 
number of D7 branes to the background. 

There are two main ingredients in our calculation of the 7]/s ratio. The calcula- 
tion of the shear viscosity in a quantum field theory relies on the Kubo formula 

r/= ^Ihnlim-^ j dtd^x e'^'-''^-^ {[7r,jit,^),r^ (0)])^^ (1.8) 

where is the spin-2 part (traceless part) of the spatial stress tensor and the 
average is taken in the fiuid rest frame. This Kubo formula relies only on the infra- 
red behavior in the linear response theory. Hence, it is valid for an arbitrarily strong 
coupling limit. 

For the entropy, we may use the Beckenstein-Hawking entropy as originally done 
by Policastro, Son and Starinets |Q. Alternatively, we can directly calculate the 
entropy by using the thermodynamic identity Ts = e + P. When local equilibrium 
is reached, the energy density and the pressure can easily be obtained once Tf^" is 
known. In the fiuid rest frame, we have 

(T°°(x))eq = e(x) and ^{Tiix)U = Fix) (1.9) 

so that the product of entropy density and the equilibrium temperature is 

Ts = {T'\,+ ^-{TiU (1.10) 

We expect that the above result matches order by order in gsNf with the Bekenstein- 
Hawking result where Qs is the string coupling and Nf is the number of fiavors (we 
will discuss this in some details in section 3.5). Once we combine these ingredients 
with the RG fiow and curvature squared corrections the rj/s ratio becomes smaller 
than the known 1/Att bound. How small the ratio becomes seems to depend on the 
details of the model as well as on the UV degrees of freedom in the corresponding 
dual gauge theory, although the latter dependence is exponentially suppressed. 

1.2 Background geometry and holographic renormalisation 

Our aim in this paper is to study aspects of large thermal QCD using the dual 
gravity picture. Of course, as we mentioned earlier, there is no known dual gravi- 
tational background to thermal QCD. What we know is a gravitational background 
dual to a theory with a RG fiow that confines in the far IR and has infinite inter- 
acting degrees of freedom at the far UV. At zero temperature and in the absence 
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of fundamental flavors, this is the Klebanov-Tseytlin background ||29| with zero de- 
formation parameter. At high temperature and in the presence of Nf fundamental 
(and M bi-fundamental) flavors the background is much more complicated that, as 
far as we know, has not been studied before (see for a different model that in- 
corporates flavors but doesn't have a good UV behavior). In sections 3.1 and 3.3 
we show that, to lowest order in QsNf and QgM'^/N, there is some analytic control 
on the background, meaning that we can derive analytic expressions for the metric 
and fluxes to the lowest order in QgNf and QsM'^/N. To higher orders in QsNf and 
QgM'^/N one can only derive the expressions for the metric and fluxes numerically. 
We clarify many previously ignored subtleties in the literature, namely the existence 
of small resolution and the effect of this on the fluxes. 

Once we have infinite interacting degrees of freedom in the far UV we naturally 
face the question of holographic renormalisation. For standard AdS background this 
has been demonstrated beautifully in the series of papers [1^, [1^, ^ |T6|, |17[. For 
the Klebanov-Strassler background without fundamental flavor, an equivalent treat- 
ment has been shown to apply in |TB[. In section 3.3 we show that, modulo some 
subtleties, such a treatment of holographic renormalisation can almost be extended 
to theories with fundamental flavors also. The subtleties have to do with the exis- 
tence of non-trivial powers of log r over and above the 1/r suppressions in various 
expressions (here r is the radial coordinate that determines the energy scale of the 
gauge theory). In the limit of small QsNf and QgM'^/N we argue that once the high- 
est integer power of r has been regularised, the theory is naturally holographically 
renormalised. On the other hand once QgNj and QsM'^/N are large, we loose all 
control on our order-by-order expansions and new methods need to be devised to 
holographically renormalise the theory. 

1.3 Organization of the paper 

The paper is organized as follows. In section 2 we summarize our main results. This 
section is meant for readers who would like to know our results without going through 
the details of the derivations. The summary section also involves a discussion of the 
subtleties of background associated with RG flows etc. In the presence of non-trivial 
RG flows and corresponding Seiberg dualities the interpretation of gauge/gravity 
duality here is not so straightforward as in the AdS/CFT case. We point this out in 
some details. 

Section 3 is the main section of our paper. As promised in the title, we perform 
five easy computations in thermal QCD. In sections 3.1 and 3.3 we give a detailed 
derivations of our background. At zero temperature and in the presence of funda- 
mental and bi-fundamental flavors the gravity dual is given by Ouyang [^. However 
inserting a black hole in the Ouyang background to generate a non-zero temperature 
in the gauge theory changes everything. We can no longer argue that the fluxes, 
warp factor etc would remain unchanged. Even the internal manifold cannot remain 
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a simple conifold any more. All the internal spheres would get squashed, and at 
r = there could be both resolution as well as deformation of the two and three 
cycles respectively. In section 3.1 we present our results to 0{gsNf, ggM"^ /N), and 
in section 3.3 we give more detailed derivations and extend this to higher orders 
in g.sNf and QsM'^/N. In the limit where the deformation parameter is small, we 
show that to 0{gsNf,gsM^/N) we can analytically derive the background taking a 
resolved conifold background. The resolution parameter depends on ggNj, ggM"^ /N 
as well as on the horizon radius Vh- 

Our background also has D7 branes that give rise to fundamental flavors in the 
gauge theory. These D7 branes are embedded via Ouyang embedding ( |3.8|) and 
therefore the strings with one ends on the D7 branes and the other ends falling 
through the horizon would be associated with Nj thermal quarks. In section 3.2 
we evaluate the mass and drag of the quarks in this background. We point out there 
how our results differ slightly from the analyses done using AdS/QCD techniques. 

In section 3.3 we give a more complete picture of the system. Although the 
section is dedicated to evaluating the wake of the quarks in a thermal medium, we 
study three related topics here. The first one is already mentioned above: we give 
a detailed derivation of the background geometry that fills up the gaps left in the 
discussions of section 3.1. We then give a detailed derivation of the holographic 
renormalisability of our theory in the limit of small gsNf and gsM'^/N. Finally all 
these analyses are combined to evaluate the finite energy-momentum-tensor of the 
background plus the quark strings. From there we evaluate the wake of the quarks 
by removing the energy-momentum-tensor contribution of the quark strings. Due to 
the complicated nature of the background we could evaluate certain formal quantities 
in this section without going into numerical details. In the appendices we provide a 
toy example with only diagonal metric perturbations, where a more direct analysis 
could be performed^. 

Another upshot of this section is the realization that we could study infinite 
number of gauge theories by cutting off the geometry at various r = and UV com- 
pleting them by inserting "UV caps" at various Vc- From gauge theory side this is like 
inserting correct relevant, marginal or irrelevant operators; and from the gravity side 
this is like inserting non-trivial geometries from r = to r = oo. We show that any 
thermal quantities evaluated in these gauge theories are completely independent of 
the cut-off scale r = Tc and only depend on the temperature (the dependences on the 
far UV degrees of freedom are exponentially suppressed), justifying the holographic 
renormalisabilities of these theories. 

The cut-off independence is also apparent in the last two sections, sections 3.4 
and 3.5 where we study viscosities and the ratios of the viscosities by their corre- 
sponding entropy densities. For both these cases there are no cut-off dependences 

^See for example Appendices A, B and C. Appendix A is not directly related to the main 
calculations of our paper, but gives the corresponding example for the AdS case. 
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but the ratios of the viscosities by their corresponding entropy densities are always 
smaller than the celebrated l/An bound 0. 

Finally in section 4 we present our conclusions. Throughout the paper, we 
allude to various future directions that will be dealt in the sequel to this paper . 

2. Summary of the results 

Before we go into the full analysis, let us summarize the main results of our paper. 
This summary is meant for readers who would want to see the main conclusions 
without going through our detailed calculations. 

In this paper, as the title suggests, we have done five concrete calculations. They 
are, in order of appearance: 

• Construction of the gravity dual of a thermal gauge theory with running coupling 
constant and fundamental flavors. 

• Mass and drag of the quark in the gauge theory from the gravity dual. 

• Wake left by the quark when it moves in the QGP medium, again from the gravity 
dual. 

• Shear viscosity of the QGP medium; and finally 

• Viscosity by entropy bound for the quark from our dual picture. 

Let us now give some brief descriptions of each of the five calculations. First is 
the gravity dual. Recall that most of the recent analysis have relied heavily on using 
AdS/CFT correspondence to study the above phenomena. However in this paper we 
would like to study the gauge theory with non-trivial RG flow, in the hope that we 
can capture some aspects of the strongly coupled QCD. Although our gravity dual 
will be far from realistic QCD, we will try to address certain generic issues that could 
lie in a class of gauge theories that are at least IR confining. 

At high temperature the situation is a little subtle, but still IR dynamics of 
QCD could not be extracted from the AdS-Black Hole (AdS-BH) picture. Here, as 
one may recall, a black hole is inserted in the AdS picture to account for a non-zero 



temperature in the gauge theory [11]. To do a better job we need another model 
that can capture the RG flow in the gauge theory. 

The model that comes to mind immediately is the so called Klebanov-Strassler 
warped conifold construction [0^. The advantage of this construction is that the 
gravity dual — which is a warped deformed conifold with three form type IIB fluxes 
— captures the RG flow of the gauge theory. The KS gauge theory is confining in the 
far IR but is not asymptotically free. The UV of the theory is a more complicated 



cascading gauge theory. The other two duals [^] and do not have a good 



four- dimensional UV description. Additionally, all these constructions are for zero 
temperature gauge theories. 

^ A somewhat equivalent constructions given around the same time are , p6[ . 
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The other issue with the KS picture is that the quarks therein are all in the 
bi-fundamental representations of the two possible UV gauge groups; and they even- 
tually cascade away in the far IR. So what we need is a dual gravity theory that 
allows fundamental quarks at high temperature. 

Before inserting fundamental quark at high temperature let us point out that 
the situation now, even at zero temperature (and without fundamental flavors), is 
much more subtle than the AdS case studied earlier. Due to renormalisation group 
flow in the gauge theory side, we need to ask precisely what aspects of the dual 
gravity picture captures the dynamics of the strongly coupled gauge theory. The key 
question to ask here is whether the weakly coupled gravity dual sees the cascading 
Seiberg dualities or it only sees a smooth RG flow in the gauge theory side. Since 
the answer to this question lies in the heart of the matter, we will spend some time 
elaborating this. We would like to caution the readers that this point has been 
misunderstood in most of the literature and the only article, in our opinion, that has 
been able to fully explain the subtleties is the one by Strassler P^. What follows 
below is an elaboration of these subtleties. 

To start off let us ask what JV = 4 AdS / CFT duality tells us? Recall that due to 
the tight constraint from supersymmetry all quantum corrections in A/" = 4 theory 
is cancelled out, leaving us with only classical amplitudes. This in particular means 
that if we choose any gauge coupling in the theory, it stays there, with possible 
small finite shifts, under any RG flow (which here means going from UV to IR). For 
example if we choose the gauge coupling to be very strong, the theory will remain at 
strong coupling from UV to IR. On the other hand once the gauge theory is at strong 
coupling we could analyze the theory from weakly coupled supergravity description 
on AdS space. However when the gauge theory is at very weak coupling there exists 
no weakly coupled gravity description, and the story there is captured by the full 
string theory on AdS space^. 

The above conclusions also mean that every value of the coupling is a RG fixed 
point in the gauge theory. There is no coupling flow and therefore the system is 
simple without any inherent subtleties. What happens now if we introduce a non- 
trivial RG flow in the theory? There are three cases to consider: 

• The RG flows lead to a non-trivial fixed point (or isolated fixed points) in the 
theory. 

• The RG flows lead to a non-trivial surface of fixed points in the theory; and 

• The RG flows lead to no fixed point (s) or fixed surface (s) in the theory. 

The first case, as far as we know, leads to no known gravity dual so will ignore this 
case. The second case is more interesting. We know of one example where non-trivial 
RG flows in the theory lead to a fixed RG surface. This is the so-called Klebanov- 
Witten model There are three couplings in the theory {gi,g2, h) corresponding 

^Alternative one could think that the weakly coupled A/" = 4 gauge theory is captured by string 
theory in twistor space. This is basically the key essence of p8|. 
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to the two gauge couplings and the couphng associated to the quartic superpotential 
of the theory respectively. The gauge group is SU{N) x SU{N) and the three beta 
functions are: 



/5, 



91 




/3, 



92 




/3, = 7^(1 + 270) (2.1) 



where rj = hfi is the dimensionless coupling of the theory for any energy scale /i. As 
discussed in details in |2^ all the fields in the theory have the same anomalous 
dimension 7o(fi'i,fi'2; h), and therefore the three beta functions vanish exactly when: 



'ioigi,g2,h) 



1 

2 



(2.2) 



which is one equation for the three couplings. Therefore the fixed points in this 
theory form a two-dimensional surface in the three-dimensional space of couplings 
(see figure 2 for details). This also means that the typical RG flow in the theory 




Figure 2: The two-dimensional RG surface in the Klebanov-Witten theory. 



will take the following simple form as illustrated in figure 3 below. Since the sign of 
the two beta functions are negative any arbitrary flow in the coupling constant space 
will bring us to the fixed point surface. This surface is IR stable. Notice also that at 
the boundary of the Klebanov-Witten fixed point surface the gauge theory is weakly 
coupled and so the gravity description is very strongly coupled. To get the weakly 
coupled supergravity description, we need to go towards the centre of the fixed point 
surface. The AdS x T^'^ description is valid here. Note also that on the surface the 
bi-fundamentals scalars all attain anomalous dimensions of — |, as obvious from the 
beta functions given earlier. 

The story changes quite a bit once we deform away from the AdS background. 
At zero temperature and in the absence of fundamental flavors, this is of course 
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the Klebanov-Strassler model [0. The gauge theory dual of the Klebanov-Strassler 
model is more complicated now because of the non-trivial RG flows of the two cou- 
plings lead to cascades of Seiberg dualities in the theory. Although the UV of the 
gauge theory has infinite degrees of freedom, the theory is holographically renormal- 
isable^ [|18] (see also 



To see how the story differs from the AdS case discussed above, let us take the 
far UV theory to be SU{kM) x SU{kM — M) with gauge couplings Qk^Qk-i for the 
two gauge groups respectively and rj to be the other dimensionless coupling defined 



above. The three beta function now are 12, 27 



gjkM 
' 167r2 



(l + 27o) + f(l-7o) 



87r2 



= 77(1 + 270) 



gU{k-l)M 
167r2 



(l + 27o)-^(l-7o) 

^ 87r2 



(2.3) 
(2.4) 



where we see now that they differ from (|2.1| ) by 0{l/k) factors. This also means 
that there will be no point in the coupling constant space where all the three beta 
functions could vanish exactly although there would be numerous points where all 
the beta functions could be very small. Two questions arise immediately: 

• Since there are infinite number of gauge theories involved here, what is the gravity 
description of the theory? 

• At the point where the theory has a weakly coupled gravity description, is the 
gauge theory described by a smooth RG flow or the description involves "choppy" 
Seiberg dualities? 



^We will discuss later that the theory with fundamental flavors is also holographically renormal- 
isable. 
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Since the answer to these two questions are rather involved, we will go in small 
steps. It is clear that the geometrical description that involves cascade of Seiberg 
dualities cannot have a weakly coupled gravity description. The RG flows in the 



theory have been succinctly presented in |27] so we will be brief here. The RG 
flows at the boundary of the Klebanov-Witten fixed point surface can be illustrated 
by figure 4. Observe that the flow takes us from one surface to another because 
we are moving from one set of gauge theory descriptions convenient at a certain 
scale to another set of descriptions convenient at a different scale. For the RG 




Figure 4: The RG flow at the boundary of the Klebanov-Witten wall for the Klebanov- 
Strassler model. The big dots are the Seiberg fixed points of the theory. The vertical 
distances between the planes don't signify anything here. 



flows illustrated in figure 4 there is no weakly coupled gravity description available. 
Therefore, contrary to popular belief, the usual cascade of Seiberg dualities do not 
have a supergravity description on a deformed conifold. In fact there is no simple 
gravity description that could capture the choppy RG flows in the dual gauge theory! 
So then where would our gravity calculations fit in? In other words, the analysis that 
we do in the gravity side captures what aspects of the gauge theory? Clearly since 
we fail to capture the cascades of Seiberg dualities using our weakly coupled gravity 
picture, do we then see any interesting aspects of the gauge theory now? 

Things are not really that bad once we realise that the strongly coupled gauge 
theory description or equivalently the weakly coupled gravity description becomes 
better as we go away from the boundary of the Klebanov-Witten fixed points surface. 
As we go towards the centre of the surface, the RG flows lose their choppy nature and 
tend to become smooth. This is illustrated in figure 5. Thus in the gauge theory 
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Figure 5: The RG flow as we go inside the boundary of the RG fixed point surface. 
Observe that the gauge theory becomes strongly coupled at the centre of the surface and 
the RG flow loses its sharp edges and subsequently becomes very smooth. 



side we can summarize our situation by three points: 

• The theory lies at the boundary of the RG surface and jumps from one Seiberg 
fixed point to another^. This is the "usual" cascading and the surface corresponds to 
the two-dimensional surface of fixed points (for the corresponding Klebanov-Witten 
theory) in a three-dimensional space of coupling constants. 

• The theory lies close to the boundary of the surface but never quite touches the 
fixed points. So this flow is parallel to the boundary. 

• The theory hovers at the centre of the two-dimensional surface and has a smooth 
RG flow from UV to JR. The degrees of freedom of this theory changes continuously 
from UV to IR (and is a well defined QFT). 

Similarly in the gravity side we can also summarise the situation by three "dual" 
points: 

• For the first case the gauge theory is weakly coupled and therefore there is only 
a strongly coupled gravity description. This amounts to saying that the tree-level 
dual is full string theory on a warped deformed conifold (which also means that we 
cannot say anything from the gravity side other than some protected quantities!). 

• For the second case a somewhat similar statement can be made. There is no weakly 
coupled gravity description available. 

we observed above, all the three beta functions do not vanish simultaneously. Therefore the 
fixed points are not the absolute Klebanov-Witten type fixed points, rather they have s small tilt 
given by 0{l/k) corrections ( |2.3D . 
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• For the third case there is a well defined weakly coupled gravity description. This 
is supergravity defined on warped deformed conifold and tells us clearly that in the 
dual QFT the colors should decrease logarithmically (although its a little tricky to 
say how many colors we have). 

Our calculations are done for 

this regime where there is 
only a smooth RG flow with a 
weakly coupled gravity description 



'k-3 



'k-1 




gk A 
— • 



'k-2 




gk-3 



gk-1 



gk-2 



Cascade of Seiberg dualities 
here, captured by strongly 

coupled gravity 



Figure 6: The complete RG flows depicting the regime of validity of our calculations. 
This figure is taken from [^^. The yellow wavy hne is the branch cut, and corresponds to 
the various planes in the earlier figures. 



All of our above discussions imply that the regime of smooth RG fiow in the gauge 
theory side can be captured by weakly coupled supergravity description. But what 
supergravity description are we looking for here? We want our supergravity descrip- 
tion to capture the RG fiow in a gauge theory with fundamental fiavor and at a 
non-zero temperature. In other words we want the sugra description for a thermal 
gauge theory with fundamental fiavors. 

The first problem of having fundamental quark may be solved by inserting Nf 
D7 branes in the KS geometry. However this is subtle as we discuss in sec (3.1). 
What we know so far is how to insert coincident D7 branes in the Klebanov-Tseytlin 
background . This is the Ouyang background |^ that has all the type IIB fiuxes 
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switched on, including the axio-dilaton. Our aim here is to insert a black hole in this 
geometry and study the resulting gravity dual. 

Inserting a black hole in the Ouyang background changes everything. We can 
no longer argue that the fluxes, warp factor etc would remain unchanged. Even 
the internal manifold cannot remain a simple conifold any more. All the internal 
spheres would get squashed, and at r = there could be both resolution as well as 
deformation of the two and three cycles respectively. Our metric then should look 
like (see the definitions of the coordinates in sec (3.1)): 



1 



gi{r)dt^ + dx'^ + dy"^ + dzA + v^L(r)~Mr2 + dMf\ (2.5) 



where we have been able to work out h{r),gi{r) and dAi^ only in the limit where 
the resolution factor of the internal conifold is a constant but the deformation factor 
is zero^. Our results are: 

• The warp factor h{r) is given by (|3.35|) . 

• The black hole factors gi{r) are given by ( |3.13|) 

• The internal manifold is given by a warped resolved conifold. 

In certain limits along with Qi = g2 = 9, that we discuss in the text, we can show 
precisely how our metric differs from the AdS case (see footnote 37 for definitions): 



ds' 



— i-gdt^ + dx'dxi) H dr"^ + dM 



(2.6) 



(a log r + B log^r) 



— r (-gdt^ + dx'dxi) dr^ + dM 

gr^ 



where the first line is the standard AdSs x T^'^ space with radius L = (AngsNy^^, 
and the second line is the deformation of the AdS as well as the internal spaces. 
Precisely because of this deformation, as we mentioned above, even the fiuxes change 
drastically from the ones given by Ouyang. Our fiuxes capture the fact that the two- 
cycles in the internal warped resolved conifold are slightly squashed. The three-form 
RR fiux is (the coordinates {Oi.cpi.ip) parametrise A^s in ( p.5|) ): 



F3 = 2MAi [ 1 + %^ log r ) A ^ (sin 9^ d9i A d(pi - 5isin ^2 d92 A 

ZTT J 2 

3gsMNf ^ dr ( ^2 . . , , 

^-A2 — A ew, A cot — sm 6*2 c/<^ 

47r r \ 2 

--^z sin 9\ sin 6^2 f cot d9\ + i?3Cot d92 ) A d(^\ A a^)2 



9 \ 
i?2Cot — sin 9\ d(f)i j 



(2.7) 



*Such details won't matter when we study only the metric. But fluxes will carry the information 
of background resolution. 
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and the three form NS-NS flux is: 



6gsA4^M 1 H log r H log sm — sm — — A - sm t^i dOi 

\ A-K 2t\ 2 2 1 r 2\ 

r. ■ n ^ ^dlMNf f dr 1 ' 

- Bism 02 dd2 A d02 + ^^5 — A - -de^ 



A 



9 
A I cot ^ d62 - B^coi dOi 



(2.8) 



where we differ from Ouyang analysis precisely by the asymmetry factors Ai^Bi. 
These asymmetry factors incorporate all order corrections and are given in ( p.83| ) 
(where we present only the first order terms). Finally, the axio-dilaton and the five 
form are given by ( p.31| ) and ( p.84| ) respectively. These are thus our first set of results 
related to the gravity dual on which we base all our subsequent calculations. Due to 
the complicated nature of our background, we have only worked things out to order 
0{gsNf, g^MNf) at small r, and provide some conjectural solution for large r from 
the full F-theory embedding of our solution. In fact for most part of this paper our 
calculations are valid in the following limits of Qs, M, Nf and A^: 



gs, gsNf, glMNj, ^ 0, 



{gsN, gsM) 



—7- oo 



(2.9) 



The precise way some of these go to zero or infinity is discussed towards the end of 
section 3.1. Note that large Nf limit is probably not realisable here because of the 
underlying F-theory constraints . 

Our second set of calculations are related to the mass M and drag z/ of a quark 
from the gravity dual ( p.6|) . The quark is identified in the dual theory with the string 
whose one end is connected to the D7 branes and the other hand falls in the horizon 
of the black hole. The drag of the quark is precisely the rate at which the string 
loses its momentum and energy to the black hole. This is given by^: 



(2.10) 



1 + 



2-kN 



log 



r 



(l-i;2)i/4 



1 + 



2lT 



log 



r 



(1-1,2)1/4 



where we see that we differ from the AdS result by the log T and log (1 — v"^) 
corrections where T = r^, the horizon radius, is the characteristic temperature and 
V is the velocity of the string (also the quark). These corrections aren't very big 



^Note that we will be working in the limit where H — c = ks ^ Mp = a' = 1. Thus all the 
relevant QCD observables would appear dimensionless throughout. Besides these choices of scales 
we will also impose tq = 1, where rg is the distance between the tip of the seven brane and the 
black hole horizon. This will be discussed more in section 3.2. 
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and occurs precisely because of the logarithmic running in our theory. The other 
coefficients appearing in (|2.1CI|) are defined in section (3.2). 

In the above calculations, we had taken the string stretching between D7 branes 
and the black hole to be a probe in the deformed AdS background. This means that 
the back reactions are completely neglected. This cannot be the full picture, because 
a moving string should create some disturbance in the surrounding media. This is the 
premise of our third set of calculations. Under some reasonable assumptions (that we 
elaborate in sec. (3.3) of the text) we can use this back reaction to calculate the wake 
of the quark in the gauge theory from our dual gravity background. This analysis is 
simple to state but involves not only a detailed series of manipulations that include, 
among other things, writing an action and then regularising and renormalising it to 
extract finite gauge theory variables, but also something much more elaborate and 
interesting. The interesting part is that not only we can do our wake analysis (or 
in principle all other thermodynamical quantities) for the dual geometry that we 
elaborated above, but also on the geometry that is cut-off aX r = and a non-trivial 
UV cap (or UV geometry) attached from r = to r = cxd. Such a procedure actually 
creates a new UV completed gauge theory at the boundary that is different from the 
original theory (that we called the parent cascading theory). This new theory and 
the original "parent" one differs by certain operators that we can define carefully at 
the cut-off. Thus the properties of the new theory are not the universal properties 
of the parent cascading theory. The parent cascading theory has infinite number 
of degrees of freedom at the far UV. Our new theory should also have very large 
degrees of freedom at the UV. We can call these degrees of freedom at the boundary 
(i.e at r = oo) as Nuv such that when Muv = e~",e — )■ 0,n >> 1 we are studying 
the parent cascading theory, and the other boundary theories have n> 1 (see figure 
10 for details). An elaborate discussion of this is given in section 3.3 which we 
would refer the readers for details. In fact the result of studying these gauge theories 
are both remarkable and instructive. The remarkable thing is that no matter what 
thermodynamical quantities we want to extract from these theories, the final results 
are completely independent of the cut-off r = Tc that we imposed to derive these 
theories. The results only depend on the characteristic temperature T (that we fix 
once and for all) and on the UV degrees of freedom via e~^^^. Since Afuv is always 
infinite for us the results are only sensitive to T. Our third set of calculations related 
to the wake demonstrates this in full details. No matter how involved are the UV 
descriptions, or the procedure to holographically renormalise these theories, the final 
answers for the energy-momentum tensors for the gauge theories are remarkably 
clean and can be stated as^°: 

rpmm _ f ^ ^ \ ^ j ( Zjrnn , rrnm \ (4) /I /' i ZV'^™'! 

medium+quark — J "^271^ ^-^ | ^ l"l 1°! ' I"! ~'~ l«l ' nn 

"'^^Repeated indices of the form a"6„, a^bnn are summed over n, but a„6„, (ab)™ etc may not be. 
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\a\ 



(5)[/3] 
nn 



oo 
j=0 



+ 0{re-^-^)} (2.11) 



from where the wake of the quark can be computed using the relation ( p.64[ ) which 
equivalently means we subtract the energy momentum tensors of the quark from 
(^TTT) above. All the functions appearing above are described in section (3.3), with a 



denoting the effects of the flavors. Here (b^"r^-y , Afuv) together specify the full boundary 



theory for a specific UV complete theory (for a discussion on the back reactions, see 
after ( |3.95|) ). Note that the result has no dependence on r = Vc- In the limit where 
■f^uv = e~", n » 1 we reproduce the result for the parent cascading theory, and for 
A/'ut, — !■ e~",n > 1 we have the UV dependences specified above. 

The instructive thing about our result is the way we differ from the AdS case. 
The energy momentum tensor of the system from AdS / CFT case is the first part of 



the above formula given by Smm*' (see (|3.115D for details). The rest is the deformation 



from the AdS result (we have to remove some constant pieces from ( p.ll| ) to account 
for the stress tensor properly). This way we can again see how different UV completed 
gauge theories can change some of these calculations. Of course all these theories do 
inherit some of the universal properties of the parent cascading theory. 

Our fourth set of calculations is related to studying the shear viscosity with and 
without curvature squared corrections. Without curvature squared corrections, but 
involving the RG flow, the result for shear viscosity r] is easy to state: 



3 r 2 



V 



N 



1 + E 



47r + i log^ (1 - r%-^""' 



(2.12) 



where ak are functions of T that are given in section 3.4. Note again that the result 
is independent of the cut-off at r = r^. In the limit where Af^v = oo the shear 
viscosity has a simple form. This is related to the AdS result also. On the other 
hand the viscosity to the entropy ratio, which is our final set of calculations taken in 
the limit where we switch on all the ingredients i.e the RG flows, curvature squared 
corrections as well as the contributions from the UV caps, has the following form: 



s 



C3K 



BojAn^ - log' Co) + 4nAo log Cp 
(^47r2-log2 a)' + 167r2 log' Co 



(2.13) 



where {Ao, Bo,Co) are given in ( |3.229|) and C3 is given in ( p.l63|) . We see that the 
ratio is again independent of any cut-off; and in the limit f/uv ^ e " and C3 — 
( p.l3|) the bound is exactly saturated. However in the limit C3 7^ we have a violation 
of the bound. The plot of the behavior of rj/s is given in sec (3.5). 
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3. Dynamics of quarks in strongly coupled plasmas 

This section contains the five elements alluded to in the title of this paper. It is 
now time to construct a more detailed scenario wherein certain aspects of QCD 
calculations could be performed. As we discussed in the introduction and in the 
summary above, most of the previous analysis relied on Anti-deSitter (AdS) spaces 
whose dual is a conformal field theory (CFT) with no running couplings. Many of the 
recent works in this field have focussed on AdS/CFT correspondence to study certain 
behaviors of QCD. Our aim in this section would be to use non-AdS backgrounds to 
study properties of QCD. 



3.1 Construction of the gravity dual 

It turns out if we embed D branes in certain geometric background, the gauge the- 
ory that lives on the D branes may become confining and exhibit logarithmic run- 
ning coupling. One of the most popular background to achieve this property is the 
Klebanov-Strassler background (see also [^5|, ^). In Klebanov-Strassler (KS) model 



D3 branes are placed at the tip of a six dimensional conifold and M D5 branes 
are placed in such a way that they wrap a 2-cycle on the conifold base T^'^. The D3 
and unwrapped part of the D5 branes extend in four spacetime directions orthogonal 
to the conifold. The SU{N) x SU{N+M) gauge theory living on these branes contain 
matter fields that transform as bifundamental color representation [N, N + M)^ and 
(A^, A^ + M)c of the group SU{N) x SU{N + M). Under a Renormalisation Group 
fiow the gauge group cascades down to SU (M) group in IR. The KG fiow and fixed 
points surface for this model have been shown earlier. In the regime where we have a 
smooth RG fiow, the theory confines in the far IR. Thus there is a small regime of the 
theory that gives us a weakly coupled gravity dual to a confining SU{M) theory. Of 
course since the UV behavior is not asymptotically free^^, KS model doesn't give us 
the full gravity dual for QCD. However it does come close in giving us at least a dual 
model that has a running coupling constant. Similar behavior can also be argued for 
[P^d and , although the UV behaviors of |^ and |^ are six-dimensional theories 
and may develop baryonic branches [B^. On the other hand, in far IR, one can show 



that many BPS quantities have one-to-one correspondences 



The original KS model (as well as |2^, |26|) do not have quarks in the fundamental 
representation of the gauge group. To introduce fundamental quarks, we need D7 
branes in the gravity dual. However this is subtle because the full global solution 
that incorporate back reactions of the D7 branes on the KS background has not yet 
been computed. What we have computed is the local metric that incorporates the 
deformations of the seven branes when these branes are moved far away from the 



11 



Plus it has infinite number of degrees of freedom. 
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regime of interest (see [34, 35 1 for more details). This is given by 



ds'^ = hi [dz + ai dx + aa dyf + h2 [rfy^ + dOj] + [rfx^ + h-^ dOl] + 

+ /i5 sin ip [dyi 6/6*2 + dy d6]\ + cos ip [d6id62 — fix dy] 
Bns = ^ dOi + h^^^ dy A d^2, Brr = -2A dy A dz, = = (3.1) 

where (x, y, z, 9i, ip, r) are the small local regime around the point {{(pi), {ip), {9i), tq) 
given in the following way: 



V^TqaT^' v^(7o"+4a2)y^ sin (6*2) 

2x 

01 = {(pi) + I , r = ro + , 

V7oV/io sin (6'i) V7oV^o 

^1 = (^1) + ^2 = (^2) + ^ (3.2) 

V7oV/io V (70 + 4a^)v^o 

where 7o(?"o) and /lo('^o) are some constant functions of ro (see ^ for details); 
and the un-tilded coordinates are used to write the standard metric of the conifold 
in the following way: 

rfs2 = ^^2 ^ ^2 I 1 \^{d9l + sin^^i d(p'f) + -(# + cos 9^ d(pi + cos ^2 d(p2)A{3.3) 

There are a few key issues that we want to point out regarding the metric ( |3.1| ): 

• First, since the seven branes are kept far away the axion-dilaton vanish for the 
background locally. Globally there will be non-zero axion-dilaton. 

• Secondly, observe that the two spheres (parametrised originally by {(pi,9i)) are 
replaced by two-tori locally. This issue has already been explained in details in [0^ . 
Furthermore the two tori do not appear with the same coefficients. In fact there 
is a squashing factor associated with the two tori. We believe that globally such 
squashing factor should also show up for the two spheres. 

• Thirdly, the local back reactions on the metric due to fluxes and seven branes 
would modify the warp factors once we go to the full global scenario. 

• Finally, the above local (and the subsequent global) picture is dual only to zero 
temperature gauge theory. What we need for our purpose is high temperature gauge 
theory. This would mean that the global extension of the above background (|3.1j ) 
should contain a black hole whose horizon size should correspond to the temperature 
in the dual gauge theory at far UV. Since the number of effective degrees of freedom 
are changing as we go from UV to IR, the entropy and temperature will depend 
crucially on what UV degrees of freedom we are considering at a given cutoff. This 
is a subtle issue and we will discuss this carefully a little later. 
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Thus putting a black hole in the global metric is non-trivial because of all the 
above considerations. However we can formally write the metric in the following way 
(although we will use Minkowski coordinates throughout, unless mentioned other- 
wise, one may write everything in Euclidean coordinates also): 



1 

7^ 



{-gidt^ + dx^ + dy'^ + dz^) + g2^dP + r'^dMj 



(3.4) 



where {x,y,z) not to be confused with (x, y,z) discussed above, gi are functions^^ 
that determine the presence of the black hole, h is the lOd warp factor that could be 
a function of all the internal coordinates and dA4l is given by: 



dMl 



hi{d^p + cos 6i dtpi + cos 62 c?02)^ + h2{d6l + sin^6'i 



+h4{h3del + sin^6'2 



/15 cos ip {d9id92 — sin 6'isin 6*2^01 



02) + ^5 cos ip {d9id92 — sin 6'isin ( 
-|- sin %l) (sin 9i d92d(f)i — sin 92 d9 







(3.5) 



with hi being the six-dimensional warp factors. One advantage of writing the back- 
ground in the above form is that it includes all possible deformations in the presence 
of seven branes and fluxes. The difficulty however is that the equations for the warp 
factors hi are coupled higher order differential equations that do not have simple 
analytical solutions. The original KS solution is in the limit 



h. 



9i 



1, hi = fixed 



(3.6) 



but has no seven branes. In the presence of seven branes we can do slightly better 
in the limit: 



/i.^ = 0, /i.ci = 1, hi - ho 



9i 



(3.7) 



which puts a seven brane in a resolved conifold background with a = constant 
using the so-called Ouyang embedding [^. A black hole could be inserted in this 
background by switching on non-trivial g^. However a naive choice of fiuxes in (|3.71 ) 



is known to break supersymmetry p^ . 

Our next choice would then be to go to the limit where the resolution parameter 
a in ( p.7| ) is vanishing. This is the Ouyang background with gi{r) = 1 i.e with 
seven branes but no black holes. Even in the absence of black holes, supersymmetry 
is an issue here as was pointed out in [^. Supersymmetry is spontaneously broken 
but could be restored by switching on appropriate gauge fiuxes on the seven branes 

my. 



The seven branes are embedded via the following equation (see also PO 



rzexp 



- 01 - 02) 



. 9i , 92 
sm — sm — = u 
2 2 ^ 



(3.8) 



12 



They would in general be functions of (r, 9i). We will discuss this later. 
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where /x is a complex quantity. In the hmit where /x — )■ 0, the seven branes are 
oriented along the two branches: 



Branch 1 : = 0, (f)i = 

Branch 2 : 02 = 0, 02 = (3.9) 

From above it is easy to see that the seven branes in Branch 1 wrap a four cycle 
{d2i<p2) and {ip^r) in the internal space and is stretched along the spacetime direc- 
tions {t,x,y,z). Similarly in Branch 2 the seven branes would wrap a four-cycle 
(6'i,0i,r,V'). 

From the above discussion, one might get a little concerned by the fact that 
the seven branes wrap a non-compact four cycle in the internal space. This would 
suggest a violation of the Gauss' law as the axion charges of the seven branes have 
no place to escape. For the time being this apparent paradox can be resolved by 
allowing the seven brane to wrap a topologically trivial cycle so that it would end 
abruptly at some r = tq when the embedding is (|3.8| ). This is similar to the seven 
brane configuration of We will discuss a better embedding later. 

Once the above issues are resolved, we can insert a black hole in the modified 
Ouyang background (|3.7|) by switching on appropriate gi{r). Clearly we do not expect 
hi to remain constant anymore. We also expect M and Nf to be given by some M^s 
and Nj^ respectively. Our first approximation would then be to take the following 
ansatze for the hi, Mes and NJ^: 

h, = l + 0{g,), h2 = h, = ]: + 0{gs), h^ = 1 + 0{gs) (3.10) 
9 o 

Meff = M + ^ a^nigsNfrigsMr, Nf = Nf + hmnigsNfTigsMY 

m>n m>n 

with amn, bmn could in principle be functions of the CY coordinates. Note that we 
have made m > n m the above expansions because the precise limits for which our 
series would be valid are: 

(^gs, gsNf, glMNf, 0, {g^N, g^M) ^ oo (3.11) 

These limits of the variables (which we will concentrate on from now on), bring us 
closer to the Ouyang solution with very little squashing of the two spheres. This also 
means that the warp factor h in (^.4|) can be written as: 



^4 



^ 3(7,M2g, f 3(7.iVf/ l\ g^Nf (.6^.62 



(3.12) 



with (iA^5 in ( p.Sp can be approximated by the angular part of the conifold metric 
(^.3|). We will however modify this further soon to get the result for large r. 
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Question now is what choices of gi are we allowed to take to insert the black hole 
in this geometry? It turns out that to solve EOMs the giS have to be functions of 
(r, 6^1, 62). Our ansatze therefore would be: 

^7i(r, ^2) = 1 - ^ + 0{glMNf), g^ir, = 1 - ^ + 0{glMNf) (3.13) 

where is the horizon, and the {61,62) dependences come from the 0{glMNf) 
corrections. We also expect a = a{rh) + 0{glMNf) as the full resolution parameter. 
Therefore in order to extract temperature from the geometry, we look at the metric 
in ( ^.41 ) in the near horizon limit r ^ and take a five dimensional slice obtained 
by setting 6i = 7i,(f)i = 0,'ip = 0. To be exact, we really need to start from the 
ten dimensional supergravity action and then integrate out the internal directions 
to obtain a five dimensional effective action. Minimization of that five dimensional 
effective action will give the five dimensional effective metric. Here we approximate 
this effective metric by taking the slice 

6, = 62 = 71, (pi = 0, 7/^ = (3.14) 

Such a choice can be justified by observing that the flavor D7 branes are all along 
this slice (see details below) 

Now, looking at the r, t direction of the metric and by change of variable, under 
the assumption that gi{r, tc, tc) ^ g2{r, vr, vr) = g{r), we can define as: 

so that the near horizon limit of five dimensional effective metric takes the following 
Rindler form: 

ds^ = -p'-f^dtl + dp^ (3.16) 
4:h{rh)g{rc) 

where prime denotes differentiation with respect to r and we only wrote the r, t 
part of the metric in terms of new variable p and tc = a/ g{rc)t. The reason behind 
rescaling time at fixed is that with this time coordinate tc, the five dimensional 
metric induces a four dimensional Minkowski metric at every r^. 



^■^Note also that all the analysis presented in the appendices are done without taking such a slice. 
What we find that taking a slice doesn't give us results very different from the full analysis. It only 



makes our calculations a little easier to handle. As an example we can quote ( 3.173 ) or ( 3.205 ) 
where the entropy is computed without taking the slice. We can see that the result differs from the 
one with slice by an 0{gsNf) term, so is very small. Similar statements can be made for the case 
with dilaton also. The dilaton profile doesn't vary too much throughout our regime of interest (i.e 
small r). So taking a constant dilaton is meaningful here. The full analysis is underway where we 
plan to take the effects of dilaton and curvatures carefully. 
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Now the temperature observed by the field theory with time coordinate tc can 
be extracted by writing the metric in (|3.16|) in the following form 



ds^ = -An^T^p^ dtl + dp^ (3.17) 
Thus comparing ( p.l6| ) and ( |3.17| ), we obtain the temperature Tc as: 



= f'^''"^ (3.18) 



J.4 

In the limit where we have gi = §2 = 9 = 1 ~ we can easily compute the 



corresponding temperature using the above formula ( p.l8 ). This is given by: 



Te = ^ + E cmn/-^^ = T, + 0(l/r.) (3.19) 



^ m,n,p 



where L is defined earlier, Cmnp is in general functions of {gs, M, N, Nj), and > 
^deconf (whcrc Tdoconf IS the deconfiuemeut temperature) is the temperature at — oo 
i.e 

Tb = Tboundary = ^ ^ ^ ^ = F(n) = T (3.20) 

where F{Tb) is the inverse transform of the above expression once we know the black 
hole factor g{rh) as well as the warp factor h{rh) to all orders in ggNf, gsM. The 
corresponding temperature function T will then be the characteristic temperature 
of the boundary theory which is fixed once and for all. This temperature function is 
a scale unique to our cascading theory and will be assumed to be greater than the 
deconfinement temperature henceforth. 

For the attentive readers, something about equations (|3.18|) and ( p.l9| ) may strike 
as odd. There seems to be an (9(l/rc) dependences. What does it mean for a physical 
variable in our theory to have an 0{l/rc) i.e the cutoff dependence? Shouldn't we 
make Tc — )■ oo to get a result that is independent of the cutoff? In other words does 
it make sense to have an explicit cutoff dependence in the physical variables? 

Of course in the standard Wilsonian renormalisation we do get results that de- 
pend on the explicit UV cutoff. In the limit where we make Tc — > oo we reproduce 
exactly the right temperature for the boundary theory. The fact that the tempera- 
ture does not have the cutoff dependence going as positive powers of should be 
pleasing: this has to do with the holographic renormalisability of the theory (that 
we elaborate in details later). However the fact that there is some cutoff dependence 
should signal something quite different from what we have in the AdS/QCD case. 

To see what is different here, let us pause for a while and ask what would the 
presence of a black hole signify in the dual gauge theory. In the standard AdS/QCD 
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case we only needed to look at the boundary of the AdS Black-Hole (AdS-BH) 
geometry to study the properties of the dual thermal QCD. Here the situation differs 
crucially because we can study the weakly coupled gravity description not only at 
r — )■ oo but also at any arbitrary r = r^. What does it mean to study the theory at 
r = Tc and not r — )■ oo? Of course, as we commented before, we can always make 
Tc — !■ oo to get the far UV results, but cutting off the theory at r = means we are 
putting a Wilsonian cutoff in the theory. Such a procedure would make sense if and 
only if we can carefully describe the degrees of freedom atr = r^. Furthermore — and 
this is one of the most crucial point — since the theory with flavor is holographically 
renormalisable (as we will show soon) all the cutoff dependences will come as (9(l/rc). 
Such a state of affairs lead to two interesting conclusions, one of which is obvious and 
the other not so obvious. The obvious point is that as long as we put our cutoff at high 
energy, the low energy dynamics remain completely unaffected by our choice of the 
cutoff. The not-so-obvious point is that the shear-viscosity etc that we will discuss 
in detail soon for geometries that are cut-off aX r = will in fact be independent 
of the cut-off once we study them from boundary point of view! This conclusion is 
surprising because we will not be making Tc — )■ oo to study the boundary theory, 
rather we will add non-trivial UV "cap" to the geometry from r = to r = oo. 

So the point that we want to emphasise here is that once we introduce a cutoff 
at r = Tc we are in principle introducing non-trivial high energy degrees of freedom 
that will in general take us away from the usual cascading dynamics of the parent 
theory! Adding such non-trivial degrees of freedom at r = is equivalent to adding 
a UV cap to the geometry. Thus it all depends on what we do at Vc i.e which 
boundary conditions we choose there. Therefore they are not universal properties 
of the cascading theory^^. Universal properties arise only as — t- oo i.e at the 
boundary^^. For example we may think of attaching A/" = 2, 4 degrees of freedom at 
r = Tc (somewhat along the lines of l^^) alongwith the remnants of the D7 brane 
degrees of freedom to have the full F-theory picture (in fact our UV completions 



""^^We thank Ofer Aharony for emphasising this point, and clarifying other details about the cutoff 
dependences. 

"'^^Alternatively this means that we are not adding any UV cap because the geometry can be 
thought of as being "cut-off" at Tc = oo. One issue here is the connection to the work of |38| . 
The question is can we write a boundary theory sX r = Vc itself instead of going to the actual 
r = CO boundary? As shown in |Q this is possible in AdS case because the theory on the surface 
of a ball in AdS space doesn't have to be a local quantum field theory. The non-local behavior in 
such a "boundary" theory is completely captured by the Wilsonian effective action at the so called 
boundary. As discussed in the text earlier, this is tricky in the Klebanov-Strassler model precisely 
because there is no unique strongly coupled gauge theory here. Again, as we saw in our RG flow 
pictures, when the dual gravity theory is weakly coupled we have a smooth RG flow in the gauge 
theory side, but the theory at any given scale can be given by infinite number of representative 
theories none of which completely capture the full dynamics. Thus it makes more sense to define 
the theories at r = oo boundary only and not at any generic r — Tc- 
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should always have 24 — Nf seven branes attached to the UV caps)^^. Clearly adding 
such degrees of freedom will give us a new theory that differs from the UV degrees of 
freedom of the parent cascading theory. In particular if is finite then this procedure 
can give us an almost free UV theory with confining IR dynamics. Also since at any 
given scale there are an infinite number of possible gauge theory descriptions available 
with different UV degrees of freedom, we can specify a particular set of degrees of 
freedom at r = to define the UV of this theory. Once this is specified the RG 
flow will take us to the IR. A sketch of the situation is depicted in figure 7 below 
where the gravity description captures the smooth RG flow. We can view the dotted 



PARENT CASCADING THEORY 




Figure 7: An oversimplified picture of the gravity description of the various gauge theories 
involved. The onion-ring model is only a caricature of a much more involved picture that we 
gave earlier. The dotted lines emanating from various points (i.e various gauge theories) 
specify the distinct UV completions of these gauge theories. Once we specify the UV 
degrees of freedom at various r = the RG flows of each theories should be assumed to 
bring us to r = points on the smooth RG flow of the original cascading theory. Note 
that at any point on the gauge theory side (i.e the vertical axis) there are infinite number 
of possible gauge-theory descriptions available (the horizontal axis). 

lines coming out from any r = Tc and going to infinity (in the gauge theory side) 
as describing the UV degrees of freedom for a theory r = r^. Thus in principle 
there is a possibility of defining infinite number of theories here by cutting off the 
geometry at various r = and then "filling" the r = r^. to r = oo spaces by UV 
caps that specify the UV degrees of freedom of these gauge theories (of course we 

^^The fact that there are infinite F-theory backgrounds gluing to our IR solutions will be discussed 
a little later. 
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expect small changes in the IR geometries every time we attach UV caps, although 
the far IR geometries wouldn't change). Under RG flows these UV theories (that are 
obviously quite different from the parent cascading theory) meet the RG flow of the 
cascading theory at various r = (where of course the degrees of freedom match 
but not the whole set of marginal, relevant and irrelevant operators). The effective 
degrees of freedom at r = are: 

iVeff = iV + log r, + ^\ „ ^ logV, 3.21 

Therefore Vc dependences in the physical result for our case should be replaced by 
the effective number of degrees of freedom at that scale, i.e by: 

TT 




SgsNf 



(3.22) 

which in turn means that the right UV degrees of freedom should flow to this value at 
r = Tc^'^. We have also defined a = ^-^V? ( ^ ^ ~^ ~ 1 I (3 — — 



3gsM \ Nf J 3gsMy/M'2-4NNf 

in the limit where N^s is small compared to the original colors and bi-fundamental 
flavors in the theory i.e > M > N^fi > Nf precisely. In general then^^ 

Tc = e^^ff (3.23) 

with A/'efr(Ac) being the effective degrees of freedom at a given scale A^,. Thus to 
summarise: 

• Making the cutoff to infinity i.e Tc — > oo gives us the property of the original 
cascading theory with a particular UV behavior and a smooth RG flow to IR. 

• Defining the theory with a cutoff at r = means that we specify the UV degrees 
of freedom of a particular gauge theory at that scale or equivalently introduce a UV 
cap in the geometry so that we are always defining the theory at the boundary^^. 
This theory also has a smooth RG flow but the UV behavior (shown by dotted lines) 



^''Even more interestingly, the final results for any physical quantities from these theories will 
only depend on the boundary degrees of freedom and not on Afcff- The AfeS degrees of freedom are 
only intermediate. However there are special cases where the boundary degrees of freedom could 
be close to MeS- We will discuss them later. 

"'^^There is another flavor-dependent degrees of freedom that one could define in this background. 
We will define this via rc(a) = ^-^'^al'^-'^io:}] ^ where the variables are defined in ( 3.35| ) and (3.36). 



^^It should be clear that unless we specify clearly the UV degrees of freedom, we cannot pin- 
point the gauge theory there. The parent cascading theory allows an infinite number of gauge 
theory description at any given scale. At the point where we have smooth RG flow, none of these 
descriptions capture the full picture there. Thus to emphasise again, the gauge theory or the UV 
behavior that we want to specify at that scale is distinct from the UV of the cascading theory as it 
could have finite (but large) number of degrees of freedom. 
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is quite different from the parent cascading theory. Only at certain IR scale (i.e 
r = Tc from the gravity side) do these theories match (wrt the degrees of freedom 
and certain set of marginal and relevant operators). 

Once this is settled, let us now introduce fundamental matters in our geome- 
try. As we discussed above, the seven branes introduce fundamental matter with Nf 
flavors. Since the above background is a deformation of both KS and Ouyang back- 
grounds, we will refer the resulting geometry as Ouyang-Klebanov-Strassler (OKS) 
background. With the black hole, this will henceforth be called as OKS-BH back- 
ground. 

In the OKS-BH background, to order 0{gsNf) and small r, let us start by 
considering the warp factor choice (|3.71) but keeping gi{r) as ( |3.13|) instead of 1. The 



20. 



background RR three and five-form fluxes can be succinctly written as' 

= dr A A (ci d9i + C2 d62) + dr A (C3 sin 6*1 dOi A d(j)i — C4 sin 6*2 d92 A d(j)2) 

ci sin Q2 d<b2 — 7: C2 sin d\ d(\)\ ) A dQ\ A d62 , 

2r 2 J 

F3 = dr A A (ci sin 61 d(j)i + C2 sin 62 d(f)2) 

9s 

H A (cs sin 9i dOi A d(j)i — Cg sin 6*2 d62 A c/02) 

9s 

1 /r t"^ -\- Qqj^ \ 

sin 61 sin 92\-C2 dOi Ci d92 \ Ad(j)iA d4>2 ■ (3.24) 

9s \2 2r y 

where is closed and F3 = F3 — CqH-^, Cq being the ten dimensional axion. The 
derivations of the coefficients appearing in (|3.24|) are rather involved, and their de- 
pendences on the resolution factor etc. will be described in section 3.3^^. For the 
present purpose, let us just quote the results: 

ci = 9sMNf ,^ _ _ 5g^2^2|Qg^ ^ aV^log(r^ + 9a^)) cot — 



^°As we mentioned earlier, at large r we expect all Ci to be finite. This means that Ci should 
become functions of inverse powers of r at large r. This is expected from F-theory considerations 
and would help us remove the Landau pole. 

^^It is instructive to note that the background EOMs cannot be trivially worked out by solving 
SUGRA EOMs with fluxes and seven branes sources. This is because, even if we know the energy 
momentum tensors for the fluxes, the energy momentum tensors for Nf coincident seven branes 
are not known in the literature! The difficulty lies in finding a non-abelian Born-Infeld action for 
Nf seven branes on a curved background. As far as we are aware of, this problem has remained 
unsolved till now. So in the absence of such direct approach, we use an alternative method to derive 
the EOMs. This method uses the ISD (imaginary self-duality) properties of the background fluxes 
and fields. Details on this have already appeared in [Q|2^, so we would refer the readers there 
for a complete analysis. Our present analysis is however more involved than [Q|2^ because we 
have a black hole and no supersymmetry. We do however find that even for this scenario, one could 
find consistent solutions to EOMs using similar arguments. For example see equations ( 3.7l| ) to 



equation (3.84) for more detailed derivations. 
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C2 
C3 



2,glMNf .2 ^ 2i / 2 ^ 2X\ ^2 

= , ^ ^ (r^ - 9a^ log + Qa^ ) cot ^ 



(3.25) 



+ 



giMNj 



36a^ - 36r^ log a + 34r^ log r 



-(lOr^ + 81a^) log(r^ + 9a^) + 12r^ log sin sin ^ 



2 2 

18a^ - 36(r^ + Ga^) loga + (34r^ + 360^) logr 



+(10r^ + 63a^) log(r^ + 9a^) + (12r^ + 72a^) log sin ^ sin ^ 



C5 = + 



C6 = + 



247r(r2 + Ga^) 



18a^ - 36(r2 + 6a^) loga + 8(2r^ - Qa^) logr 

+ (10r2 + 63a2) log(r^ + Qa^)] 



247rr2 



36a^ - 36r^ log a + 16r^ log r + (lOr^ + 81a^) log(r^ + 9a^)] 



with K = ''otc 2 ■ AH the above coefficients have further corrections that we will 
discuss later. Finally, this allows us to write the NS 2-form potential: 



9 9 
B2 = ( 61 (r) cot ^ d9i + 62 (r) cot ^ (i^2 ) A 



(3.26) 



+ 



?>glMN^^ 



f 



in 
127rr2 



9 9 

(1 + log(r^ + 9a^)) log I sin — sin — J + 63 (r) 



sin 9i d9i A d(j)i 



(-36a^ + 9r^ + 16r^ logr + r^ log(r^ + 9a^)) log ^sin ^ sin ^ + 64( 



X sin 92 d92 A (i02 
with the r-dependent functions 



^1 (0 = ^^f^i 9x (18q' + (16r' - 72a2) logr + (r^ + 9a') log(r2 + 9a^)) 
z47r(r^ + ba^j 



and 63 (r) and 64 (r) are given by the first order differential equations 



Kir) 



Sg^Mr , gsMNf 



+ 



r2 + 9a2 87rr(r2 + 9a2) 



- 36a' - 36a' log a + Mr'^ log r 



+(10r' + 81a') log(r' + 9a' 



3gsM{r^ + 6a') g^MNj 



18a'-36(r' + 6a')loga 



(3.27) 



(3.28) 



(3.29) 



+(34r' + 36a') logr + (lOr' + 63a') log(r' + 9a') 
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where M denote the remnant of the number of fractional three branes (in the gauge 
theory side) and Nf denote the number of flavors or seven branes in the dual gravity 
side. Therefore once we know B2 and the string coupling e~* then it is easy to 



determine the two couphngs at the UV of our dual gauge theory (see also P2| ): 



9l 



TT 
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(3.30) 



The string coupling can be determined easily from the monodromy around the seven 
brane to be: 



gs Sir 
which immediately gives us: 
d 



-f = 1 _ ^ log + 9aV^) - ^\og 



sm — sm — 
2 2 



(3.31) 



Slog A 
d 

dlog A 



47r2 47r2 
+ 



47r2 



gl 

If 



Nf /^6r6 + 36aV 
, ^gsNf 



3M 1 



An 



(3.32) 



which to the leading order is consistent with the Shifman-Vainshtein /3-function [^^^. 
The behavior at subleading order will tell us how the color changes as we cascade 
down from UV to IR. For example, in the presence of Nf flavors the SU {N + M) 
gauge group has 2A^ + Nf effective flavors. Under RG flow, Seiberg duality will tell 
us that the weakly coupled gauge group will become SU{N — M + Nf) x SU{N). 
This also means that the cascade will slow down quite a bit as we approach the IR 
and therefore the end point of the cascade could either be a conformal theory or a 
confining theory. As pointed out also in |^ if in the end of the cascade A^ decreases 
to zero with finite M left over, we would have SU (Af ) SYM with Nf flavors in the IR. 
Extending this to the centre of the RG flxed points surface will allow us to analyse 
this using weakly coupled supergravity. This is of course the theory we are aiming 
for. Finally, the flve-form flux is as usual given by 



{l + *w){dh-' Ad'^x) 



(3.33) 



with h being the warp factor described above. 

Before we end this section we want to point out few subtleties about the back- 
ground. First as we mentioned before, to maintain Gauss law constraint, we need to 



^^Note that the LHS of the equations ( 3.32 ) involve only gauge theory variables whereas the 
RHS of the equations involve gravity variables. In particular A should be identified with the radial 



coordinate r of the geometry. Thus the equations ( 5.32 ) capture the essence of gauge/gravity duality 
here. 
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embed our model in the full F-theory setup. This means that the background 
configuration that we presented above should be understood as an F-theory on a 
four-fold where all but Nf of the seven branes have been moved to infinity. This 
means that Nf < 24, and the four- fold is a non-trivial torus fibration over a resolved 
conifold base. Such a four-fold has already been constructed in p3, E2[ and so we can 



direct the readers to those papers for details. What is interesting here is that due to 
F-theory embedding, and with modified three-form fiuxes, we expect the background 
to not have any naked singularities or Landau poles that are often associated with 
these backgrounds (see for example [^, |0 for some details). To 0{gsNf) our result 



that we presented may indicate the presence of log r type singularities. Additionally 
at: 

f An \ 

r ^ exp - — — , for a << gsNf (3.34) 



we might think that the dilaton is blowing up on the given slice ( p.l4| ) leading to some 



kind of naked singularity. However note that we are seeing such behavior because 
we have evaluated the background locally near the seven branes, and upto 0{gsNf). 
Clearly the metric, dilaton and the fiuxes have to have a good behavior at infinity 
to be a F-theory solution. One way to show this would be to rewrite the warp factor 
( p.l2|) completely in terms of power series in r in the following way: 



h = + - 4^ + ^ (3.35) 

where ej,r(Q) etc are defined as: 

ei = — - H log sm — sm — , eo 



gsM FlNf 



TT \ N 



na) =r e(i) = -, 6(2) = -, e(3) = 6(4) = - 

r(±„) = ri^^('^), = L(2) = i:(4) = i:^ L^^) = -2L^ (3.36) 

which makes sense because we can make to be very small. The angles 6i take fixed 
values on the given slice (|3.14j) . Note that the choice of doesn't require us to have 
gsNf small (although we consider it here). In fact we can have all {N,M,Nf) large 
but ei small. A simple way to achieve this would be to have the following scaling 
behaviors of {gs, N, M,Nf): 

gs ^ e", M ^ e"^, Nf e~^ ^ e"^ (3.37) 

where e — )■ is the tunable parameter. Therefore all we require to achieve that is to 
allow: 

a + 7>2/3 + /t, a > K, 7>a (3.38) 
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where the last inequahty can keep QsNf smalL Thus QsN, QgM are very large, but 
Qs, ^^^,gsNj: are all very small to justify our expansions (and the choice of super- 
gravity background) 



The warp factor ( 3.35 ) has a good behavior at infinity and reproduces the 
0{gsNf) result locally. Our conjecture then would be the complete form of the 
warp factor at large r will be given by sum over a as in ( ^.36] ) but now a can take 
values 1 < a < oo. We will use this conjecture to justify the holographic renormal- 
isability of our boundary theory. We will also discuss the case when e(Q,) becomes an 
integer at large r. 

Similarly one would also expect the dilaton to behave in an identical way. Asymp- 
totically the dilaton should go to a constant. So on the given slice (|3.14| ), and near 
one of the seven brane, we expect: 



2gs 



1 3eaa^ 

— h constant 



(3.39) 



where the constant could be determined from the full F-theory picture. Somewhat 
similar discussion has been given in [^. Our conjecture for the large r behavior 
stems from the finiteness of F-theory. In fact this also gives us an argument to 
realise the possibility of infinite F-theory backgrounds gluing to our IR solutions. In 
the presence of 24 — Nf seven branes there are infinite possible ways by which we 
can adjust the positions of the seven branes. For every possible configurations of 
seven branes there would be non-trivial background axio-dilaton r, realised from the 
degree eight and degree twelve polynomials / and g respectively that appear in the 
Weierstrass equation (see [BOl), via: 



where 2; is a complex coordinate orthogonal to the seven branes, and 7^ bj in 
general. Since ai,bj take continuous values, there are infinite possible configurations 
of T here (modulo SL{2,Z) transformations). For certain choices of ai,bj, j{T) is 
a constant implying that there are possible configurations of seven branes that give 
rise to zero axio-dilaton (see the second and the third references of 0]). Thus 
for a generic configuration of seven branes we expect axio-dilaton r to behave as 
= 12i '^■i where {Cj} take a particular set of values for a given configuration 
of seven branes. Plugging these values of axio-dilaton in sugra equations of motion 



^■^For example we can have gs going to zero as gs e^^^ and {N,M,Nf ) going to infinities as 
(e~^, e~^, e~^) respectively. This means {gsN, gsM) go to infinities as (e~^^/^, e~^/^) respectively, 
and {gsNf, g^MNf, gsM^/N) go to zero as (e'^^^, e, e^/^) respectively. This is one limit where we can 
have well defined UV completed gauge theories. Note however that for the kind of background that 
we have been studying one cannot make Nf large because of the underlying F-theory constraints 
pO| . Since we only require gsNf small, large or small Nf choices do not change any of our results. 
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alongwith a similar configuration of fluxes, we can easily argue the existence of infinite 
configurations of warp factors of the form (|3.35| ), justifying the IR changes from UV 
caps. A particular configuration of F-theory background will glue to our IR solution 
to give us the required UV completion^^. In fact the configuration with constant 
coupling (like the last two references of [0) will give rise to AdS completions of our 
IR backgrounds! More details will be presented elsewhere. 

One caveat is that a full analysis incorporating non-perturbative effects still 
needs to be performed to justify the whole scenario. However we expect this to 
be very involved because at large r we have to consider not only the effects of all 
the (p, q) seven branes (as discussed above) but also the back-reactions from fluxes 
etc.^^. We will address this in the sequel [^. Therefore with this assumption, all 
the log r dependences of the fluxes should also be replaced by inverse powers of r 
at large r. We will however, for the purpose of concrete calculations, only work to 
0{gsNf) locally for many representative examples unless mentioned otherwise. The 
singularities appearing in these examples at large r should then be considered as an 
artifact of the order at which we do the analysis. This will also be clear from the 
holographic renormalisabilities of these theories^^. 

3.2 Quark mass and drag coefRcient 

We now embed a string in OKS-BH background with one end attached to one of the 
D7 branes and other end going into the black hole (see figures 8 and 9 below). 
If X*(cr, r) is a map from world sheet coordinates a, t to 10 dimensional space time, 
then string action or fundamental string Born Infeld action is (see for example pl|): 

^string = ToJ rfadr [yQ<^7+^^W) + l^^'^-b + J{4>) 



^^More precisely, the non-trivial geometry of the UV will capture the effects of the operators 
defined at the cut-off; and the F-theory seven-branes will capture the effects of the ultra massive 
fundamental quarks that help us remove the Landau poles. 

^^We thank Ingo Kirsch and Diana Vaman for discussion on this issue, and for pointing out the 
reference pof . 

^^One might wonder whether the procedure of holographic renormalisability should work for the 
full F-theory picture. It is easy to see why this should proceed without any complication: F-theory 
background simply gives us the non-perturbative completion of a given IIB background. These 
non-perturbative objects are the seven branes and they contribute to the bulk lagrangian as some 
matter multiplets. As long as the effects of the matter multiplets are not too large (as discussed 



after ( 3.97 ) and before ( 3.98 )) the procedure of holographic renormalisability should proceed as 
in the usual supergravity case. However if we incorporate these multiplets they still don't change 
anything because the procedure of holographic renormalisability requires derivative interactions that 
come exclusively from \/—GR term of the lagrangian. All the fluxes etc contribute to polynomial 
interactions, as discussed around ( p. 117 ). Of course the effective values of the metric fluctuations 



are affected by the background seven branes. 
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D3 Branes 

D5 Branes 

Figure 8: Brane construction that we used to describe the zero temperature gauge theory. 
The fundamental flavors come from the wrapped D7 branes whereas the bi-fundamental 
flavors come from the wrapped D5 branes. To introduce temperature all we need to do is 
to Euclideanise and identify the time coordinate in the gauge theory. 



D 7 Branes 



Figure 9: The dual gravity picture for the high temperature gauge theory. The branes 
(except the D7 branes) are replaced by fluxes and the quark string has one end on the D7 
branes and the other end going into the black hole. The conifold is replaced by a non-trivial 
geometry that we discussed in the previous section. 



String (x) =Toj dadr y-detiU^ + 5,09^0) + l^^'B^b + J{4>) 

+ e r^r'^'"=-T„ e Fabc... + c(e^)j ^^^(x - %) (3.41) 

where J(0) is the additional coupling of the dilaton to the string world-sheet, Tq is 
the string tension, X" are the ten bosonic coordinates, B is a 32 component spinor, 
Fahc... = [dC]abc.. with Cabc. being the background RR form potentials, and f^s is 
world sheet metric, given by the standard pull-back of the spacetime metric on the 
world- sheet: 

^^[x-X'X'-X')^[ ^^L^ ^ + ^ I ^^-^^^ 

with J or 6 = 0,1, 1]^ = T,r]^ = a and Bab is the pull-back of the NS two form field. 
In the ensuing analysis we will keep both B^b as well as da4> zero. The former will 
be imposed (see discussion after (|3.75|) ). The latter case will be addressed soon. The 
interesting thing however is to do with the background RR forms. Note that the 
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RR forms always couple to the 32 component spinor. Therefore once we switch-off 
the fermionic parts in (|3.41| ), the fundamental string is completely unaffected by the 
background RR forms^''. Thus in the ensuing analysis, for the mass and drag of the 
quark, we can safely ignore the RR fields. We have also defined: 

det / = -G.jX'X'^ + {G,,X''X'^){GmX''X') 

X-.f: (3.43, 

where Gij is more generic than the background metric, and could involve the back 
reaction of the fundamental string on the geometry. The analysis is very similar to 
the AdS case discussed in so we will be brief in the following. However, since our 



background involves running couplings, the results will differ from the ones of . 

At this point it might be interesting to point out an equivalent classical calcula- 
tion to evaluate the effect of a drag force on a point charge particle moving through 
a media. One can consider the point charge to be the endpoint of the fundamental 
string on the D7 brane. The mass of the point charge is given by the length of the 
string / as m = TqI where Tq is the tension of the string. In the presence of a constant 
electric flux J-'oi one can show that the total drag force Fi on the point charge in the 
dual gauge theory is given by: 

= -gJ-Q. (e-*/^ - 1) (3.44) 

where r is the so-called collision time. This can be easily determined from the 
mean free path of the collision. Viewing the end points of the open strings on the 
D7 brane as a gas of charged particles, we can study the number of collisions at 
any given point in the gas of particles at a particular temperature, and from there 
evaluate the collision length between them. The number of collisions Z happening 
per second per unit volume at a point r (with r not to be confused with the radial 
coordinate r in the gravity picture) in the gas can be written as: 

Z = j d^pi jd^P2(y\-Vi - V2U(r,Vi,t) ^(r,V2,t) (3.45) 

where ?, v are vectors and g(r^,\,t) is the distribution function. The quantity a 
is the collision cross-section, and is a constant for our purpose. The distribution 
function g(r,\,t) can be anything generic, and only becomes Maxwell-Boltzmann 
when the gas attains equilibrium. If there are n strings per unit volume of the gas, 
one can show that the mean free path Im in the dual theory is given by the following 
expression: 

71 

Im = (3.46) 



^''This is of course the familiar statement that the RR fields do not couple in a simple way to the 
fundamental string. 
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where Vp is the most probable velocity which differs from the average velocity (v) by 
a numerical constant. 

Finally when the system attains equilibrium the distribution function, as dis- 
cussed above, becomes MB. In that case there are no r,t dependences in g{r,\,t), 
and we have 

where m is the mass of a single particle in the gas, and k the Boltzmann constant. 
Using this distribution it is straightforward to show that the mean free path Im is 
independent of temperature and is given by: 

^r. = ./I ■ — (3.48) 
V 8 na 

from which the collision time r can be determined to be r = Plugging this in 
( p. 441 ) gives us the drag force on a particle in the dual theory. 



The above analysis is clean and simple, but doesn't completely give us the full 
answer. This is because, the above analysis ignores higher order quantum corrections 
on the collision process. Such corrections are in fact captured by the classical back 
reactions of the underlying geometry in the gravity side^^. The important fact that 
the D7 brane is wrapped on a curved manifold, changes much of the above analysis. 
To see the effect of the background geometry on the process, let us evaluate using 



the technique of 



To begin, we need the embedding of the D7 brane in our set-up. This has 
already appeared above as ( p.8|) . We can use the embedding equation to determine 
tq, the distance upto which the D7 ends in the throat. The embedding equation (|3.8| ) 
implies: 



1/^ 



2 



sm^ ^ sin^ f 



(3.49) 



from where, by minimising r with respect to the angular coordinates 6i, we obtain 
9i = 02 = TV and from here one can see that rg = l/ip''^ = 1.^^ 

As we said before, a fundamental quark will be a string starting from tq on the 
D7 brane to the horizon rh of the black hole. For simplicity of the calculation, we 



^^This is of course pure supergravity analysis valid for gauge theory at very strong coupling. 
There would be stringy corrections over and above this result once we go closer to the boundary of 
the RG surface. This will additionally introduce 0{gs) corrections to the existing result. For the 
time being we will ignore these effects as we will always be in the regime shown in figure 6. 

^^This is done to make = f dimensionless. Note that this is another scale in our problem, 
which we identify to 1. As alluded to before, this will help us to make most of our QCD variables 
dimensionless. 
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will then restrict to the case when 



X° = t, X^ = r, X^ = x{a,T), = (A; = 2, 3, 5, 6, 7) 

(X^ X^) = [Oi, 62) = 71, e = = (3.50) 

and we choose parametization r = t, a = r also known as the static gauge. Thus 
we are only considering the case when the string extends in the r direction, does 
not interact with the RR fields, and moves in the x direction of our manifold. More 
general string profile, while being computationally challenging, does not introduce 
any new physics and hence our simplification is a reasonable one. 

Before moving further, let us consider two points. First is the effect of the black 
hole on the shape of the D7 brane. We expect due to gravitational effects the D7 
brane will sag towards the black hole and eventually the string would come very close 
to the horizon. In fact, as is well known, putting a point charge on the D-brane tends 
to create a long thin tube on the D-brane that in general extends to infinity. The 
end point of the string being a source of point charge should show similar effects (see 
3 for a discussion of a somewhat similar scenario) In this paper we will ignore 



this effect altogether and discuss it in details in the sequel to this paper 



The second point is to see how the background varying dilaton effects the string. 
Near the local region around the string, we expect the following behavior of the 
dilaton: 



log Qs + , log r + — (3.51) 



where we could insert /i to make r dimensionless. In the limit where the resolution 
parameter a goes to zero, it is easy to see that near r/j and tq = 1, the dilaton behaves 
respectively as (inserting the scale /i): 

(p log gs H log 



An |/i|2/3 

^ log - ^ log |/i| (3.52) 
An 

which means that near the seven brane the dilaton behavior is almost a constant 
although there is a log T dependence of the dilaton near the horizon, where T is 
defined in (|3.20|) . This behavior is quite different from the AdS case where there is no 



profile of the dilaton. As we will show later, there will be additional log T behavior 
coming from the warp factor also, which for small QsNf will have dependence on 
Thus the log T dilaton dependence can only contribute to order 0{gsNf) if 
we are close to the horizon. Away from the horizon and near the seven brane, the 



■^'^Even in the supersymnietric case, putting a point charge on a D-brane tends to create a long 
tube that extends to infinity. One can then view an open string to he at the end of the thin tube. 
This effect is somewhat similar to the one discussed in [p6| . 
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dilaton is approximately constant. Therefore to simplify our ensuing analysis, we 
will take the dilaton to be a constant and only consider the detail implication in 
the sequel. In the following analysis we want to convince the reader that even with 
these simplifications: Bat = J{4>) = in ( |3.41|) , our system will have interesting new 



physics. For more details on the dilaton behavior see l^Tj . 

Now introducing a fundamental string in the geometry will change the IR geom- 
etry. We expect the back reaction to be small, i.e we can specify the back reaction 
via the perturbed metric Gij = Qij + nlij where Qij is the original metric and nlij is 
the back reaction. With our choice of parametization and string profile, it is easy to 
verify: 

- det / = ^ + TT^i^x 2 - g,{r)-'x' + + 0{k') (3.53) 
g2{r) h{r, 7t,7t) 

where the warp factor h{r, 6i, 62) = h{r, vr, vr) with (flmn; &mn) ~ in (|3.12|) henceforth 
(unless mentioned otherwise); and the back reaction of the string appear as 0{k,) 
effect. 



With this, the rest of the analysis is a straightforward extension of |^^. The 
Euler-Lagrangian equation for = x{t,r) derived from the action ( ^.41| ) and the 
associated canonical momenta are: 

Id/ X \ d / gix' \ 
^rft Vy^diFTv' ^ d^Vh^^Sj) ~ 
{X ■ X'){Xi)' - {X'f{X^) 



v/=dir/ 



n° — — ToGjj- 

n, - -ToG., -^=== (3.54) 

If we consider a static string configuration, i.e. x(cr, r) = b = constant, then energy 
can be interpreted as the thermal mass of the quark in the dual gauge theory. Using 
the static solution in ( p.54| ), we obtain the mass using E = — J da 11°, as 

m{T) = To(ro - rj,) = Tq -T)=n{l- T) (3.55) 

Now recounting our earlier analysis of the drag force (|3.44|) we see that the velocity 
of a particle in the gas (given by the end point of the string) will eventually approach 
a constant velocity v as 

V(t) = V (1 - e"*/") (3.56) 

where Vi = with m derived at a constant temperature. Thus once the string (or 
the fundamental quark in the dual gauge theory) moves with a constant velocity, we 
need to apply force constantly to keep it in that state. This gives us the drag force. 
To obtain the drag coefficient z/, we consider strings moving with constant velocity: 

x{t,r) = x{r) + vt (3.57) 
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where v = vi from (|3.56|) , as after large enough time V{t) — )■ f and the string moves 
at a uniform speed. Then from ( p.54|) , noting that / is independent of time, we can 
solve the equation of motion to get: 



X 



12 



2„,2 



9i 



ur2 2 (3-58) 

where C is a constant of integration that can be determined by demanding that 
— det / is always positive. Using the value of x'"^ from ( p. 5^ ) we can give an explicit 
expression for the determinant of / as: 

,2 

(3.59) 



det/ = ^ 



9i - ^ 



92 9i - hC^v^ 

For —det / to remain positive for all r, we need both numerator and denominator to 



change sign at same value of r. This is the same argument as in [44|. The numerator 
changes sign at^^ 



+ 0{gsNf,gsM) 



Requiring that denominator also change sign at that value fixes C to be: 

1 



C 



r2T-2 



1 + 



2nN 



log 



(1-^2)1/4 



1 + 



2tt 



log 



(1-1)2)1/4 



+ 



(3.60) 



<3.61) 



where M and Nf differs from M,Nf due to the 0{gsNf, QsM) terms in ( p. 601) . The 



first part of C is the one derived in |^ . The next part is new. Now the rate at which 
momentum is lost to the black hole is given by the momentum density at horizon 

UUr = rh) = -ToCv (3.62) 

while the force quark experiences due to friction with the plasma is ^ = —up with 
p = mvj\/\ — . To keep the quark moving at constant velocity, an external field Ei 



does work and the equivalent energy is dumped into the medium This external 
field Ei is exactly the flux T^i discussed above. Thus the rate at which a quark dumps 
energy and momentum into the thermal medium is precisely the rate at which the 
string loses energy and momentum to the black hole. Thus upto 0{gsNf, QsM) we 



have 



—Il^{r = rh) and 



ToCVT 



m 



(3.63) 



Tn 



1 + 



2ttN 



log 



r 



(l_t,2)l/4 



1 + 



3gsNf 

2-K 



log 



r 



(l__„2)l/4 



+ 



3iNote that by 0{gsNf,g^M) we will always mean O^g^N f , g^M N f , g^M"^ / N) unless mentioned 
otherwise. 
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which should now be compared with the AdS resuh |44|. In the AdS case the drag 
coefficient v is proportional to T^. For our case, when we incorporate RG flow 
in the gravity dual, we obtain O {xj \l A log T + B log^T^ correction to the drag 



coefficient computed using AdS/CFT correspondence 
3.3 Wake created by the moving quark 

In the previous section we computed the drag force on the quark. Clearly a moving 
quark should leave some disturbance in the surrounding media. This disturbance is 
called the wake of the quark. Thus, in order to compute the wake left behind by a 
fast moving quark in the Quark Gluon Plasma, we need to compute the stress tensor 
T^'^, p, g = 0, 1, 2, 3 of the entire system. Our goal therefore would be to compute^^: 

medium+quark quark 

where the ffist term is basically the energy-momentum tensor of OKS-BH background 
plus string i.e ^background +string restricted to four-dimensional space-time. Similarly 
the second term is the energy momentum tensor of the string i.e Tj^^j^^g restricted 
to four- dimensional space-time. This is similar to the analysis done in for the 



AdS case. For our case the above idea, although very simple to state, will be rather 
technical because of the underlying RG flow in the dual gauge theory side. Our 
second goal would then be to see how much we differ from the AdS results once we 
go from CFT to theories with running coupling constants. 

For a strongly coupled QGP, we will apply the gauge/gravity duality to compute 
T^'^ of QGP using the supergravity action S'totai- The supergravity action will be 
defined as a functional of the perturbation Ipq from the string on the background 
metric gpg. Making use of the duality, we expect that the Hilbert space of strongly 
coupled QCD to be mostly contained in the Hilbert space of low energy weakly 
coupled classical Supergravity i.e. the OKS-BH geometry of ( p. 41) , or alternatively, 
the full Hilbert space of QCD should be contained in the Hilbert space of string 
theory in the OKS-BH background. However there is a subtlety here. The standard 
supergravity analysis in this theory will lead to actions that blow up at the boundary 
(i.e taking r = Tc — > oo). The reason for this is rather simple to state (see also |^). 



The UV completion of cascading type theories require infinite degrees of freedom — 
much like string theories. This is of course another reason why the dual of cascading 
theories are given by string theories. Once we require infinite degrees of freedom at 
the UV, we no longer expect a finite boundary action from supergravity analysis! 
What we need is to regularise and renormalise the supergravity boundary action so 
that finite correlation functions could be extracted. This would also mean that the 
usual Witten type proposal |^ for the AdS/CFT correspondence can be re-expressed 



^^In the published version this section is sHghtly abridged by emphasising mostly on the holo- 
graphic renormalisability of the theory with fundamental flavors and very little on the actual wake. 
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in terms of the boundary variables to give us the complete picture. Therefore we can 
rewrite the ansatze proposed by Witten et al for our OKS-BH geometry to take 



the following Wilsonian form (also as mentioned in footnote 2, by QCD we will only 
mean a theory that approximates large N IR QCD, but with a strongly coupled and 
almost conformal UV): 

2^QCd[0o] = (exp / 000) = ^total[0o] 

= exp(S'total[0o] + 5'gH + 5'countcrtemi) (3.65) 

where is Minkowski manifold, Sgh is the Gibbons- Hawking boundary term [^] , 
00 should be understood as a fluctuation over a given configuration of field, and 
S'counterterm IS the couutcr-term action added to renormalise the action. Observe that 
in the usual AdS / CFT case we consider the action at the boundary to map it directly 
to the dual gauge theory side. For the OKS-BH background, as we discussed above, 
there are many possibilities of defining different gauge theories at the boundary 
depending on how we cut-off the geometry and add UV caps. Taking the radial 
coordinate as setting the energy scale A the gauge theory side would make sense at 
that scale once we define the UV degrees of freedom there. This is what we called 
A/'cff(A) in (|3.23|) . Therefore in general the action at any point r = rc in the OKS-BH 
geometry will map to the dual gauge theory with A/'eff(A) degrees of freedom at that 
energy scale. The properties of this dual gauge theory may not necessarily coincide 
with the universal properties of the parent cascading theory when both are studied 
from the boundary. However the RG flow of this theory will eventually catch up with 
the RG flow of the smooth cascading theory at that scale fixed by our choice r = Tc 
(see figure 7). For large enough Tc the above correspondence ( |3.65D should give us 
finite boundary action. Our procedure then would be to fix the boundary action for 
large rc (typically Tc — )■ oo) by adding the corresponding 5'gh and S'counterterm, and 
then extrapolate this to smaller radii'^^. Thus, for computing stress tensor of QCD, 



we have O = T^'^ and 0o = i^lpq- K follows that 



K 6lpq 



(3.66) 



■^■^Once the boundary theory is holographically renormahsed, all rc dependences would go as r~^. 
In the dual gauge theories this would mean that in addition to the universal properties inherited by 
each of these gauge theories, there would be corrections associated with the UV degrees of freedom 
that we needed to provide to define these theories dX r — rc- These are precisely the corrections 
that take us away from the parent cascading theory and give us results associated with new gauge 
theories. Needless to say, such extravagant richness of physical theories are not available to us in 
the AdS/QCD picture. Note however that even for the AdS/QCD case there might arise situations 
where we would require Af = 2 degrees of freedom to UV complete a Af = 1 configuration at IR. 
The issue of holographic renormalisation is not much affected by this because such a UV completion 
only affects the internal space (here it may change T^'^ geometry to /7i2). For our case the full 
F-theory completion of the parent cascading theory would also be UV complete. 
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where S = S'totai + <S'GH + <S'counterterm and S'totai IS the low energy Type IIB supergravity 
action in ten dimensions defined in string frame as: 



total 



SUGRA 



+ s, 



string 



s 



D7 



~l~ 'S'solitons 

(3.67) 



where Ssugra is the background supergravity action and S'soutons is the action of the 
sohtonic objects in our theory, namely the fundamental string and the D7 brane. 
Recall that the OKS-BH background is constructed by inserting the D7 brane in 
the supergravity action with an additional black hole singularity. Therefore we will 
define 



OKS-BH 



SUGRA 



+ s- 



D7 



(3.68) 



so that the background solutions that we gave in the previous subsection will corre- 
spond to S'oKs-BH- Once we introduce an additional fundamental string we expect 
some of the background values to change. The change in the metric will take the 
following form: 



lij — lij (r, X, y, z, t) 



(3.69) 



where gij is the OKS-BH metric and Uj = 0, .., 9) denote the perturbation from 
the moving string source (with /t — )■ 0). We also expect the NS two form B2 defined 
in ( p.26| ) to pick up an additional component along the (Or) direction. The three and 
five forms RR field strengths defined as: 



F5 + A F3 



(3.70) 



would change from the values defined earlier because changes. In the absence of 
the string, F3 is the three form sourced by D5 branes and F5 is the five form sourced 
by D3 branes. 

To proceed, let us first figure out the possible changes in H3, the NS three form 
field strength from ( p.24|) given earlier. To order 0{gsNf) locally the background 
value of in the absence of the fundamental string for a Klebanov-Tseytlin type 
geometry 

Oi 
2 

9sNf 



Cl 



^ can be given by the following values of Cj pO 



1 52 



C3 = C4 



Airr 
3g,M 



C2 



coi j + 0{glN]) (3.71) 



1 + 



9 log r + 2 log ( sin sin 



+ 0{giN 



-43- 



where the O^g^N^) local corrections are typically of the following form (see also 
111, 11,11,111): 



n>m;p 



(3.72) 



with ttmnp are not in general constants. For certain examples studied in [51, 52, 53, 18 



without D7 branes, amnp are functions of the radial coordinate r. With the D7 branes 
these corrections have not been computed. 

In addition to the 0{glN'j) corrections we have another set of corrections already 
to order 0{gsNf) that appear because of our choice of background ( ^.71 ). These 
0{gsNf) corrections are in general difficult to work out if in { \i.7\ j the parameter a is 
not a constant. When a is a constant then this is a small resolution of the conifold 
and changes the coefficients q of ( p.71|) in the following way: 



Aci 

Cl 
AC2 

AC3 

AC4 
C4 



6a^ 



-18a^ 



log r'^ — 2 
log r 



9a2 ^gsNftt"^ 



Oia'log a) + 0{g',N]) 

+ 0{a') + 0{glN'f) 

8 + 9 log r — ^ log a 



47^2 



1 + ^ (9 log r + 2 log sin f sin f ) 



(3.73) 



+ 0{aHog a,glNh 



3a2 



1 + ^ (3 + log sin I- sin |) 
1 + ^ (9 log r + 2 log sin | sin f ) 



+ ^(a^log a) + 0{glN}] 



Putting everything together we see that the NS three form changes a bit from what 
is given in |51|, |52|, |T8|. To order gsNf the local three form is given by: 



i/3 = QgsM 1 



9gsNf 



9sNf 



. log r + 
4tt 2tt 
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00 



01 
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ISgsMa^ 
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9gsNf 
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27r 2 2 I r \ 
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00 



91 



27T 



— A e^, de^ A cot — d92 + cot — d9i 



+ 



0{a'g,Nf) + 0{glN]) 



2 

A doj'i 



(3.74) 



where we only pointed out those 0{o?gsNf) terms that have similar forms as the 
terms of with a = 0. There are a few more 0{a?gsN f) terms described by generic 
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three form du^ in the OKS-BH geometry, that we do not write here but could be 
easily worked out. In addition to that we have O^g^Nj) terms of the form ( p.72| ) 
that we need to incorporate. The two form U2 is defined above as: 



W2 = - (sin 9i dOi A 



sin 6*2 d02 A d(p2) 



(3.75) 



Notice also that the fourth term in (|3.74|) has a relative sign difference from the second 
term of the angular forms. Similarly the third term has an unequal distribution of the 
three form on the base two spheres. In additional to that — for duo^ given exclusively 
by the internal three forms — we can always define a B2 field that lie only along the 
angular directions (6'j, (piiip). In that case the dynamics of a fundamental string will 
get influenced by the background NS field (so to avoid issues like we will take 
the slice ( |3.14| ) where the pull-back Bah in ( |3.41| ) is zero). Also in this paper we will 
only consider the case where B2 can be made to lie in the internal angular directions 
by appropriate gauge transformations. We should however remind the readers that 
this is not generic. In the presence of branes and fluxes such procedure cannot 
always be done without generating noncommutative theories on the branes. 

Let us now determine the RR three form F3 which is a combination of F3 and 
H^. For a Klebanov-Tseytlin kind of background we expect, in addition to the 
Cj defined earlier in ( |3.71|) , there would be two more c, given by: 

^gsNf 



c-, = cq = gsM 1 + 



2tt 



log r ) + Oig^N]) 



(3.76) 



As mentioned earlier, the above choices of q are still incomplete. For constant small 
resolution a we expect Aq to be given by: 
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+ 0{a^\og a) + 0{glN^f) 



(3.77) 



Combining everything together, we then get the following value of the background 
RR three form F^: 
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27gsMNf aHog r 



4tt 



01 



sin 6*1 sin 62 cot — dOi — cot — d92 A d(f)i A u(/;2 
\ 2 2 



+ 



A (iw.s 



(3.78) 



where as before we considered only the 0{a'^gsNf) terms that are proportional to the 



existing resolution free terms that appeared in ||2y, |T8]. However our form is more 
involved than the ones considered earlier even if we ignore the 0{glN'j) corrections 
of the form ( p. 72 ). The fact that the background has a resolution changes much 
of the details. Notice also the fact that these corrections cannot be absorbed as a 
renormalisation of the a = terms because the form structures have relative sign 
differences. 

In the presence of a non-trivial dilaton, these three forms combine together to 
give us G3 = F^—ie'^^H^. It is instructive to construct G3 for our background because 
most of the details can be expressed directly in terms of G^. The full expression for 
G3 is involved, but could be combined succinctly using certain one and three forms 
to take the following background value: 
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where we have used two new two-forms uj2 and 0J2 defined in terms of uj2 in the 
following way: 

_ 9 

(X'2 = (2 — Slog r)uj2 — - log r sin 62 d62 A d(j)2 

U2 = UJ2 + sin 02 d62 A (i02 (3.80) 

These forms help us to express the deformation of G3 once we take the resolution etc 
into account. The other effective M and Qs are then defined in terms of {M, Qs, Nf) 
in the following way'^'*: 

~ 9gsMNf ^ Ma^ ^ 87r _ , , 

M=^. M = -—, 9s = 9s = l2gsr 3.81 

Using these definitions we can see how the back reactions effect the three forms 
in our setup. Note that the deformations appear in G3 almost exactly like the 
undeformed forms, but there are crucial relative signs and extra (r, 6i) dependent 
factors. However the way we have written the backgrounds is not very productive 
because of many complicated terms. But there exist an alternative way to rewrite 
the above background which would tell us exactly how the black hole modifies the 
original Ouyang setup. This can be presented in the following way: 

F3 = 2MAi (1 + %^ log r ) A ^ (sin 9i d9i A d(pi - Bi sin 62 A rf^s) 
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where we see that the background is exactly of the form presented in |2D| except 



that there are asymmetry factors Ai, Bi. These asymmetry factors contain all the 



■^"^Since upto 0{gsNf) our metric remains similar to the Ouyang metric, we then expect the 

effective number of fiveform flux to change as A'cff = N + '^^f,^^ ^log + ~^^log^''^ ■ This would 

imply that under radial rescaling r — >■ g-^^/^gsM ^.^ jy^^ decreases by M — Nf units, exactly as we 
discussed earlier (see also pOt). 
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informations of the black hole etc in our backgrounds^. To order 0{gsNf) these 
asymmetry factors are given by: 
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These asymmetry factors tell us that corrections to Ouyang background |20 



come 



from 0{a^ /r"^) onwards. Thus to complete the picture all we now need are the values 
for the axion Cq and the five form F5. They are given by: 

Co = ^(^-01-02) 

Air 

F5 = —[d'^x Adh'^ + *{d^x Adh'^)] (3.84) 
9s 

with the dilaton to be taken as approximately a constant near the D7 brane and 
h is the ten dimensional warp factor discussed above. Thus combining ( p.82|) and 
( p.84|) our background can be written almost like the Ouyang background |^ with 
deviations given by ( p.83| ). 

So far we got the background without taking the back reaction of the string. 
The back reaction of the string can be computed from its energy momentum tensor. 
Using the action ( 3.41| ) we can obtain the energy- momentum tensor of the string as: 
6S.- 



rpij / X ^ "'-^String 
string V'^/ SG - 



'2X ■ X'X'X'^ - X'^X'^X^ - X'X^X 



dcrdr 



2V-det/ 




■^^One can easily see from these asymmetry factors that one of the two spheres is squashed. As we 
mentioned before, this squashing factor is of order 0{gsNf) and therefore could have a perturbative 
expansion. Note also that although the resolution factor in the metric is hidden behind the horizon 
of the black hole the effect of this shows up in the fluxes. As far as we know, these details have not 
been considered previously. 
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where all the variables have been defined earlier. 

To study the back reaction without incorporating the backreactions from other 
F-theory seven-branes, we need to analyse the geometry close to the string. This will 
be rather involved because our metric ( p.4| ) with the choices of warp factor ( p.l2| ) 
and the black hole functions ( |3.13| ) are complicated. However we can impose two 
immediate simplifications: 



h{r, 9i, 6*2) = h{r, vr, vr), gi{r) ^ g2{r) = g{r) 



(3.86) 



which is motivated from the fact that we are close to the string and the Olg^MNf) 
corrections to the black hole factor are subleading (as we saw in the previous sec- 
tion). Henceforth we will stick with these choices throughout our analysis. With this 
therefore we see that the local metric is given from ( p.4|) and ( 3.12| ) (or ( 3.35|) ) as: 

ds^ = — {l - A\og r - B logV) {-gdt^ + dx'dx^) + - — 



+B logV) + (L^ + ALHog r + BLHog\) dMl 
where A and B are defined in terms of gs, N"^^, N and Meg as: 
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(3.87) 



(3. 



with Meff and Nf are defined in ( |3.10D ; and = y/AngJ^ being the usual definition. 
For supergravity to be valid we require weak string coupling but large such that: 



gs ^ 0, 



Nf 



M 



» 1, « 1, — < 1, 



» 1, gsM » 1 (3.89) 



Using these limits (see also footnote 23 for more precise parametrisation) one can 
easily show that both A and B can be made very small, and 
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(AL^ BL"^) « 1 



(3.90) 



The above equation is very useful for us because we can recast our metric using 
( p.90|) to show how much we deform from the AdS black hole metric. Since » 
{AL"^, BL"^) we can rewrite ( p^STj) as^^ 



ds' 
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(3.91) 
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Using the parametrisation given in footnote 23, observe that A — ^ e^/^, B — >• and L"^ ^ t '^^1'^. 



Thus AL"^ 6^/4, ^ £13/4^ _^ g29/4 B _^ ^35/4 fo^. ^ ^ rpj^jg would justify the hmits 



considered above. 
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where the first hne is the AdS^ x A^5 black hole solution, and the second line is the 
deformation of the AdS^ and the internal geometries. Observe that (a) the internal 
space and the radial direction are very mildly deformed in our limit for regions 
close to the string^^, and (b) the deformation of the AdS geometry is not another 
AdS space. These conclusions also imply that we can integrate over the internal 
coordinates 0i, 02, 6*2 in ( |3.67D , to obtain the five dimensional effective action 
S^^^^i in the following way: 

-^tftai = / d^xV^le-^'f'RiG) - 2A(r) + G.^d^ ■ d^^ - m,,$^$^' + 

SttGn J L 

+T5 / d^x F A*F + scalars + ^ c„i?" + ^ 6„tr {R A R) 

n>2 n>l 

— 'S'oKS-BH ~ 'S's^'l^j (3.92) 



•-^string 



where On is the Newton's constant, T5 is the effective tension, A(r) is the cosmolog- 
ical "constant" coming from the contributions of the background fields and are 
the scalars that we get by dimensionally reducing IIB supergravity fields over the 
internal manifold. Some of these scalars come from the metric fluctuations of the 
internal five- manifold gmniv) as: 

G^^^{x,y):=G,,ix) © J2'^'{x)niM © gUy) (3.93) 

i 

where (x, y) are the five-dimensional spacetime and the internal indices respec- 
tively, Gmn{x, y) would incorporate the background geometry ( p.91| ) as well as string 
contribution (to be discussed below), Vt^ are the normalisable p-form satisfying 
J fi* A Vt^ = 6^^ . Similarly the other scalars (excluding the ones that come directly 
from the axio-dilaton and the five-form) come from the two two-form fields B in the 
following way: 

$i = j {b- (B)^ a (3.94) 

where (B) are the background B^s and Bur fields that we derived earlier. The rest 
of the scalars and the gauge fields are from the wrapped D7 brane. We have also 
included the higher derivative R^ and tr [R A R) type terms that come from the 
D7 brane back reaction and world volume Cherns-Simons terms respectively. These 
terms will have important implications that we will discuss in the next section. For 
the time being we want to point out that 0{R'^) terms have been studied recently in 



'^''This is because the coefficients of dr^ and tlie internal space are proportional to which go 
to infinity as e~^^^'^, whereas the other parts are proportional to which go to zero. Clearly 
then the dr^ and the internal space are strongly dominated by whereas the other parts can have 
small log r deviations. Additionally any UV modifications will not affect the IR geometry. 
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the context of the rf/S bound, namely the viscosity over entropy bound in |58|. It 



was found therein that the rj/S receive quantum corrections that doesn't lower the 
known ^ bound of Kovtun-Son-Starinets 0. However an 0{R^) term should lower 
the bound as discussed in , although one needs to be careful about two issues: the 
sign of the term and the relative strength of the i?^ and terms^^. 

Once we have our effective five dimensional action, we can derive linearized 
Einstein equation using this effective action. At this point we can parametrise the 
string contribution as hil^v with k — )■ 0. This means that the perturbations are 
small, which is a reasonable assumption for our case. Thus the total metric is G^p = 
Qfiu + ^lu with g^i, given by (|3.91|) . The equation of motion for Gfj^u would be: 



R,. (gap + nLf^) - I {g,u + kI,,) R {ga^ + nl^p) = T^r^ + T^^' + t(^''^)^(3.95) 

where the T^'^™^ come from the five- form fluxes F(5) (that give rise to the AdS^ part) 
and the remnant of the H^s^ Hrr and the axio-dilaton along the radial r direction 
(that give rise to the deformation of the AdS^ part). In five-dimensional space these 
fluxes would appear as one- forms fJ*'' with i = 1, ...,4. The effect of Tjfj'''^'^ will not 
be substantial if we take it as a probe in this background. 

At this point we can approach the problem in two ways. The first way is to 
assume that the fluxes contribute to the five dimensional cosmological constant A(r) 
as given in ( p. 92 ). In this set-up we have an effective five- dimensional theory ( p.92| ) 



and we put a string in this background to study the perturbation in the metric. Here 
the assumption is that the string do not back-react on the cosmological constant 
A(r). This seems to be the general approach in the literature. 

The second way is to actually consider the back reaction of the string on the 
five dimensional cosmological constant. This can be worked out if one considers the 
effects of all the background fluxes in our theory. The result of such an analysis can 
be presented in powers of k. For our case we are only interested in back reactions that 
are linear in n. To this order the equation of motion satisfied by lap is determined 
by expanding ( p.95| ) in the following way: 



K 



- Bf, - A'^J) = Tfr' (3.96) 



where A^^ is an operator whereas B'^^ and A'^^ are functions of r, the radial co- 
ordinate^^. We have been able to determine the form for the operator A|^^ for any 
generic perturbation /q,^. The resulting equations are rather long and involved; and 
we give them in the Appendix B. For the functions A^^ and B'^^ we have worked 
out a toy example in Appendix C with only diagonal perturbations. For off diag- 
onal perturbations we need to take an inverse of a 5 x 5 matrix to determine the 



■^^We thank Aninda Sinha for pointing this out to us. 

■^^There will be another contribution from the (p, q) seven branes in the background, although 
for small gsNf these are subleading. 
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functional form. We shall provide details of this in the following. The variables 
defined in ( p.96|) are given as: 



•^iw = ^ 2^ ^(5)fj.bcda-r(5)ub'c'd'a'g 9 9 9 9 



b,c,d,. 

n^)nabcdFi^,)n'a'b'c'd'9'''''9'''9''9'''9^''\9,u9''^-5''X) 



4 



\ E a.uF^^Fi^^g'^^g'P + ^ (3.97) 



4 

1=1 1=1 



where we have given the most generic form in ( p.97 ) above. Furthermore, and 
are operators and not functions. As usual Qab is the metric of the OKS-BH 
background. 

At this point note that the first approach of ignoring the back reaction of the 
string on the cosmological constant A(r) can only be achieved in the limit where 
there are no additional metric-derivatives and the corrections to F(5) and fJ*^ are 
very small. Otherwise this assumption would clearly fail. In this paper we will work 
out one concrete example (given in Appendix C) using the first approach only and 
give the full analysis in the sequel |]55[. Therefore in the limit where the background 
fluxes, including the effects of the D7 brane, are very small the above equation ( p.96| ) 
can be presented as an operator equation of the following form: 

n Afj^pix) ^ T^;r\x) (3.98) 

where x is a generic five- dimensional coordinate. In this final form, the operator 

is a second order differential operator derived from ( |3.95|) , with /x, z/, a, /3 = 0, 1, 2, 3, 4. 



Before we solve the above equation, we need to determine the relevant number of 
free components of Z^^. Notice that as is a symmetric tensor, it has fifteen degrees 
of freedom. We can use coordinate transformations to fix ten degrees of freedom. As 
an example consider the coordinate transformation: x"" — )■ x'"" = x"" + 9sG°', ^^lu will 
transform as: 

liMu ^ ^iMh> D^Cy DyC^ (3.99) 

where is the covariant derivative. Using this we can fix five components, say: 
= 0, /i = 0, 1, 2, 3, 4. As one might have expected for a similar example in electro- 
magnetism, this does not completely fix the gauge. The residual gauge transforma- 
tion allow us to eliminate another five degrees of freedom and we end up with five 
physical degrees of freedom i.e. five independent metric perturbation using which all 
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other components could be expressed. Alternatively one can use certain combina- 
tions of the fifteen components to write five independent degrees of freedom for the 
metric fiuctuations. 

The above result is easy to demonstrate for the AdS space, as has already been 



discussed in [^. For non-AdS spaces this is not so easy to construct. Therefore in 
the following we will take the complete set of ten components and using them we 
will determine the triangle operator /\'^ in ( |3.9(i| ). If the ten components that can 
be labelled as a set: 

In = {IqO, ^01; l02, lo3, hi, ^12, ^13; ^22, ^23; ^83} (3.100) 



then the operator AJJ^ in ( |3.96| ) gives rise to 77 equations that we present in Ap- 



pendix B. The warp factor h appearing in these equations can be taken to be 
h = h{r,n,7i) i.e (|3.12| ) with the slice condition (|3.14|) , because we will analyse fiuc- 



tuations close to the string. 

Before analysing ( |3.96 ), we make a coordinate transformation r = 1/u which will 



be useful in the wake analysis at zeroth order in Qg. Since the internal space is now 
independent of warping or more appropriately, has very mild warping (see ( |3.91| )), 
the local OKS-BH metric under this transformation will become: 



2 —g{u)dt^ + dx^ + dy"^ + dz^ \/h{u, vr, tt) ^ 

h{u, vr, tt) = u'^L^ [l — A log u — B log^w] , g{u) = 1 — 
r.-l ^^^^1 + ^^ ' (3101) 



We will now solve ( |3.96| ) order by order in QgNj and QsM /N . At zeroth order in 
QsNfjQsM'^/N, the warp factor becomes h{u) = L'^u^ and the metric p.l01| ) reduces 



to that of AdS^. Hence at zeroth order in g.sNf, QsM'^/N our analysis will be similar 



to the AdS/CFT calculations but with certain crucial differences that we will 
elaborate below. As (|3.96|) is a second order non linear partial differential equation, 
we can solve it by Fourier decomposing x, y, z, t dependence of l^^, and writing it as 
a Taylor series in u. To start off then let us decompose Z^,^ as: 

= + (3.102) 

where the subscript [0], [1] refer to the zeroth and the first order in [qsN ggM"^ / N). 
To express the zeroth order fiuctuation ffi in a Fourier series one has to be careful 
about the fact that our underlying space is not fiat. In a curved space the Fourier 
series is expressed in terms of the corresponding harmonic forms in the space. These 
forms satisfy a Klein-Gordon type equation that one could easily determine for a 
given choice of the background metrics. Although in general the spatial parts of 
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these forms could be complicated, the temporal part however is always of the form 
g-iojT Tf^i^Qj-e T is the time in the curved space. On the other hand, once we take the 
radial coordinate very large i.e r = — )■ oo, then we can take the harmonic forms 
on the boundary to be approximately wave hke with r defined as: 



r=Vgi^t (3.103) 

Of course this is not a generic picture at all points away from the boundary, but 
will nevertheless suffice for us because we will eventually address the theory on the 
boundary only. Such a procedure will help us extract the universal features of the 
cascading theory. For all other theories that we will define by specifying UV degrees 
of freedom will all inherit some aspects of the universal features (plus additional 
dependences on the UV degrees of freedom). Therefore in this limit, the zeroth 
order can be succinctly presented as a Fourier series in the following way: 



k=0 



d^qduj 



(27r)4y^ 



(3.104) 



where s^t'''°' are expansion coefficients of the solution The precise deviations 
from the wave-like behavior for ffl will not be very relevant for the present analysis. 
Similarly, we can also write the source in Fourier space as: 

jistring _ rp[0]stTing _|_ jilljstring 105) 

where as before, [0, 1] refer to the zeroth and first orders in {gsNf,gsM'^/N) respec- 
tively. The zeroth order can then be written as: 

Tj^}^'^'-^it,u,x,y,z)= I "^y^ e-(-v/^*-^^--'^"^)tM(u;,«,gi, 92,93) 

(3.106) 

where tf^l are expansion coefficients of source j'jj^J^^^^^^ ^t zeroth order in gsNf, ggM"^ /N . 
These coefficients are obtained by using explicit expressions for T^^^^^{x^) that can 
be extracted from ( ^.41| ). One may also note that since g is defined at r = rc, we will 
henceforth be analysing the theory there with appropriate degrees of freedom to be 
added later so as to have the full UV description. As mentioned before, this theory 
will correspond to a certain gauge theory with a RG fiow that would inherit some 
properties of the — )■ oo cascading theory, but in general will be different from the 
parent cascading theory (see figure 10)"^°. 



^'^Note however that even if these theories differ significantly in the UV, they all confine in the 
far IR at zero temperature. This is illustrated in figure 10 by the existence of rmin for both the 
gravity duals. 
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Figure 10: The full ten-dimensional picture of the gravity dual. The figure on the left 
is the actual dual to the cascading theory with infinite (gauge theory) degrees of freedom 
(dof) at the boundary. The figure on the right is the one with a cutoff at r = rc where we 
add a finite (but large) number of dof to define the UV of that theory. Once the UV dof 
are specified we can describe this theory at the boundary. The IR physics changes a bit 
by inclusion of irrelevant operators. Clearly the gravity dual on the right can capture only 
some universal properties of the parent theory (on the left). Nevertheless the theory on the 
right is an interesting theory (with good UV behavior) that can be studied directly from our 
analysis. The RG flows of both these theories have been discussed earlier. The existence 
of 

^min signify the confining nature of these theories in the far IR at zero temperature. At 
high temperature we should view the blackhole covering the region near r^i^. 



One may also analyse similarly the metric perturbation at linear order in gsNf 
and gsM'^/N where it is easiest to switch to coordinate u = ^ , so that the entire 
manifold (from r = to r = rc) is now described hj < ( < 1 with Tc arbitrarily 
large and we can get a meaningful Taylor series expansion of the logarithms and 
other functions appearing in the equation. As we did for the zeroth order cases, we 
can decompose the first order in QsNf, ggM'^/N fluctuations via the following Fourier 
series: 



oo „ 

/W(t,C,x,i/,z) = J2 

1 n ^ 



d^qdu) 



"J (27r)4y^ 



(3.107) 



where Sjit''^^' are the corresponding Fourier modes. These Fourier modes will eventu- 
ally appear in the final equations for 

However, as we discussed above, these modes are not all independent. There 
are only five independent metric fiuctuations which could be written as some linear 
combinations of the ten components ( |3.100| ). The zeroth and first order independent 
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fluctuations can then be expressed as: 



= Va = V / 

^ t^J (27r)V^?(^c) 

r[i] _ Vft /W - V /" ^^^^^ 

n=l k=0 



,-i{iUy^t-q,x^) {k)[0]( 



(3.108) 



where a^n and hmn are functions of all the coordinates in general with m running 
from 1 to 5. Notice that the way we are expressing the independent modes is different 



from the way it is presented in [45]. It should also be clear that 



{km 

mm 



n ?<'=)M ^ n ^'^^M Jr- n ^^^^l^l ^ n ^^^^ ^ 



(3.109) 



Since there are five independent metric components, we expect that there should also 
be five independent sources. These sources, coming from the string, should again be 
some linear combination of the ten possible components. We can write them as Tmm 
and Tmm to denote the zeroth and first orders in {gsNf,gsM'^/N) respectively. The 
above considerations also imply that the operator equation (|3.96|) will become: 

'mm (3.110) 



/ J '—'mm-^nn 
n=l 



where the form of can be easily derived from ( p.96| ). 

Knowing the metric perturbations Lmm, we can now compute the effective super- 
gravity action in five dimensions and use this action to compute the four dimensional 
boundary action. Recall that our space is a deformation of the AdS space (|3.91|) and 
from all our previous discussions, we need to be careful when we want to give a bound- 
ary description. One of the simplest boundary description that appears naturally in 
this framework is the one with an infinite degrees of freedom. The other alternative 
descriptions are defined through the cutoffs imposed at various rc- All these infinite 
possible descriptions may have finite (but very large) UV degrees of freedom and 
flow to A^eff degrees of freedom at r = Tc (see ( ^.21] )). Once we specify this we expect 
that the boundary will capture the UV of the corresponding gauge theory. In the 
limit where N » M one can even compute the spectrum of the operators from 
the boundary description (see for example ^7] for a zero-temperature example with- 
out flavor). To describe the boundary action we will in principle take a "functional 
derivative" of our action with respect to the perturbation L^m — meaning that we 
will extract only the linear coefficient of the perturbation Lmm — obtained by a fixing 
a value for the radial coordinate (u=fixed or (^=fixed). Using the boundary action 
our aim then is to compute (|3.64|) which is the wake left behind by a fast moving 
quark. 
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For comparison with AdS/CFT resuh, a simple way would be to split the five 
dimensional effective action into two parts: 



ceff _ Qeff(AdS+string) Qeff(OKS+mn) i i i \ 

'-'total — '-'total + '^total (^O.iiiJ 

where the first part is the vanilla AdS with a quark string and the second part is 
the deformation that captures the running of the gauge theory in the dual picture. 
The AdS-string part is measured at zeroth order in QsNf and QgM'^/N whereas the 
OKS-run part is measured at first order in QsNj and QsM'^/N. It should also be 
clear that in the limit {gsNf,gsM'^/N) — )■ 0, our geometry is AdS^ x T^^ Black-Hole 
with back reaction from the string^^. More concretely, for the AdS-string part we 
will use Lmm as our metric perturbation computed using the zeroth order energy 
momentum tensor of the string (|3.105| ) whereas for the OKS-run part we will use 
Lmm as the metric perturbation coming from the first order energy momentum tensor 
of the string ( p.l05| ). To avoid cluttering of formulae, we will henceforth label: 

^eff(AdS+strmg) _ ^(l) ^eff(OKS+mn) _ ^(2) /g 112) 

trOtctl ' tot Cll \ / 

Thus we can represent the wake of the quark by the following formal expressions^^: 



OhO, 



cS 



/rpmn\ _ 1- " "fe'-' total rrr 

\-L /wake — ■ ^7 J-c 



-imn 



u=0 K 6hL. 



- quark 



mm 



= hm — + - r":,k (3.113) 

which would give us the result for the UV of the dual gauge theory because C = 
would take us to the boundary of our space^^, and the operation '^'''^^ extracts the 
linear coefficient from the boundary action. The energy momentum tensor of the 
quark can then be written as: 

T^,,^{x,y,z,t) = m{T)U^U-Vl^' 5\-i - t - t) (3.114) 

where [/™ = (?7°, ^) and is the energy of the quark and is the three-velocity. 
Now at zeroth order, as mentioned before, the five dimensional effective action that 
we have is that of AdS^ plus back reaction of the string. For an infinitely massive 



^^We will discuss the back reaction from the D7 brane in the next section when we compute the 
shear viscosity. In this section we will only consider the back reaction of the string on the geometry. 

Although the energy momentum tensor in the following equation may seem diagonal in terms 
of i„i„i variables, they are in fact not diagonal in terms of the Z^j/ variables. We simply find it 
convenient to express T™" in terms of the five independent metric fiuctuations. Furthermore both 



terms of (|3.113D are restricted to four dimensional space-time. 

"^^To be more precise C = takes us to r = Tc] and then once we add appropriate UV degrees of 
freedom we can define our theory at the boundary r — ^ oo. 
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string i.e. uo = — = 0, just like in [|^, using exactly the regularisation procedure 



of [|T3] [jl4| [|T5| and [jl6[, the stress tensor evaluated at m = is known to be: 



[0] ) ~ '^^^rki^y 't) 

" ^ OhLmmJ FT 

= / ^a-^-^Hc^, ^) - T-^{^, t) (3.115) 

where FT is the Fourier transform and T^l^^y^{oj, are the quark energy-momentum 
Fourier modes and the value of Co will be discussed later. Note however that for a 
finite length string i.e. the string stretching from uq = = I to Uh the metric 

perturbation due to the string goes to zero for u < uq. This means that the previous 
analysis done in will not be vahd for our case and we would require to regularise 



our effective action. This is where it may be more convenient to regularise the system 
to 0{gsNf,gsM^/N) without splitting up the action. 

In the following therefore we will give a brief discussion of our regularisation 



scheme using the warp factor (|3.35| ). We will start by analysing the background 
string configuration by plugging in the If^i, in the Einstein as well as the flux terms of 
the five dimensional action. Once the action is expressed in terms of l^i, we expect 
that it will be equivalently rewritten in terms of Lmm- An example of this is given 
by and we provide another example in Appendix C. Once the dust settles, the 



result in the most symmetric form is given by: 

J (27r)V^W J [2 L 
+i?r(r, g)$'W(r, qWi'\r. -q) + q) [$'W(r, g)$W(r, -q) + <|.W(r, g)<|.'W(r, -q) 



+Drir, g)$W(^, -q) + rrir, g)$W(r, q) + i?rC(^, Q) + F^^f^r, q) 



(3.116) 



where m,n = 1,...,5, prime denotes differentiation with respect to r,^^; the script 
[1] denote the total background to 0{gsNf, gsM^ /N); and the explicit expressions 
for A™", 5™", C™", for a specific case are given in Appendix D. We have also 
defined $m (r, q) in the following way (with go = ^\/g(jcj as before): 

$W(r,g) = / ,f ^ e^(^°*-'?^^-«^)L^^(t,r,a;,y,;^) (3.117) 
J (27r)V^(^c) 



^"'Recall that for our case the string ended at a finite distance uq = (see (3.49) for details). 

Once the string is infinite, we are effectively putting the D7 brane at infinity i.e u — uq — 0. 

^^This is for simplicity. Prime could be more generic and denote derivatives wrt r as well as 
We will rectify this below. Notice also that we have shifted to r variable, but the action could 
equivalently be expressed by u variable. 
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We will see that the effective four dimensional boundary action is independent of D^" 
and and hence we do not list their explicit expressions in appendix. Furthermore, 



note that the derivative terms in ( |3.116| ) all come exclusively from y/—GR, whereas 
fluxes contribute powers of $^^1 but no derivative interactions (this is also one of the 
reasons why F-theory seven branes will not change our conclusions provided we ignore 
the anomalous gravitational couplings ( |3.156D or ( ^.92| )). In the following we keep 
upto quadratic orders, and therefore the contributions from the fluxes will appear 
in L)™" and E^. Taking all these affects into account we compute the coefficients 
A^"^ , B'^^ , C]^'^ , D^"- , E'^ , and for a specific example, whose values are given 
in Appendix D. However, as we mentioned before, we will be ignoring the back 
reaction of the string on the fluxes because we expect any additional components 
of fluxes coming from the string to be very small^®. The equation of motion for 
(r, —q) is given by: 



1 



Ar(r,g)$W(r,-g) 



Brir,q)^'l'\r,-q) 



Cr{r,q)^^^\r,-q) 



+ Dr{r,q)^^^\r,-q) + ^A^r, q)K^'\r, -q) - icr(r, (r, -g) 



rrir,q) - E'r{r,q) + Fr(r,g) = 



(3.118) 



The next few steps are rather standard and so we will quote the results. The vari- 
ation of the action ( p.ll6| ) can be written in terms of the variations 6^m {r, q) and 
5$n^(r, —q) in the following way"^^: 



Ct^q 



[1] 



dr 



[1] 

m 



(3.119) 



^^This is easy to motivate from the fact that the size of our quark string is very small. One end 
of the string goes into the horizon of the black hole and the other end is on the D7 brane placed 
at |^|2/3 = 1. Thus the Bns sources from the string will be negligible. Furthermore fundamental 
string will not effect any of the background RR forms, at least in the limit ( 3.50 ) that we are 
considering. So it makes sense to ignore the effects of the string on fluxes. See section 3.2 for more 
details. 

^''Henceforth, unless mentioned otherwise, $m , ' will always mean {r, q) and <i?L^' (r, —q) 
respectively. Similar definitions go for the variations (5$m and (5$L^' . 
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which includes the equations of motion as well as the boundary term. We can then 
write the variation of the action 5S^^^ in the following way: 



1 

+ 2 




^ I (EOMfor<|.W)^$W+(EOMfor<|.W)^$W 

(27r)V^(^c 



[1] 

m 



boundary 



(3.120) 



where by an abuse of notation by the "boundary" here, and the next couple of pages 
(unless mentioned otherwise), we mean that the functions are all measured at Vh and 
Tc i.e the horizon and the cut-off respectively^*. It is now easy to see why D^"^, 'J'^^ 
and E'™" etc do not appear in the boundary action. Finally, we need to add another 
boundary term to ( p.l20|) to cancel of the term proportional to 5$^. This is precisely 



the Gibbons-Hawking term EH 



(27r)4y^ 



4™"$W(T)'[1] I /1™"<T)[1]<T)'[1] I 2F-,"$'W 
1 m n '1 n m 'In 



(3.121) 



boundary 



Taking the variation of ( p.l21| ) SJCi we get terms proportional to 5^'^^^ as well as 
(5$[il Adding SJCi to 56^'^^ we can get rid of all the S^'^^^ terms from ( ^.12U| ). This 
means we can alternately state that the boundary theory should have the following 
constraints^^: 



5$'W(r„g) = S^'l'\r,,-q) = 



(3.122) 



Next, we can make a slight modification of the boundary integral ( ^.119| ) by spec- 
ifying the coordinate r to be from to certain r^ax = ''"c(l — C) instead of just 
specified earlier. For different choices of ( we can allow different UV completions at 
the boundary such that their degrees of freedom would match with the parent cas- 
cading theory only at rc(l — C)? while would be associated with the characteristic 
temperature T of the parent cascading theory. 



^^Note that we have not carefully described the degrees of freedom at the boundary as yet. For 
large enough rc we expect large degrees of freedom at the UV because of (3.23). This would mean 



that the contributions to various gauge theories from these degrees of freedom would go like e~^°" , 
which would be negligible. Thus unless we cut-off the geometry at r = Tc and add UV caps with 
specified degrees of freedom we are in principle only describing the parent cascading theory. For 
this theory of course ACg is infinite at the boundary, which amounts to saying that UV degrees of 
freedom don't contribute anything here. We will, however, give a more precise description a little 
later. 

"^^One can impose similar constraints at the horizon also. 
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With all the above considerations we can present our final result for the bound- 
ary action. Putting the equations of motion constraints on (|3.12CI|) , as well as the 
derivative constraints ( |3.122| ), we can show that the variation ( p . 1 20| ) can come from 
the following boundary 3 + 1 dimensional action: 



5(1) 



(27r)4y^ 



Cr(r, q) - A'rir, q)] $W(r, g)<|.W(r, -q) 



Brir, q) - Arir, q)\ [$'W(r, g)$W(r, -q) + $W(r, q)^':^\r, -q) 



+ [Er-Fry^{r,q) 



(3.123) 



However the above action diverges, as one can easily check from the explicit expres- 
sions for A™", 5^", C{"", E"™ and for the specific case worked in Appendix D. 
Indeed, comparing it to the known AdS results, we observe that from the boundary 
there are terms proportional to in each of C^"' , A'(^"' , and F["^ and propor- 
tional to in A^^, B^^. As Tc — > oo the action diverges so as it stands Vc cannot 
completely specify the UV degrees of freedom at the QCD scale Ac^°. Thus we need 
to regularise/renormalise it before taking functional derivative of it. This renormal- 
isation procedure will give us a finite boundary theory from which one could get 
meaningful results of the dual gauge theory. 

Once we express the warp factors in terms of power series in r(^a) ( p.35[) the renor- 
malisability procedure becomes much simpler. Note however that this renormalisa- 
tion is only in classical sense, as the procedure will involve removing the infinities in 
( p.l23| ) by adding counter-terms to it. Comparing with the known AdS results, and 
the specific example presented in Appendix C, one can argue that the infinities in 
( P.123|) arise from the following three sources: 

1. Cr{rc,q)-A'r{r,,q) = Y.HvA{q) rt^^.) + finite terms 

a 

2. Sr(rc,g)-Ar(r.,g) = 5^i^|':r(?)^cV) + fimte terms 

a 

3. Er(rc,g)-Fr(r„g) = Y.W'^ti.) + finite terms (3.124) 



where in the above expressions we are keeping a arbitrary so that it can in general 
take both positive and negative values; and the finite terms above are of the form 
^c(a) "with n > 1. Therefore to regularise, first we write the metric perturbation also 



50 



We are a little sloppy here. Of course our theory is large N QCD like only in the far IR. 



Nevertheless should specify the UV degrees of freedom as (3.23) at scale Aj, for the resulting 
theory. 
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as a series in 



oo (k)[oi] 



(3.125) 



fc=0 a 



where the above relation could be easily derived using ( 3.117|) , ( 3.108 ) taking the 
background warp factor correctly; and tq = l/Up/^ is given by (|3.49|) . 

Plugging in ( |3.125| ) and ( p.l24| ) in ( p.l23[ ) we can easily extract the divergent 
parts of the action ( |3.123|) . Thus the counter-terms are given by: 



S. 



(1) 



counter 



2(27r)4^ 



Mm Mb] ^4 



mm nn {a) 



r7^^ + s 



'mm ^nn ^(ai) 



ct,p,-f 



(0)[/3]«(l)[7]^3 



mm '^nn {0.2) 



J2)[/31 (0)[7U2 



mm nn (0^3) mm nn (cf4) mm nn (c«5) 



(0)[/3]^(2)[7]^2 



(1)[/3]«(1)[7V2 



J3)[/3] (0)[7] 
"'mm "-^nn 



' ("ej ^ ^mm •^nn ' (a?) ^ ^mm ^nn ' (ag) ^ '^mm ^nn } ' (agj 



.(2)[/3],(l)[7L 



,(1)[/3]J2)[7] 



+ir|7^(ro - r) 



_2„(0)[/3] (2)[7] 2 , 2s(°)[^ls(2)Wr? ^ 

mm nn (ais) nn mm (014) 



('J0)[/3] (1)[7]^3 , „(l)[/3] (0)[7]^3 \ _ f 9 .(1) Wr^ . 
V mm "'nn (ckio) mm "'nn I "-^mm ^nn (cn2) 

'2s(l)[/3]„(2)[7U, , I .(2)[/3]„(l)[7U, , 

I Qo(0)[/3]„(3)[7U, , I 2s(l)[^ls(2)[7]^, , I ,(2)[/3]„(l)[7U, , 1 Q.(0)[/3] „(3)[7]^^ \" 
~'~"^m.m ^nn ' {ai7) ^ "^^nn ^mm ' {ais) ^ ^nn ^mm ' {aig) ^ "^^nn ^mm ' (a2o) J 

_L /-."^/^r^^ _ ^\/^o(0)[/3V4 I Jl)[/3] 3 , „(2)[/3] 2 , J3)[/3] \\ |"? 1 9fi^ 

+ -'|arVO '^j^^'Smm ^(a) + "^^m ^(ai) + "^mm ^(0:3) + -^mm l^O.iZDJ 

with an equal set of terms with In the above expression r(aj) = r^"^'"^'; and 

as before, the integrand is defined at the horizon rh and the cutoff rc(l — C); with the 
string stretching between rh and tq. The other variables namely, s^m\ H'l^^ , 
and /|™| are independent of r but functions of g*. For one specific case their values are 
given in Appendix C. Finally the e(Q.) can be defined by the following procedure. 
Lets start with the expression: 

4 

rrmn (a)[/3] (6)[7] p ^ rrmn [aW] {b)h\ ^(«) 
-'-'|q| ^mm ^nn ' (aj.) — ^mm ^nn a „b 



T^mn (cm {dM q ^ I^mn (c) (d) [7] (") 
^^\ol\ ''mm ^nn ' (ai) — ^^\a\ ''mm "nn c 



(3.127) 



from where one can easily infer: 

p = 4-a-b, e(«,) 
q = 5- c - d, e(a,) 



4e(a) — ae(/3) — 6e(^) 

A - a-b 
5e(a) — ce(^) — (ie(^) 

5 — c — (i 



(3.128) 



Using this procedure we can determine all the r^^fi) in the counterterm expression 

(ra- 
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At this point the analysis of the theory falls into two possible classes. 

• The first class is to analyse the theory right at the usual boundary where — )• 
oo. This is the standard picture where there are infinite degrees of freedom at the 
boundary, and the theory has a smooth RG flow from UV to IR till it confines (at 
least from the weakly coupled gravity dual). 

• The second class is to analyse the theory by specifying the degrees of freedom at 
generic and then defining the theories at the boundary. All these theories would 
meet the cascading theory at certain scales under RG flows. The gravity duals of 
these theories are the usual deformed conifold geometries cutoff at various Vc with 
appropriate U V caps added (of course for r < Vc the geometries change accordingly) . 

The former is more relevant for the pure AdS/CFT case whereas the latter is 
more relevant for the present case^^. 

For the pure AdS/CFT case without flavors e^ai) = (so that the subscript a^'s 
can be ignored from all variables), we can subtract the counter-terms ( p.l26| ) from 
the action (|3.123|) to get the following renormalised action: 



c{i) ^ c(i) _ c( 

^rcn ^ '-^counter 



+s. 



e(rn - r) [if"^"!'^^^) + I o(2) „(2) , (1) (3) 

'^V U ' / y'-'mm'-'nn ' '-'mm'-'nn ' '•'mm'-'nn ' '•'mm'-'nn 

mm nn I \ mm nn ' mm nn ' mm nn ' mm nn 

A.W^i^) I Q.(l)o(3) I „(3) (1) \ , jm (4) 1 /o loq^ 

^nn^mm ' '-''^nn'^mm ' ^nn mm I ' "^mm ytj.^^^uj 



where we have made all the 0{l/rc) terms vanishing, and in the limit Vh small the 
small shifts to Snn given by Snl + 0{rf^) can also be ignored. Furthermore for more 
generic case where primes in ( p.ll6| ) denote derivatives wrt r as well as we can 
reinterpret the renormalised action ( |3.129| ) as the following new action: 



^(1) / 



Z-"<|.^(g)<l>„(-g) + t/™"(<l>^(g)$'„(-g) + <l>'^(g)$„(-g)) 



(2vr)4 

+ F"<l>™(g) + F"$„(-g) + V"'<^'^{q) + l^"<(-g) + X 



(3.130) 



where $m are now only functions of ±g with prime denoting derivatives wrt and 
we isolated all the r dependences so that X, Y, Z, U, V could be functions of r and 
'f. We can determine their functional form by comparing ( p.l30| ) with ( |3.129| ). For 



us however the most relevant part is the energy momentum tensors which we could 
determine from (|3.130|) by finding the coefficients and Y^. One can easily show 



^^In both cases of course we need to add appropriate number of seven branes to get the finite 
F-theory picture. The holographic renormahsation procedure remains unchanged and the far IR 
physics remains unaltered. The UV caps affect mostly geometries close to r^, as expected. 



- 63 - 



that, upto a possible additive constant, F™, are given by: 

from where the wake of the quark can be shown to be exactly given by ( |3.115| ) 
mentioned earlier. We also note that Cq = H'^^6{rQ — r) there. 

Now let us come to second class of theories wherein we take any arbitrary r = 
Tc, with appropriate UV degrees of freedom such that they have good boundary 
descriptions satisfying all the necessary constraints. For these cases, once we subtract 
the counter-terms ( p.l26| ), the renormalised action (specified by r) takes the following 
form: 

c(i) = c(i) _ c(i) 

^ren ^ "-^counter 



2(27r)4 



I „(2)[/3]„(2)[7V2e(^)+2e(^)-4e(„) , „(l)[/3] „(3)[7Lf(^()+3e(^)-4e(„) , „(0)[/3] (4)[7] 4e(^)-4e(„) 

_AT^mn( (0)[p] Ji)['y] he^^^)-bei^a) i „(3)[7] r^e(/3)+4e(^) -5e(„) , 2e(^)+3e(^)-5e((:,)l 

_|_4o(2)[/3] e(2)[7] r^2e(^)+3e(^)-5e(a) , 3e(0)+2e(^) -5e(„) 1 , 4„(0)[/3] (4)[7] 5e(^)-5e(„) 
I o(3)[/3] „(1)[7] r^4e(^)+e(^)-5e(a) i 3e(^)+2e(^) -5e(Q) 1 ^ , rm (4) [/3] 5e(^) -5e(Q) 1 j-g i go\ 



defined at the cut-off and the horizon radii as usual. Notice now the appearance of 
j,me(„)+ne(^)+pe(^) fg^Q^Q^g^ Quc cau casily show that: 

+^-™^(.)-n.(«-peMj _ I ^ 0[e(c.,iir,)\^ (3.133) 

Since the warp factor h is defined only for small values of QgNf, QgM'^/N, g'^NjM^ /N 
we don't know the background (and hence the warp factor) for finite values of these 
quantities. Therefore for our case we can put O [€(^a,i3,'y)\ to zero so that the value in 
( p.l33| ) is identically 1. For finite values of these quantities both the warp factor and 
the background would change drastically and so new analysis need to be performed 
to holographically renormalise the theory. Our conjecture would be that once we 
know the background for finite values of QsNf, QsM'^/N, the terms like ( |3.132| ) would 
come out automatically renormalised by choice of our counterterms. We will discuss 
this in some details later. More elaborate exposition will be given in the sequel 



Once this is settled, the renormalised action at the cut-off radius would only 
go as powers of r'^^y Thus we can express the total action as: 




oo \ / oo ~{") 

mn{j) \ p^mn^ ^ , ( \^ rnn{j) 
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T^mn 



I ^rnrn ^ ri 



+ 



E 



m(j) 



0=0 '(q) 



4s(4)[/3],T, _4s(4)[/3]$ I s(5)[/3U' 



where are independent of r with prime denoting derivatives wrt ~^ henceforth; 
and the radial coordinate is measured at the horizon r/j and the cutoff rc(l — () 
as before. The exphcit expressions for the other coefficients hsted above, namely, 
Qmn j\^"in ~(") jm. b^"} cau be worked out easily from our earlier anal- 

ysis (see Appendix C for one specific example). Note that X[r(„)] is a function 
independent of $m and appears for generic renormalised action. 

Now the generic form for the energy momentum tensor is evident from looking 
at the linear terms in the above action ( p.l34| ). This is then given by: 



rpmm 

-'o 



(2^ 



( TTmn _|_ TTnm,\ (i)[l3] _ A/jy-mn , T^nm\ (i)[l3] , ( T^mn , ]^nm\ 



o da) 



1«l 



H 



■H 



■\a\ 



,(5)[/3] 



PS 



(3.135) 



at r = Tc (we ignore the result at the horizon) and sum over [a, (5) is implied. This 
result should be compared to the ones derived in [jl3| |]14| [jl5| [|16[ which doesn't 
have any dependence. This is of course the first line of the above result. For 
the case studied in 



^ [|T^ ||T5| [|T6] the boundary theory is defined with infinite 

degrees of freedom at UV (which we have been calling as the parent cascading theory) . 
How do we then reproduce the results of those papers? Before we go about elucidating 
this, notice that the second line also has a independent additive constant. This 
additive term is irrelevant for our purpose because the energy-momentum can always 
be shifted by a constant to absorb this factor. Thus once we specify the cutoff Tc and 
the UV degrees of freedom, then our result shows that the energy-momentum tensor 
not only inherits the universal behavior of the parent cascading theory but there are 
additional corrections coming precisely from the added UV degrees of freedom at 
r = Tc- These corrections go as r^^ or as e~-^'=** with being the UV degrees of 
freedom at the cut-off (see ( p.23| )). As long as these UV degrees of freedom are not 
infinite, they contribute to small corrections to the energy- momentum tensors of the 
various gauge theories. 

The above description is one of the key points of our paper, and tells us how we 
can distinguish our results from the standard AdS/QCD answers. Therefore let us 
elaborate this in little more details. This will justify what we have been saying so far 
about UV caps, and put everything in a rigorous mathematical framework at least 
from the action point-of-view. In the process we will also be able to reconcile with 
the results of [|T^ ||T3| |TB|. The first important issue here is that we can study 
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infinite number of UV completed tlieories in our full F-theory set-up. All of these 
theories have good boundary descriptions and have same degrees of freedom as the 
parent cascading theory at certain specified scales. The simplest UV complete theory 
(however with infinite degrees of freedom at UV) is of course the parent cascading 
theory that we will discuss in a moment. The question now is to construct other 
possible theories by defining the degrees of freedom at the boundary. To do this 
observe that we defined the boundary theory using the identification: 

where the boundary is at r — )■ oo. For the boundary cascading theory the above 
expression simply means that 



(271) 



J 



+ — factors (3.137) 

00 

where the sign is crucial and sum over a is again implied (note (a) the cut-off de- 
pendence, and (b) rc(^) is some function of Tc^a) that one can determine easily). We 
now see that the contributions from the UV cap give the following values for Bf 

-°n(0) ~ -°n(l) — -°n(2) " -°n(3) ~ ^ 



1 

n(4) ~ 2 



b%ymr, - r.) + 0{H1^^, Kl^^, .M) 



(3.138) 



this means that in the limit of small rh we can ignore the contributions coming from 
. For this section such an assumption will not change any of the results, so we 
will stick with this. From the next section onwards, we will restore back the 0{rh) 
dependences. 

Taking all of the above considerations, this implies that the contribution from 
r > Tc exactly cancels the (9(l/rc) contributions coming in from the action measured 
from r/j < r < Tc, giving us the boundary energy-momentum tensor 



(2^ 



(3.139) 



which is the result derived in |45] [lilllSl |16|. The above way of reinterpreting 



the boundary contribution should tell us precisely how we could modify the boundary 
degrees of freedom to construct distinct UV completed theories. There are two 
possible ways we can achieve this: 
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• From the geometrical perspective we can cutoff the deformed conifold background 
at r = Tc and attach an appropriate UV "cap" from r = rc to r oo by carefully 
modifying the geometry at the neighborhood of the junction point. As an example, 
this UV cap could as well be another AdS background from to r — oo. There are 
of course numerous other choices available from the F-theory limit. Each of these 
caps would give rise to distinct UV completed gauge theories. 

• From the action perspective we could specify precisely the value of the action 
measured from Tc to r — )■ oo, i.e r5renl°°. The simplest case where this is zero gives 

L J Tc 

rise to a boundary theory whose degrees of freedom at the UV is the one given in 
( |3.21|) . We will give an example of this towards the end of this section. To study 



more generic cases, we need to see how much constraints we can put on our integral. 
One immediate constraint is the holographic renormalisability of our theory. This 
tells us that the value of the integral can only go as powers of l/rc{a) otherwise we 
will not have finite actions. This in turn implies^^: 



(a) 



x(<l>^<l>'„ + <l>:,<l>„) + -j^ J"'^^ + ^h")] r(^o - re) + finite terms 

where by specifying the coefficients -^mn(j) ^^"^ '"^^ specify the pre- 

cise UV degrees of freedom! The finite terms are r^. independent and therefore would 
only provide finite shifts to our observables. They could therefore be scaled to zero. 
Notice also that the contributions from (|3.140|) only renormalises the coefficients 



a 



mn(j)' ^n{j) ^^"^ ^m{j) (|3.134|) , and therefore the final expressions for all the phys- 
ical variables for various UV completed theories could be written directly from ( p. 134] ) 
simply by replacing the l/r^ dependent coefficients by their renormalised values. This 
is thus our precise description of how to specify the UV degrees of freedom for various 
gauge theories in our setup (see figure 11 below). 

Once the UV descriptions are properly laid out, we can determine the form for 
y4(")^^(a) Qjid E^""^ by writing Callan-Symanzik type equations for them. These are 
classical equations and do not capture any quantum behavior. Nevertheless they tell 
us how A^°'\ 3^°"^ and E^"'^ would behave with the scale Tc or equivalently Hc- For A 
the equation is^^: 

^^(°) 

/^c^f^ = e^] [iLI,, + aLloJ (3.141) 



^^Rcmcmber that there are additional 0{rh) contributions to the coefficients. For small rn they 
will only make the UV contributions more involved without changing any of the underlying physics 
as we saw above. We will therefore ignore them. 

^^The following equation is derived from the scale-invariance of [iSron] ^ . 
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s- X 



X 



finite UV degrees of 




Figure 11: The equivalence between two different ways of viewing the boundary theory 
at zero temperature. To the left we add finite UV degrees of freedom at r = rc of the 
deformed conifold geometry. Such a process is equivalent to the figure on the right where 
we cut-off the deformed conifold geometry at r = rc and add a UV cap from r = Tc to 
r = oo. The boundary theory on the right has A^uv degrees of freedom at r = oo and all 
physical quantities computed in either of these two pictures would only depend on A^uv 
but not on r = Vc- At non-zero temperature the UV descriptions remain unchanged. Of 
course, for both cases one need to put in the corrections to IR from UV deformations. 



with similar equations for i?^"-* and E^"\ These equations tell us that physical 
quantities are independent of scales. The parent cascading theory is defined as the 
scale- invariant limits of ( |3.14ip , i.e: 



A^"\., = ^^"V^ = B^'^l, = -6^"?., (3.142) 

The above relation gives us a hint how to express and in terms of TVeg, 

the effective degrees of freedom at r = Tc and Mm, the effective degrees of freedom 
at r = oo i.e the boundary: 

mn(j) mn(j) mn(j) 

mn(j) mn(j) mn(j) 

Btl = -b':l, + ^(V'f^"^"^"^^"'^^"^ (3-143) 

where the actual boundary degrees of freedom are specified by knowing amn{j), ^mn{j) 
and bm(j) as well as M^. Since j goes from to oo, there are infinite possible UV 
complete boundary theories possible^^. For very large Muv (i.e Afuv — ^ e^",n >> 1) 



^^The connection of 7 with UV completions come from the coefficients a^"\ e^'^\ and s 
etc. that depend on j. For different choices of these coefficients we can have different UV comple- 
tions. In this sense j and UV completions are related. See also the F-theory discussion presented 
towards the end of section 3.1. 
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the boundary theories are similar to the original cascading theory. The various 

choices of {d^ml{j)C^)y ^ml{j)C^)'^^m{j)C^)) ^^11 ^o"^ degrees of freedom change 
from A/uv to A/'efT under RG flow. The ~^ dependence of all the quantities will tell us 
how the UV degrees of freedom affect IR physics. This is to be expected: addition 
of irrelevant operators do change IR physics, but not the far IR^^. 

Therefore with this understanding of the boundary theories we can express the 
energy-momentum tensor at the boundary with Afuv degrees of freedom at the bound- 
ary purely in terms of gauge theory variables, as: 



rpmm 

-'o — 



oo 



+ E %(^)'^"e"^-^"^H (3-144) 
i=o 

where sum over a is again implied, and the first line is the universal property of the 
parent cascading theory inherited by our gauge theory. This part, common to all the 
theories, will be numerically different. The second line specifies the precise degrees 
of freedom that we add at r = rc to describe the UV behavior of our theory at the 
boundary r — )■ oo. Using this procedure, the final results of any physical quantities 
should be expressed only in terms of Mm i-e the UV degrees of freedom^^. 

The above analysis of holographic renormalisation tells us something very inter- 
esting. The final result of the renormalisation procedure is almost identical to the 
renormalisation procedure done by taking the highest positive integer power of r and 
ignoring the log r part! The only difference is the presence of coefficients like Hl^^l^ 
that depend explicitly on flavor degrees of freedom. This will be useful for us when 
we study shear viscosity in the next section. 

Another question to ask regarding holographic renormalisation is the issue of 
corrections for finite values of QsNf^ggM'^/N. This is probably more difficult to 
tackle. However an easier question would be to ask what happens if the UV behavior 
is governed by instead of for a certain choice of an UV cap? 

To evaluate this observe that we can take the following ansatze for the back- 



^^The ~^ dependences of the UV caps are also one-to-one correspondence to the changes in the 
local geometries near the cut-off radius as we discussed before. All in all this fits nicely with 
what one would have expected from UV degrees of freedom. Of course it still remains to verify the 
story from an actual supergravity calculation. We need to analyse the metric near the junction by 
studying the continuity and differentiability of the metric and see how far below r — rc we expect 
deformations from the UV caps. Various types of deformations will signal various set of irrelevant 
operators. Needless to say, the far IR physics remain completely unaltered. 

^^Restoring back the 0{rh) contributions would mean that there should be an additional contri- 
bution to ( p4^ of the form J2T=o '^(^""^ ' -^"P ' ' -S"" )T'^^' e-^^"" where G is a function whose 



functional form could be inferred from the UV integral ( f3.140 ). 
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ground warp factor for large r on the slice ( |3.14 )^'': 



r4 

Mr) = -r 



m,n,p 



N 



^4 



Po 

p=i 



Cp 



(3.145) 



where po is some given integer, Cp are functions of QsNj, QgM'^ /N and the arrow is 
motivated by the fact that, in the presence of a suitable UV cap the background is 



given by F-theory construction that guarantees that the warp factor goes from 
log r at small r to inverse powers of r at large r to avoid the Landau poles. 
In this background we expect the fluctuation to go like: 



(3.146) 



which is similar to the definition of $m {r, q) given earlier in (|3.117|) with go = 
u^^/gifc) as before. Our next step would be to write the above action as an EOM 
part and a boundary part like as in (|3.12CI|) with appropriate changes. Once we do 



this we would require the Gibbons-Hawking boundary terms to cancel any unwanted 
parts from the action. 
Taking all these into considerations the boundary action takes the usual form 
( |3.123[ ). Now the coefficients C^^ , A^"- , B^"- , and are again divergent, and 
from our previous consideration we know that the divergences will be controlled by 
the highest integer power of r. Therefore we need to regularise and renormalise the 
action again. The divergences are: 



po 



1. Cr{Cq)-A'r{Cq) 



y i/7" (g) rP(l - C)^ + finite terms 



p=i 

po 



2. 5r(C, q) - AnC, q) = Yl ^r(l) - O'^' + fimte terms 



p=i 

po 



3. ETiC, q) - EnC, q) = ^p"(^) ^"(^ " 0' + finite terms (3.147) 



To start off then, we write the metric perturbation $ 
( is defined with the equation r = rc(l — Q: 



[1] 



as a power series in ( where 



$[;.'(c) = ^(c-co)E^ 



(fe) 

Smm 



k=0 



cy 



(3.148) 



^^The equation (3.145) could be motivated by adding a suitable UV cap that allows an inverse r 
dependences in (3.35) which in turn stems from the finiteness of F-theory. At large r the axio-dilaton 



contributions will come from (3.40) that will back react on the geometry to give a finite asymptotic 
geometry. There are infinite such geometries possible. For a given geometry the coefficients of r~P 
terms would be fixed accordingly (by solving the sugra equations with axio-dilaton and three and 
five-fluxes as sources) provided these fluxes are also asymptotically finite. For more details see the 
F-theory arguments presented earlier. 



- 70 - 



which can again be easily derived by plugging (|3.108|) in (|3.146|) . 

Once we have the mode expansions of the field, we can use this expansion directly 
in the action ( |3.116| ) to determine the possible counter-terms. The infinities come 
from ( 3.147|) , much like the situation with ( 3.124 ). Therefore its no surprise that we 
get the following counter-terms: 



S. 



(2) 

counter 



2(27r) 



E {^r^(C - Co) - C)^ 



' » mm nn ' 



P 

(0) ^(i)Up-Vi 



mm"nn I c 







P-1 + ( ,(2) ,(0) , „(0) (2) 



I „(l)A^p-2/i _ A\p-2 , / (3) (0) , (0) (3) , (2) (1) , (1) (2) 

' mm nn I c V S7 ' \^mm nn ' mm nn ' mm nn ' mm nn 

x(i - c)^-^ + ....] + i^r^(c - Co) [ - {s^±s^^i + s^^s'^;)) rr\i - cy-' 

- ( 2s(^) s^^) + 2s(°) s(2) -I- 2s(°)s(2) V^-^fl - CY-^ - ( 2s(^) s^^^ + s^^^ s^^^ 

\ '^mm'^nn ' '^mm'^nn ' -^nn'-^mm J c \ S/ \^'^mm'^nn ' '-^mm'^nn 



,p-3 



+ 3,(0) .(3) I 2s(^)s(2) +3s(°)s(3) IrP-^n - + 

' '-''^mm^nn ' '^nn'^mm ' '^nn'^mm ' '-''^nn'^mm j c \ SJ ' 

+ /r^(C - Co) [si^Lr^(l - C)^ + ^^^Lrri(l - C)^^^ 

~^ ^rnm''"c ~ C)^ "I" ^rnm''^c ~ C)^ ' "I" •••• ^ 



(3.149) 



with of course another set of counter-terms defined at the horizon Vh that we do not 
write here. Observe that we have written the counterterms wrt the highest power of 
p. Finally the renormalised action after some manipulations can be written, similar 
to ( p.l34| ) but with different coefficients (whose explicit values will not be relevant 
for the present 



5(2) 



Cmn{j) 



fmn{j) 



0=0 



1 

ri 



^m^nj -'-'p l^nn^m ^ ^mm^n 



4e(P)<T) _ 4e(P) $ + s(P+l)<|)' 

nn ''^ mm 'n ' nn m 



' "mm n 



+ 



■m{j) 



Q™<l>„ + r(r)^0(C-Co) 



(3.150) 



from here we can extract the energy-momentum tensor T"^'" once we add the right 
boundary degrees of freedom. This is exactly as we did before. The contribution 
from UV cap now will be: 



oo 
renj ^ 



dS 



oo ^ 



(V 



0=0 



+ ^'m^n) 



D 



m(j) 



Q*"*™ ^^(C - Co) + finite term£(3.151) 
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where by the range C to oo we mean the coordinate range rc(l — C) to r — )■ oo. As 
before, specifying the quantities C^°'\ F^"^ and D^°'^ will in turn specify the boundary 
degrees of freedom. We also expect similar Callan-Symanzik type equations for them. 
Indeed: 



f^C = J 



(3.152) 



with similar relations for the others. This way we can compute the the total energy- 
momentum tensor of the system once we know the precise value of po, Cp for the 



exact background. On the other hand, for the case that we know very well (|3.11|) , 
the energy-momentum tensor that we should substitute in ( p.64| ) to compute the 
wake will be: 



Mm 



rpmm f ^ 1 \ ^ J ( Tjmn rrnm \ (4) A(K'"^'"'-l- T<r'^"^\Q^ 

medium+quark 

oo "\ 

+ (i^i^p + + bl^J^'Snme-'^-^ + OiTe-^-^) \ (3.153) 

i=o J 

where (^i"]); t^i°]))-Miv) together will specify the full boundary theory for a specific 
UV complete theory. Observe that the result is completely independent of the cut-off 
that we imposed to do our analysis. In the limit A/'uv e~" with n >> 1 i.e when 
the boundary degrees of freedom go to infinity as this way, we reproduce precisely 
the result of the parent cascading theory. 

Before we end this section let us make the following observation. This is a slightly 
more non-trivial example where A/'uv ~ -^cS with A/'eg being the degrees of freedom at 
r = Tc given earlier in ( |3.21| ). This means that the degrees of freedom don't change 



significantly as we go from far UV to the scale Ac. From the gravity dual perspective 
this is like adding a UV cap given by an AdS geometry. For such a theory ( |3.153| ) 
can be exactly determined (both d^n^j) can be scaled to identity), and so precise 
prediction can be made^^. 

In the next two sections where we deal with shear viscosity and viscosity-to- 
entropy ratio, we will find that all the techniques that we discussed in this section 
will become very useful. 

3.4 Shear Viscosity 

In this section we will compute the shear viscosity of the four dimensional theory 
following some of the recent works ^ ^ . Our basic idea would be to use the Kubo 

^® Interestingly for small enough Ac we are almost in the IR where we expect the gauge theory to 
confine (at least at zero temperature) . The UV of the gauge theory has very slow running, which 



can be constructed from our F- theory model by suitably choosing the coefficients ai,&j in ( 3.40 ) 
such that axio-dilaton vanish. With vanishing axio-dilaton the geometry at UV approaches AdS. 
Therefore this set-up is almost like large N flavored QCD with vanishing Beta function at UV and 
confinement at IR! Our method could then be used to evaluate the thermal quantities of this theory. 
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formula |60|1 : 

V=\im ^ [ dtd'x e*-*([T23(x),T23(0)]) = - lim ^''(^' (3.154) 
t^->o Zuj J t<j-5>o UJ 

where G^{u, 1^) is the momentum space retarded propagator for the operator T23 at 
finite temperature, defined by 

^'{uj.t) = -^J dtd^x e'("*-^-^)^^(t)([T23(x),T23(0)]) (3.155) 

In the following, we will compute the Minkowski propagator (T23(t, "a^)T23(0, a^) 
using gauge/gravity duality^^. However before we compute this explicitly, let us take 
a small detour to evaluate the higher order corrections to the effective action from 
the wrapped D7 brane in our theory. We have already given a brief discussion of this 
in (3.921). It is now time to deal with this in some more details. 



In the case of a D7 brane, the disc level action contains the term 61 



where C4 is the four-form, R is the curvature two-form, (f) is the dilaton and is 
a non-trivial eight manifold forming the world-volume of the D7 brane. The above 
action is SL{2,Z) invariant which one can show by doing an explicit analysis [|61|] . 
Since for our case the D7 wrap a non-trivial four-cycle, we can dimensionally reduce 
it over the four-cycle to get the following action: 

S^'^^ = ^—[ Re\logr](T) ti (R A *R-iR A R)] (3.157) 
lovr^ J Mi 

where //(r) is the Dedekind function, and r is the modular parameter defined by: 

^ = // AO ( / ^4 + ^V4 ) =-{n+ IT2) (3.158) 

with V4 being the volume of the four-cycle on which we have the wrapped D7 brane. 
An interesting point here is that the above action can be derived from the following 
action that has two parts, one CP-even and the other CP-odd 



I loMr)\hT {RA*R)--^ [ log^tr (i?Ai?) (3.159) 
Jm4 327r^ Jm4 r][T) 



^^Note that our prescription ( 3.65 ) computes a path integral {0(j))j- ^ J [VO]j' e^M o<t>+c ^ (where 
T is the time ordering) with a classical action S'sugra, unaware of the time ordering. Therefore 
computing any commutator is rather subtle here. For this reason, we compute only the correlator 
(r23(T, ~af )T23(0, ~af )) and not it's commutator. Using this correlator, we will eventually obtain 
information about commutators and finally the viscosity 77. 
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where the first part is CP-even and the second part is CP-odd. To compare ( |3.159| ) 



with (|3.157|) note that the Dedekind t] function has the following expansion in terms 
of g = e^""'^: 



log|r/(r) 



TT 



3g2 4g3 

— + — + 
2 3 



rj{f) 6 



3g2 4g3 



c.c 



.. — c.c 



(3.160) 



Combining everything together we see that, upto powers of q ( |3.159D and ( ^.157] ) 
are equivalents^. However writing the action in terms of (|3.159|) instead of ( p. 157 ) 
has a distinct advantage: from D7 point of view ( p. 1591) captures the D3 instanton 
corrections in the system |^2|, But there is a even deeper reason for writing the 
action as (|3.159 ). The CP-even and CP-odd terms can be — in the case where the 
space is not a direct product of D7 world- volume times and normal space — expanded 
further to incorporate Gauss-Bonet type interactions in the following way |]62|: 



CP— even 



24 



(3.161) 



where ai is a constant, and £gb and Cr are respectively the Gauss-Bonnet and the 
curvature terms defined in the following way^^: 



= g^^2 {^RaM^R"^'^'^^ ~ 2-Ra/3-R"'^ " Rah-fSR"'^'^^ + 2RabR"^''^ (3.162) 

In the above note that the three curvature terms RapR"'^, RabR'^^ and i?^ are not the 
pull-backs of the bulk Ricci tensor. Furthermore we have used the notations (a, /3) 
to denote the world- volume coordinates, and (a, h) to denote the normal bundle. 
Thus from the CP-even terms, the coefficient of R^p^sR^^^^ is given by: 



C3 



327r2 



(3.163) 



which has an overall plus sign because ai in many interesting cases tend to be zero 
(see ^ ^ 

If we now compare this to 



62| for details on this). However in general for certain exotic compactifications 

47r4Q'2 



we expect ai « 



we see that C3, which is the 



^°Note that q can be made small because both C4 and V4 in ( 3.158 ) can be made small. 

^^The Gauss-Bonnet term is in general a topological invariant in four dimensions, but it is a total 
derivative at quadratic order in all dimensions. Therefore it wouldn't contribute to the equation of 
motion. 
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coefficient of Rap-ysR"'^^^ in [0, is positive. This would clearly mean that adding 
fundamental flavors lowers the viscosity to entropy bound!^^ 

The CP-odd term on the other hand has a standard expansion of the following 



form for the D7 brane [|^, 61, 62 



^CP-odd = Tr {Cs + MiRAB^ (3.164) 

where the ffist term gives the standard dual axionic charge of the D7 brane. Combin- 
ing ( p.l61| ) and ( |3.164| ) we get the full back reactions of the D7 brane upto 0{a'^). 

Having computed the back reactions of the embedded D7 brane, we can use this 
result to compute the shear viscosity. There are three scenarios from which we can 
compute the viscosity using the gravity dual now: 

• Allow a gravity dual that has no running (i.e no RG running in the gauge theory 
side), but incorporates the back reaction of the D7 brane ( |3.161| ). 



• Allow a gravity dual that shows the RG running of the gauge theory, but does not 
incorporate the back reaction of the D7 brane (although may allow the CP-odd part 

(EH)- 

• Allow a gravity dual that not only shows the RG running of the gauge theory, but 
also incorporates both the CP-even as well as CP-odd parts ( p.l61|) and (|3.164|) . 



The ffist part has recently been done in [0, H, |[. The background remains AdS, but 
now there would be terms from the D7 brane ( |3.161| ) that would lower the viscosity 
to entropy bound. However the second and the third part is not yet been addressed 
in the literature. As we show below, the second part is rather easy to do following 
the calculations of because we have already constructed the background in the 
previous section. 

The third part on the other hand is rather subtle because incorporating the 
higher order (curvature)^ corrections would change all the results that we derived 
earlier. In particular the mass and drag of the quark would need to be modified 
alongwith the wake of the quark. Plus the holographic renormalisability would also 
get modified because of the extra derivative terms from these corrections. In this 
paper we will deal mostly with the second part of the above list and some calculations 
of the third part, but leave a more detailed analysis of the third part for the sequel 

To start off the analysis, we need the correlation function of T2z{x) and T23(0) 



to use it in the Kubo formula ( |3.154 ). From gauge/gravity duality we know that 



switching on T23 in the gauge theory is equivalent to switching on graviton mode 
along x"^ = X and x^ = y directions. Thus in the given OKS-BH background ( |3.91| ), 



^^We were motivated to carry out the above analysis from a comment by Aninda Sinha. His 
paper Q dealing with the violation of viscosity to entropy bound has appeared recently and has 
some overlap with this section. 
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we switch on the following off-diagonal part: 



fi'oo 9ox goy goz gor 

gxo gxx gxy gxz gxr 

gyO gyx gyy gyz gyr 

gzo gzx gzy gzz gzr 

yfl'r'O grx gry grz grr j 



(-gi.f) 



1 

7^ 
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(r, t) 
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1 








(3.165) 



9{r)/ 



where h = h{r, vr, tt) and the internal space metric is independent of the radial coor- 
dinate r as shown in ( |3.91| )^^. Note that the off-diagonal gravitons are propagating 
in OKS-BH geometry with energy u. 

Since our goal is to compute the Fourier transform of (T23(t, 1?")T23(0, a?')), we 
can do this by first writing the supergravity action in momentum space, treating it 

functional of Fourier modes for 0(r, t) where: 



(r,t) = (j){r,t)(f){t) = I due' 



(3.166) 



where as before, we defined the Fourier transform using the curved space time r = 
\/ g{rc) t and not simply t. Although this definition is precise for the theory at 
the cut-off r = only, we will use it also for r < because in the end we will 
only provide description at the boundary (i.e r — )■ oo) where the results would be 
independent of the choice of the cut-off. 

The way we proceed now is the following^'^. We consider the metric fluctuation 
as in (|3.165| ) and plug this in the effective action ( |3.92|) but with c„ = 6„ = 
(we will modify this soon). We will also take a convention where g{rc) = 1 in the 
subsequent analysis to avoid clutter. In the final result we will substitute the exact 

(2) 

value of g{rc). Finally, we will call this resulting action as Sq^q where the subscript 
(2) involves writing the action in terms of quadratic 0(r, w). The reason for doing 



^•^Note that the above choice of the warp factor means that we have instinctively chosen the 
minima where 9i = 02 = tt. Needless to say, since we are analysing the effect of the flavor on the 
viscosity, we took the point close to where the quark string originally was. Again a more detailed 
analysis could be performed directly from the ten-dimensional point of view, but such a analysis do 
not reveal any new physics. 

^''This is almost similar to the procedure of [Q whom we refer the readers for more details. Notice 
however that the theory considered by |^ has no running but contains higher curvature-squared 
corrections, as mentioned above. 
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this is because there exists a very useful relation for computing the shear viscosity 
(see for example [|65l): 



lim Im Glf(u,t] 



2T 

lim — Im G^itu, ) 



(3.167) 



where Gf^ is the Schwinger-Keldysh propagator [|65l , |66| . Comparing this with our 
earlier Kubo formula ( p.l54| ) we see that the shear viscosity is nothing but: 



21 <^^o 

Thus if we can write our effective supergravity action in the following way: 



(3.168) 



(2)r 
SGI 



{ro,uj)] 



(3.169) 



we 



where ro is a specified point, then taking the G^f (w, ) part and using (|3.168 
can easily get our shear viscosity^^. Saying it a little differently, we will be taking 

(2) 

two functional derivatives of S'gQ[0(ro, w)] with respect to (^(ro,a;) and thus are only 
interested in terms quadratic in 0(ro, u) in the action. Of course as mentioned above, 
in real time formalism, we are concerned with the Schwinger-Keldysh propagator 
Gfj^ of the doublet fields (f)i{r,t), (j)j{r,t). In the context of gauge/gravity duality, 
we follow the procedure outhned by for AdS/CFT correspondence and treat 
(j)i{r,t), (f)2{r,t) as the perturbation 4>{r,t) and it's doublet in the four dimensional 
Minkowski space^^. In ten dimensional gravity theory, 0i(r) = (f){r) is the field in 
the R quadrant of the Penrose diagram while 02('") is the field in the L quadrant. 
For more details see m J7| ]7ll |7| |73| and |68|. 

Let us make this a bit more precise. Our aim is to get the effective action in 
the form (|3.169|) . To this effect we take our metric (|3.165|) and plug it in the five 
dimensional effective action ( p.92|) . The net result is the following action: 



q 



(2) 
SG 



dud^q 



dr 



rh 



yl(r)0(r, —u)(p"{r, u) + B{r)(f)'{r, —u)(f)'{r, u) 



+ C{r)(f){r, —uj)(f)'{r, u) + D{r)(f){r, —ijj)(t>{r, u) 



(3.170) 



^^Notice that there would be an overah volume factor of T^'^ that would appear with the effective 
action. This factor is harmless and just modifies the Newton's constant in five dimensions. 

^^We are a little sloppy here. Our background is a deformation of the AdS space and therefore 
one might not expect the arguments of [ |66| to carry over exactly as above. However although 
the gauge/gravity duality argument for our case is more involved, we can still consider the 0i 
perturbations to compute the Schwinger-Keldysh propagator because by definition a propagator 



always appears sandwiched between the fields. See also |68 for more generic approach 
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where prime denotes derivative with respect to r and the exphcit expressions for 
A, B, C, D are given in Appendix E. The five dimensional Newton's constant is 
given by: 

^ (3.171) 

where volume of T^'^ i.e Vtli is dimensionful and Kiq is proportional to ten dimen- 
sional Newton's constant. 

The fiuctuation 0(r, u) is not anything arbitrary of course. It satisfies the fol- 
lowing Euler- Lagrange equation of motion: 

,/// X A'(r)-B'(r) , 2D{r) - C ir) ^ A' ir) , ^ , ^ 

which we can derive from ( ^.170 ) above. Once we plug in the values of A, B etc., the 
above Euler-Lagrange equation takes the following form: 



_g{r) r 



h{r) = 



1 + 



2'kN 



1 + 



(r, w) + 



9[r) 



+ J{r) 



[r,uj) 







27r 



logr + 



Air 



logr 



(3.173) 



where i7(r) and Ai are added to allow for most generic conditions on the scalar 
fields^^. As before, primes in (|3.172|) and (|3.173|) denote derivatives wrt the five 
dimensional radial coordinate. 

Now as we mentioned above in ( p.l66|) , 0(r, cj) can be decomposed in terms of 
0(r, |c<;|) and 0(a;). Then as a trial solution, just like in 0, we first try 0(r, u) = g{r)"' 
and look at ( |3.173| ) for r near rh where g{r) — > 0. Plugging this in (p. 1731) with 
g{r) = we obtain the following expression for 7: 



7 



±i 



h{rh)g{r^ 
16 



-rh 



(3.174) 



where in the last step we have used the definition of temperature Tc as in ( p.l8| ). 

To get the solution with g{r) 7^ we propose the following ansatze for the 
solution to ( pTf3| )^^: 



r, w 



F{r,\uMco) 



(3.175) 



^^In special cases we expect J7'(r) and Ai to vanish (see for example |£^). However when this 
is not the case (as may arise for non-trivial UV completions of our model) we could expect a 
non-minimally coupled scalar field. 

^^Incidentally, this form for (j){r,uj) can be shown to be exactly like (3.125) discussed earlier. This 
will become clearer as we go along. 
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Plugging this in (|3.173|) we see that the equation satisfied by F(r, |a;|) can be ex- 
pressed in terms of 7 and 7^ in the following way: 



F'\rM) + 
2g'{r) 



+ 7< 



g{r) 



g'jr) ^ 5 
gir) r 

F\r,\co\) + 



M F\r, \uj\) + 



+ \ - + M 



9"{r) 
. 9{r) 



\uj\^g{rc)h 
g'{r) 



+ J{r))F{r,\u\) (3.176) 



9[r) 



F{r,\u;\)}+^'^-^F{r,\uj\) = 



where the 7^ terms come from both the last term in the above equation as well 
as the Icup term above. Furthermore, note that the source J7'(r) ~ 0{gs) + O^gf), 
so in the limit gs ^ we find that (|3.176|) has a solution of the form F{r, \u\) = 
Ci + C2g{r)~'^"' with Ci, C2 constants. Then we expect the complete solution for gs ^ 
to be F(r, \uj\) = Ci + C2g{r)~'^'^ + f {r, Demanding that F{r, \co\) be regular at the 
horizon r = Vh forces C2 = as g{rh) = 0. We choose ci = 1 and / = ^+77^+7^/0+... 
as a series solution in 7. Then our ansatze for the solution to (|3.176|) becomes 



F(r, \co\) = 1 + g{r) + -fn{r) + -f^}C{r) + 



(3.177) 



Once we plug in the ansatze (|3.177|) in (|3.176| ) we see that the resulting equation can 
be expressed as a series in 7: 



g" + 



- + - + M]G' + jii + g) 

g r 



+ 7{n" + { - + ^ + M ] n' + jn + '^g' h 



\g 

— + 


(-- 




g 




J g. 



-,/2 



+ K0 + 



+ 7^ 



g 



+ 



M 



-,/2 



U + 



+ 0(7^ 



(3.178) 



where we have avoided showing the explicit r dependences of the various parameters 
to avoid clutter. We have also defined in terms of the variables of ( p.l74| ) in the 
following way: 



Kq = - 



16 
TV 



(3.179) 



Although the above equation ( p.l78| ) may look formidable there is one immediate 
simplification that could be imposed, namely, putting the coefficients of 7^,7,7^, ... 
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individually to zero. This is possible because one can view 7 to be an arbitrary 
parameter that can be tuned by choosing the graviton energy u or the temperature 
Tc. This means that the zeroth order in 7 we will have the following equation: 



g"{r) 



g{r) r 



g'{r) + J{r)[l + g{r)] = 



(3.180) 



In the above equation observe that the source J{r), for cases where it is non-zero, has 
a complicated structure with logarithms and powers of r. To simplify the subsequent 
expressions, let us choose to work near the cut-off r = r^.. This is similar to the spirit 
of the previous section where we eventually analysed the system from the boundary 
point of view. Then to solve ( ^.18U| ) near r ~ Tc we can switch to following coordinate 
system 



r = rc(l - C) 



(3.181) 



Taylor expanding all the terms J^{r),g{r), in (|3.180|) about C = and using 

similar arguments as the previous section, we obtain a power series solution for Q as: 



g{r) 



00 

EE 

a i=0 



a) 



c{a) 



;i-c) 



4i 



i=0 



(3.182) 



Since ( p.l80| ) is a second order differential equation, we can fix two coefficients and 
we choose = 0. Then the rest of a[s are determined by equating coefficients 

of on both sides of equation (|3.180|) . The exact solutions are listed in Appendix 
E. Note that all are proportional to Qs and in the limit Qs ^ 0, Q ^ 0. 

To next order in 7 we have an equation for "H that also depends on the solution 
that we got for Q. The equation for 7i{r) can be taken from ( |3.178|) as: 



n"(r) + 
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+ M{r) 


n'{r) 
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^ + - 




[ai 


r) 


r 
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'5 


+ M{r) 


9'{r) 1 


I gir) 




r 




9{r) J 



2g\r] 
9{r) 



G'i 



(3.183) 



To solve this we make the coordinate transformation ( p.l81|) and plug in the series 
solution for Q{r) given above. The final result for T-L can again be expressed as a 
series solution in ( in the following way: 



•H(r) 



00 

EE 

a i=0 



lie) 



C) 



(3.184) 



We again set bo = bi = and following similar ideas used to solve for Q, we determine 
all b'j^s by equating coefficients in( p.l83| ). The exact solution is given in Appendix 
E. Again note that all bi are of at least 0{gs) and thus with 0, H —> 0. 
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Finally the second order in 7 is a much more involved equation that uses results 
of the previous two equations to determine fC. This is given by: 



/C"(r) + 



9'{r) 
9{r) 



9"{r) 
9{r) 



r J 

g'{r) 



+ -+M(r] 



9{r) 



Uir) 



1^0 + 



___2g\r) 

9{r) 
9'\r) 



n'{r) (3.185) 



9'^{r) 



which could also be solved using another series expansion in (we haven't attempted 
it here) . Therefore combining (|3.182|) and ( p.l84| ) we finally have the solution for the 
metric perturbation: 

|2 



0(r, |a;|) 



l + g{r)±i- 



\UJ\ 



-nir) - 



\UJ\ 



,/C(r) + .... 



(3.186) 



We can analyse this in the regime where the gravitons have very small energy, i.e 
— 7- or equivalently 7—7-0. In this limit we can Taylor expand 0(r, |a;|) about 
7 = to give us the two possible solutions: 

\uj\ 



[rM)± = l + Q{r)±i 



A-nT,. 



nir) + [1 + ^(r)]log g{r] 



\UJ\ 



/C(r) + nir) log g{r) + [1 + g{r)]\og^g{r] 



(3.187) 

O{\co\') 



which consequently means that to the first order in u the off diagonal gravitational 
perturbation at low energy is given by two possible solutions corresponding to positive 
and negative frequencies as: 



[1 + Qir)] ^(oj) ± i^{nr) + [1 + Gir)]\og gir)]^{u) (3.188) 



As is well known following, say, [74| [36 , we can define field on the right R and left 
L quadrant of the Kruskal plane in terms of 0+(r, w) and 0_(r, w) in the following 
way: 



(f)±{uj,r) in R 
in L 
(f)±{u,r) in L 
in R 



(3.189) 



Now 0R,±,0L,± contain positive and negative frequency modes but a certain linear 
combination of 0r,±,0l,± gives purely positive or purely negative frequency modes 
in the entire Kruskal plane |^ ||66|. Furthermore imposing that positive frequency 
modes are infalling at the horizon in R quadrant and negative frequency modes are 
outgoing at the horizon in R fixes two combinations : 



-■pos 

0ncg = 0R,+ (w,r) + e'^/^'^'=0L,+ (c^,'^) 



(u;,r) + e'^/2^>L 



0L-(w,r) 



(3.190) 
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With ( 3.190 ) we see that we can define fields in R(L) quadrant as hnear combination 
of positive and negative frequency modes 

(l>n.{(^,r) = ao[0R,+ (a;,r) - e'^^^'^^R _(a;, r)] = 



1 

>2 (3.191) 
2) defined on the complex 



r) = aoe-' — [0l,+(w, r) - 0l 

where we have identified 0r(0l) with the thermal field 0i 
time contour which familiarly appears in the Schwinger-Keldysh propagators of real 
time thermal field theory. Here is a constant. The final physical quantity that we 
will extract from here will only depend on T, as we will show soon. 

Having got the graviton fiuctuations 0(r, w) = 0R(co',r), we are almost there 
to compute the viscosity rj using ( p.l68|) . Our next step would be to compute the 
Schwinger-Keldysh propag ator ^^^(O,^). All we now need is to write the action 
( P.170|) as ( p.l69|) and from there extract the Schwinger-Keldysh propagator. This 
analysis is almost similar to the one that we did in the previous section, so we could 
be brief (see also [^). The action ( p.lTOj ) can be used to get the boundary action 
once we shift 0(r, uj) to 0(r, uj) + 50(r, u) in the following way: 



5, 



(2), 



+ 



9{rc) 



-1/2 



SttG 



N 



dud^q 



rh 



dr<^ A{r)(f){r, —u)(j)"{r, u) + B{r)(f)'{r, —u)(f)'{r, u) 

(3.192) 



-C(r)0(r, —u)(j)'{r, u) + D{r)(j){r, —ijj)(j){r, u) + 2A(r)0"(r, u 
-2fi(r)0"(r, u) - 25'(r)0'(r, co) - C\r)(f){r, u) + 2D{r)^{r, u[ 



-A"(r)0(r, to) + 2A'(r)0'(r, to) 



-U) + dr 



2B{r)(f)'{r,uj)6(j){r,-Lo) 
+C{r)(f){r, u)6(f){r, —u) + A{r)(f){r, u)6(f)'{r, —u) — dr (A(r)0(r, u)) S(f){r, —u) 

Plugging in the background value of 0(r, u) will tell us that only the boundary term 
survives. And as before, to cancel the A{r)(f){r,u)6(f)'{r, —u) we will have to add the 
Gibbons-Hawking term to the action The net result is the following boundary 
action: 



(2) 
SG 



9 [re 



-1/2 



SttGn 
1 



dud^q 

dud^q 



r, -cj 



C{r)-A'{r) + B{r)-A{r) 



[r, -u) 



r, -u 



■0(r,w) 



rh 



?>TiGN\/g{rc) 



f27r) 



(3.193) 



rh 



Now comparing ( p.l69| ) with ( p.l93| ) we see that the terms between the braces com- 
bine to give us the required Schwinger-Keldysh propagator: 



Gff(0,t) 



lim 

47rGArA/fi'(rc) 0i(r,a;)</)i(r, -w) 

lim = 

-^^0 47rG'jvvM 



(3.194) 



rh 




C{r)-A'{r) + B{r-)-A{r) 



0i(r, -w) 



rh 
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where we assume^^ that 0i(r/i,a;) = ao [0R,+ (r/i, w) — e'^/"^'=0R^_(r/j, w)] . Now to 
evaluate the shear viscosity from the above result we need to perform two more 
steps: 

• Evaluate the contributions from the UV cap that we attach from r = Tc to r = oo. 

• Take the imaginary part of the resulting total Schwinger-Keldysh propagator. This 
should give us result independent of the cut-off. 

To evaluate the first step i.e contributions from the UV cap, we need to see precisely 
the singularity structure of S^q. The second step would then be to extract the 
imaginary part of SK propagator from there. Since the imaginary part can only 
come from the second term of ( p.l93| ), we only need to evaluate: 



lim 



1 



A-kGn^M^ 



B{r) - A{r) 



(r, -uj) 



(3.195) 



with 0(r, —uj) being the graviton fluctuation in the regime r > r^. To analyse this 
let us first consider a case where (^s — > and (^(r), 'H(r), /C(r), ...) — )■ 0. In this limit 
we expect for Vh < r < Vc- 



B{r) - A{r) = 
0i(r,a;) = do 



^9. 



g{ry + 0{g,Nf) 



(3.196) 



271 



<pi[r,-uj) g[r) \A-n^ + \og^g[r 

The above considerations would mean that the contribution to the viscosity, r^i, for 
this simple case without incorporating the UV cap will be: 



1 




3 r 2 



1 



'^qIGn U7r + Mog2^(re 



(3.197) 



where we have used the relations TTTc^/g(j 



rh) 



and h{rh] 



^ ^ in this 



limit. This helps us to write everything in terms of T and not the scale dependent 
temperature T^.. In fact as we show below, once we incorporate the contributions 
from the UV cap, the Tc dependence of the above formula will also go away and 
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At this point one might worry that the solution for (jji is only known around rc. That this is 



not the case can be seen in the following way: Integration by parts gives (3.194) which says one 
only needs to know the value of the field (pi at Vc and r^. The solution for 0i = ^ is given in 
( 3.186 ) from which it is clear that (j>i{rh) = as g{rh) — 0. Furthermore to know rj we only need to 
know the imaginary part of (3.194), which is evaluated using ( |3l96| ) in ( |3l9^ ) and using boundary 
values of 0i(rc) and (f>i{rfi). 
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the final result will be completely independent of the cut-off. Note that in the limit 
Tc — )■ oo we recover the result for the cascading theory. 

Combining all the ingredients together, the contribution to the viscosity in the 
limit where (^(r), 'H(r), /C(r), ...) etc are non-zero can now be presented succinctly 
as (although rji below doesn't have any real meaning on the gauge theory side as this 
is an intermediate quantity): 



1 + 



4:ri 



H' 

1+g 



2glG^ 1 47^ + i[log^7(r,) + ^]■ 



(3.198) 



Note that the above expression is exact for our background at least in the limit where 
we take the leading order singularity of the background. This is motivated from 
our detailed discussion that we gave in the previous section. Note that the second 



term in the action ( |3.193| ) is exactly the second equation of the set (|3.124|) whose 
singularity structure has been shown to be renormalisable. Thus taking the leading 
order singularity instead of the actual r^^^ will not change anything if we carefully 
compensate the coefficients with appropriate gsNf,gsM'^/N factors! 

But this is still not the complete expression as we haven't added the contributions 
from the UV cap. Before we do that, we want to re-address the singularity structure 
of the above expression. The worrisome aspect is the existence of factor in (|3.198|) . 
Does that create a problem for our case? 

The answer turns out to be miraculously no, because of the form of H and Q 



given in (|3.184|) and ( p.l82|) . This, taking only the leading powers of Vc, yields: 



n' 



862 



+ 



G' 



4ai 



802 



+ .... 



(3.199) 



killing the dependence in ( p.l98| ) . This would make r]i completely finite and all 
the Tc dependences would go as 0(l/rc). Therefore we expect the contribution to 
the viscosity from the UV cap to go like: 



00 ^ 

i=0 c 



(3.200) 



where the total viscosity will be defined as r/ = r^i + 772- As this is a physical quantity 
we expect it to be independent of the scale. Therefore 



dr] 







(3.201) 



70 



It is now easy to see why bi 



additional log r terms from 



(a)- 



cii = is consistent. For non-zero 61, ai there would have been 
These would have made the theory non-renormalisable. Thus 



holographic renormalisability would demand hi = ai 
with what we choose earlier. 



from the very beginning — consistent 
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which will give us similar Callan-Symanzik type equations, as discussed in the previ- 
ous section, from where we could derive the precise forms for Gi in ( p.200|) . Finally 
when the dust settles, the result for shear viscosity can be expressed as: 



1] 



N 



1 + E 



oo 



-4fcA/'uv 



47r + i log' (1 - r%- 



(3.202) 



where ak are functions of T that can be easily determined from the coefficients 
{di,bi) in (|3.182|) and ( p. 184 ) or {ai,bi) worked out in Appendix E; and h{T) = 
+ 0{gs, Nf, M). Observe that the final result for shear viscosity is completely 
independent of Vc and T^, and only depend on T and the degrees of freedom at the 
UV i.e through e"-^"". Needless to say, for large enough Afu^ (which is always the case 
for our case because Muv — ^ e~"',n > 1), the shear viscosity is only sensitive to the 
characteristic temperature T of the cascading theory. The interesting thing however 
is that the shear viscosity with finite but large enough Afuv can be smaller than or 
equal to the shear viscosity with Afuv e~"',n » 1 i.e for the parent cascading 
theory provided: 

^ n-Ak 

~ — ^' n<k, keZ (3.203) 



< 



47J-2 



n] 



in the limit of small characteristic temperature T. This will have effect on the 
viscosity to entropy ratio, to which we turn next. 

3.5 Viscosity to entropy ratio 

Our final set of analysis will be to calculate the viscosity to the entropy ratio for the 
above two cases i.e one with only RG fiow, and the other with both RG fiow and 
curvature squared corrections. As usual the former is easier to handle so we discuss 
this first. 

Starting with the type IIB supergravity action in ten dimension i.e the S'oks-bh 
in ( p:B71 ) the entropy is given by the Wald's formula J73| , |7§ , [[77) , |7S[ 



S 



-2ti i dxdydzd^MW 

J oRabcd 



(3.204) 



where the integral is over the eight dimensional surface of the horizon at r = Vh, Ciq 
is the lagrangian density of the action in ( |3.67| ), Vab, a,b = 1..8 is the induced 8x8 
metric at horizon, e^fe is the binormal normalized to tab^""^ = —2. Finally using 
explicit expression for the metric (|3.4|) and ( p.5|) , we have 
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(3.205) 
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where V3 is the infinte three dimensional volume and we have used the definition 
of five dimensional Newton's constant Gn introduced in ( p.l71| ) as well as h{rh) 
introduced in ( |3.173| ). The relation above is consistent with the notion that the 



effective five dimensional warp factor is well approximated with taking a slice 9i = 
Ti^ipi = ip = 0. The results would only differ by a 0{gsNf) term which, in our 
approximation, is very small. Note that the definiton of temperature depends on 
the effective five dimensional warp factor and as the approximation in (|3.2(J5|) is 



consistent with it, our computation of entropy and temperature are consistent. 

Once we replace by the characteristic temperature T, we see that the entropy 
is only sensitive to the temperature and is independent of any other scale of the 
theory. Since the above result is also independent of J\fuv it would seem that the 
Wald formula only gives the entropy for the theory with M'uv = 00 i.e for the parent 
cascading theory^^. The interesting question now would be to ask what is the entropy 
for the theory whose UV description is different from the parent cascading theory? 
In other words, what is the effect of the UV cap attached at r = on the entropy? 

To evaluate this, observe first that in finite temperature gauge theory, entropy 
density of a thermalized medium having stress tensor {T'^'^) = diagonal(e, P, P, P) is 
given by 

e + P , , 

s = -jr- (3.206) 

where e is the energy density, P = Px = Py = Pz is the pressure of the medium and 
T being the temperature. With our gravity dual we can compute the stress tensor 
(-^med) {^^^ thus the energy e = (T^ed) pressure P = (T^ed)) medium 

through equation of the form (|3.66|) , i.e 

(Cod) = (3-207) 

Obgpq 

where again p,q = 0,1,2,3 and gpg is the four dimensional metric obtained from 
the ten dimensional OKS-BH metric gij,i,j = 0,1.., 9; and 6b operation has been 
defined earlier. There are two ways by which we could get a four-dimensional metric 
from the corresponding ten-dimensional one. The first way is to integrate out the 
6i, (pi directions to get the four-dimensional effective theory. This is because the warp 
factor for our case is dependent on the di directions. The second way is to work on a 
slice in the internal space. The slice is coordinated by choosing some specific values 
for the internal angular coordinates. Such a choice is of course ambiguous, and we 
can only rely on it if the physical quantities that we want to extract from our theory is 
not very sensitive to the choice of the slice. Clearly the first way is much more robust 
but unfortunately not very easy to implement. We will therefore follow the second 



^^This can be argued by observing that fact that in a renormahsable theory, hke ours, the de- 
pendences on degrees of freedom go hke O (e~^"") corrections as we saw in the previous sections. 
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way by choosing the the five dimensional shce as 6i = 62 = 7r,ip = = (j)2 = and 
thus obtaining 



g^.u{0^ = 7l,i^ = (pi = O) (3.208) 



with /i, z/ = 0,1,2,3,4. The next step would be to evaluate all the fiuxes and the 
axio-dilaton on the slice. To do this we define: 



IH3P = 


\H,\\9, = n,^ = 


<P^ = 0); 


IF3P = 
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k = 0); 


|Fir = 


l^il 



2/ 



2/ 



TT,^ = = 0) 



= vr, = 0j = 0) 

$ = $(^. =n,ip = ^i = 0) (3.209) 

Once the fiuxes have been defined, we need the description for Stotai in ( p.208| ). This 
is easily obtained from ( p.67|) as: 



24 



d'x./^\ |Fi|2 + IF3P + llF^l'' 1 (3.210) 



with R being the Ricci-scalar for g^^, and g = det gf^^. Note that in the definition 
for the slice sources H3, Fi, F3, F5 and R, we still have gij, i,j > 5 which we evaluate 
at 9i = 7i,ip = = 0, treating them simply as functions and not metric degrees of 
freedom. 

To complete the background we need the line element. Here we will encounter 
some subtleties regarding the choice of the black-hole factors and the corresponding 
QsNf type corrections to them. With the definition of g^^^ the line element is: 



ds' = - -^M=de + V^E^I^rfr^ + - d-^' (3.211) 



^Jh(r, TT, tt) g2{r) ^/h{r~K,TT] 



i,j=0 j,a (a) 



i,j=0 j,a (a) 

where aij, Pij are all of 0{gsNf, QsM) and only involve the parameters of the theory 
namely, rh, L and fi from the embedding equation ( |3.49| ) and guarantees that ^ 
are dimensionless. On the other hand Cj""*, Kj""* can incorporate zeroth orders in ggNf. 
However note that so far we have been assuming gi{r) ^ 92(1^) = g{r), ignoring their 
inherent 6i dependences, and also the inequality stemming from the choices of aij 
and Pij. This will be crucial in what follows, so we will try to keep the black hole 
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factors unequal. These considerations do not change any of our previous results of 
course. 

Now looking at the form of the metric, knowing the warp factor /i(r, vr, vr) and 
gi{r), just like before we can expand the line element as AdS^ line element plus 
0{gsNf, QsM) corrections. We can then rewrite the line element ( |3.211| ) as: 

ds' = - 4 [g{r) + h] de + VSw) + ;!(! + 

g2[r) U 

i,i=0 



oo 



log^ 



h{r) = ^Q^^^ (3.212) 



ij=0 

where again 7jj, C,ij are of 0{gsNf, gsM) and we are taking /i(r, vr, vr) = ^+0{gsN f, gsM). 
Such a way of writing the local line element tells us that there are two induced four- 
dimensional metrics at any point r along the radial direction: 

g^J) = diagonal (-5- (r), 1, 1, 1), = diagonal (-/i, /s, /2, ^2) (3.213) 

where we haven't shown the r^/L^ dependences. The reason for specifically isolating 
the four- dimensional part is to show that we can study the system from boundary 
point of view where the dynamics will be governed by our choice of the boundary 
degrees of freedom. It should also be clear, from four- dimensional point of view, the 
metric choice g^pj is directly related to the AdS geometry whereas the other choice 
g^pq is the deformation due to extra fluxes and seven branes. This decomposition is 
similar to the decomposition that we studied earlier. 

The above decomposition also has the effect of simplifying our calculations of 
the energy momentum tensor (T^^^). We can rewrite the total energy momentum 
tensor as the sum of two parts, one coming from the AdS space and the other coming 
from the deformations, in the following way: 

/T^pg \ = '^6^^°^ total 

= (^mcd)AdS + (^med)def 



'total 



Sll, + SlL (3.214) 

where S|°'^^[ is zeroth order in gsNf,gsM and Sl^^g^j is higher order in gsNf,gsM. 

Note that {T!^j.)Ads = r ^'f^ is the well known AdS/CFT result obtained from the 
analysis of [|l3l I3[|l^[|l6] in the limit Vc — t- 00. With the 0{l/r) series expansion 
of our metric gQQ = 1 — r'^/r^, gn = g22 = §33 = 1; the result at the boundary is 

\ med/Ads 2^2c^ 2g^Gr, 

(^red)Ads = m,n = l,2,3 (3.215) 
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This only gives the CFT stress tensor as we evaluate the tensor on the AdS bound- 
ary at infinity, reproducing the expected first term of ( p.205|) . How do we then 
evaluate the 0{gsNf, QsM) contributions from the deformed AdS part i.e the energy 
momentum tensor (T^ed)def at any r = cut-off in the geometry? 

In fact the procedure to evaluate exactly such a result has already been discussed 
in the last two sections: namely the wake analysis in section 3.3 and shear viscosity 
analysis in section 3.4. Therefore without going into any details, the final answer 
after integrating by parts, adding appropriate Gibbons-Hawkings terms and then 
using the equation of motion for gpj, we have 

1 r d^q 



'total 



Crir, q) - (r, q) <|.W(r, g)$W(r, -q) 



$:i^l(r, g)$W(r, -q) + <|.W(r, g)<|.'W(r, -q) 



+ {E^-FrY^^{T,q) 



(3.216) 



where rmax = ^c{^ ~C)- The values of the coefficients are given in Appendix F. The 
above form is exactly as we had before, and so all we now need is to get the mode 



expansion for $m . Note however that the subscript m can take only two values, 
namely m = 0,1 as there are only two distinct fields ggg and gfj = = 833- 
Therefore our proposed mode expansion is: 

$W = gi^J = > > 

m Oram / j / j 



00 



a 1=0 c{«)V 



(3.217) 



Just like the wake and shear viscosity analysis, the action in ( |3.216| ) is divergent 
due to terms of 0{r'^),0{r^) and hence we need to renormalise the action. The 

equations for renormalisation are identical to the set of equations (|3.126|) (|3.131|) , 

and therefore we analogously subtract the counter terms to obtain the following 
renormalised action: 



d^q 



1 - 



rtil - 0^, 




Qmn[l]^ $ 



(a) 

mn(i)[l] 



, i=0 



r 



(a) 



mn[l] 
\a\ 



_4g(4)[/3]$ _4s(4)[/3]$ is(5)[/3]$/ , ^(S)!/?] ^/ 1 , 



^m(i)[l] 



i=0 ^{") 



(3.218) 

where the radial coordinate is measured at the two boundaries Vh and rc(l — C) and 
$m are independent of r as before. Note that X[r(Q,)] is a function independent of 
and appears for generic renormalised action. 
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Now the generic form for the energy momentum tensor is evident from looking 
at the hnear terms in the above action ( p.218| ). This is again the same as before. 
However now we also need the entropy from the energy-momentum tensor as in 
( p. 206 ). The result for the energy-momentum tensor at r = Vc is given by: 



med /def 



d^q 1 ^ 



,TTmn[l\ rT"ni[lK„(4)[/3] 
^-"|a| |a| '^rm 



^nm[lK (4)[/3] 



, 1=0 



(3.219) 



mn[l] T^mn[l] ,(a) 

'-^^In-I ' n(i)[l] 



The explicit expressions for the coefficients listed above, namely, H^^^ , 

and Snn^^ are given in Appendix F. 

To complete the story we need the contribution from the UV cap. This is similar 
to our earlier results. The final expression for the ratio of the energy-momentum 
tensor to the temperature takes the simple form: 



med /def 



8nG 



N 



a,(3 



(4)[/3] _ /I r 



I /^'^'»W^„(4)[/3] 



/^mn[l] TV-rimfllx (5)[/3] 



-r "n(j)[l]"nm>- 



-jMuv 



(3.220) 



We would like to make a few comments here: First, observe that the final result is 
independent of our choice of cut-off. Secondly, in the string frame there should be a 
l/gl dependence. Finally, we can pull out a term because the coefficients have 
an explicit dependences (see Appendix F). This means that both from the AdS 
and the deformed calculations performed above we can show that the entropy is of 
the form: 



(3.221) 



where the first part is from ( 3.215 ) and the second part is from ( p.220| ). The result 
for the parent cascading theory is ( p.205| ), and so we should regard ( p.221| ) as the 
entropy for the theory with Afuv degrees of freedom at the boundary. Of course in 
the limit Afuv e~",n >> 1 we should recover the entropy formula ( 3.205 ) for the 
parent theory. All in all we see that the correction due to JV^v degrees of freedom 
only goes as e~^^^, so in practice this is always small for the type of Afuv that we 
consider here. This means that we can use the entropy for the parent cascading 
theory to estimate the viscosity by entropy ratio for a system with Afuv degrees of 
freedom at the UV as: 



s 



1 + Er=i ake-''"^- 
47r + i log^ (1 - r%-^-) 



(3.222) 
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where we see that the boundary entropy term ( 3.205|) neatly cancels the coefficient 



in the viscosity ( p.202| ) to give us the precise bound of ^ when A/'uv — )■ oo. Of course 
from our other analysis ( |3.221| ) we might expect a O {^ggNf, ggM, e"-^"") contribution 
that would make (|3.222|) saturate the celebrated bound j- if the total entropy density 



factors compensate the factors coming from the viscosity. This would seem consistent 
with, for example, |]7U|^^. In fact our conjecture would be for non-zero M,Nf and 



Afuv — ^ €'"^,11 » 1, the bound is exactly saturated''^ i.e 7 = 4^- 

Our second and final step would be to incorporate both the RG flow as well as 
curvature square corrections. As we discussed before the curvature squared correc- 
tions are typically of the form c^R^upaR^"^'^ with C3 being the coefficient ( |3.163D that 
we computed before. 

The crucial point here is that (see where this has also been recently empha- 
sised) in the presence of curvature squared corrections the five dimensional metric 
itself changes to: 

ds^ = ^ dt^ + ^ \' ' ' dr^ + -^=^= (3.223) 

^Jh{r,'K,'K) g2\f) A//?.(r, 7r,7r) 

where the black hole factors gi are no longer given by (|3.13|) or its simplified version 
( p.86|) . They take the following forms: 

4 8 1 m 

, , ^ r? r? _ log r 

^72(r) = l-^ + a + 7^ + /3.„^ (3.224) 

where are all of 0{gsM, gsNj) and can be worked out with some effort (we 

will not derive their explicit forms here). Similarly we could also express ( ^.2241 ) in 
terms of inverse powers of r to have good asymptotic behavior. Observe that we can 
still impose gi ~ g2 because the corrections are to 0{gsNf, ggM), although all our 
previous analysis have to be changed in the presence of curvature corrections because 
the explicit values of gi{r) have changed. We will address these issues in the sequel 
Finally («,7) are given by 

« = 7 = ^ (3.225) 

At this point one might get worried that the metric perturbation on this background 
would become very complicated. On the contrary our analysis becomes rather simple 



^^Provided of course if we assume that afc's are more general now, being functions of 7~, gsM, gsNf. 
This way even for non-zero M, Nj, whenever we have A/'„„ — > 00 the bound is exactly 

^•'Note that any possible deviations from ^ due to ( 3.19E ) in ( ^.198 ) cannot happen because 



the underlying holographic renormalisability will make di = bi — 0, as discussed earlier. Thus the 
bound in itself is a rather strong result. 
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once we ignore terms of 0{c3gsM, c^ggNf) (which is a vahd approximation with 
C3 << 1). In this hmit the metric perturbation can be written simply as a linear 
combination of the terms proportional to C3 in $ which appears in ^ and our solution 
( p.l86| , p.l87| ) derived for RG flow. The final result is: 



l±i- 



\UJ\ 



nir) + [l + g{r)]\og g{r) + 



r^g{r) 

/C(r) + H{r) log g{r) + [1 + Q{r)]\o^g{r 



a + 47^ 



167r2T2 

+ + 0{c:,gsNf) + 0{o,gsM) 



(3.226) 



where we have written ( p.l3|) as; 



g{r) + 0{c^g,Nf) + 0{c^gM) 
g{r) + 0(c3(?,iV/) + 0{c^gM) 



(3.227) 



with g{r) being the usual black hole factor defined in (|3.86 ). Of course as emphasised 
above, this is valid only in the limit C3 << 1, which at least for our background seems 
to be the case (see ( p.l63| )). 

The above corrections are not the only changes. The entropy computed earlier 
also gets corrected and therefore the horizon can no longer be at r = r/^. To evaluate 
the correction to entropy we again ignore the terms of 0{c3gsM, c^ggNf). In this limit 
the correction terms are precisely given by the analysis of [0| and are proportional to 
the C3 factor (|3.163|) as expected. This means that the final result for 77/5 including all 
the ingredients i.e RG flows, Riemann square corrections as well as the contributions 
from the UV caps; is given by: 



7/ 

s 



47r + i log^(l-r%-^-^ 



C3K 



3L2(1 -r^e-^A^-)^/^ 



Bo{AtI^ - log' Co) + ATiAo log Co 

(^47r2-log2 Co)\w log' Co 



(3.228) 



where we see two things: one, the bound is completely independent of the cut-off 
r = Tc in the geometry, and two, the bound decreases in the presence of curvature 
square corrections even when Muv e " with n = 0{1)7'^ The constants appearing 



^■^In p4[ , non-relativistic systems that appear to have no lower bound were constructed. However, 
these are systems which necessarily require large chemical potentials and low temperature. In highly 
relativistic system created at high energy colliders such as the RHIC or the LHC, chemical potentials 
are small and temperature is high. Our discussion here assumes that the system under discussion 
has such properties so that the use of thermodynamic identity e + P = Ts is valid. Hence, our 
discussion here is in no direct conflict with the models constructed in [E4l. 
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in ( p.228| ) are defined as: 



Co = l- T^e-^- 

Ao = -isr^e-^-^"" + (sr^e-^-^"" - 47r^e-^"'') log Co + 2QT^e-^-^ 

+ 24 (1 + T^e-^-^) log Co 
Bo = -887rr*e-*^"" + 487rr^e"^"" + 48 (3.229) 

This is consistent with 0, and the only violation of rj/s may be entirely from the C3 
factor provided the increase in bound from the first term of ( p.228| ) is negligible, as 
we discussed earlier for ( ^.222 ). This means in particular: 



- = - ribCs + OiTe-^-^) (3.230) 
s 47r ^ ' 



where rii, can be extracted from (|3.228| ). In this paper however we will not study the 
subsequent implication of this result, for example whether there exists a causality 
violation in our theory due to the curvature corrections as in . We hope to address 
this in the sequel . 



To see the explicit behavior of rf/s we can plot the function (|3.228[) (see figure 



12). Of course our result is valid for infinitely large UV degrees of freedom i.e 
■^uv e^",n > 1, but we can extrapolate the result to see what properties we get 
for small enough UV degrees of freedom. Incidentally with the geometry cut-off at 
small Tc and an AdS cap attached from r = Tc to r = 00 would start resembling 
standard QCD as we briefiy mentioned earlier. For such a scenario we expect the 
UV degrees of freedom to be smaller than the UV degrees of freedom for the parent 
cascading theory i.e Muv would approach infinity at a smaller rate. 



4. Conclusions 

In this work we have studied the dynamic response of a strongly-coupled, strongly- 
interacting medium to a fast quark. The drag coefficient and the wake resulting from 
the parton-medium interaction was evaluated. The calculations were performed by 
constructing a gravity dual inspired by the Klebanov-Strassler model with D7 branes 
to take into account fundamental quarks. The gauge dual then has a running coupling 
constant, unlike models stemming from a pure AdS geometry, with features close to 
that of QCD. This procedure in fact allows the consideration of a family of gauge 
theories which have a well-defined completion in the UV, beyond a cut-off scale. We 
could also show that physical results were independent of the choice of this cut-off. 
We have applied this model to the calculation of the drag force and the wake left by 
a moving quark in the strongly interacting medium were computed. In addition, we 
have evaluated the ratio of shear viscosity to entropy density - rj/s - and shown the 
violation of the bound conjectured in Ref. [Q, for a range of parameter values. It 
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is important to test and verify the robustness of this hmit, and we view the current 
work as contributing to this effort. More work is needed in order to identify the size 
and extent of this violation, at the moment this violation is only parametric and our 
parameters need to have better defined physical origins. In the end, one may need 
to rely on the empirical identification of key quantities, like transport coefficients for 
example [^. In this regard, the role of heavy ion experiments at RHIC and at the 
LHC can't be overestimated. 
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Figure 12: Plot for tj/s with and without Riemann Square term. The x-axis is defined as 



k 



where A^uv is the UV degrees of freedom and T is the characteristic temperature 
of the cascading theories. For the parent cascading theory fc — )• oo and we see a violation of 
the bound (the solid blue hne). As k decreases (assuming this is possible) the red dashed 
line dips slightly below the I/Att axis, but the solid blue line remains considerable below 
the l/47r axis. For k sufficiently small the bound is not violated. However all the models 
that we studied in this paper can only realise the large k limit. 



Note added: In our recent paper |81|, we have given a concrete example of a UV 
complete theory that consistently proves all the statements that we made here. Our 
new geometry has a UV cap given by an asymptotic AdS-Schwarzchild geometry. The 
IR dynamics of our theory is captured by the OKS-BH geometry. The deformation 
of the OKS-BH geometry at the boundary is indeed captured by an interpolating 
geometry that ties smoothly, in the presence of sources and fluxes, to the asymptotic 
AdS-Schwarzchild geometry. All the relevant details like changes in stress-tensor. 
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entropy etc due to the UV cap can now be calculated. In addition to that we have 
been able to study zero temperature linear confinement as well as high temperature 
quarkonium suppression and melting from our UV complete theory. Our results 
are also consistent with many of the earlier predictions made in the literature using 
completely different techniques than ours [0, Q ^ . 
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A. Back reaction effects in the AdS Black-Hole background: 
A toy example 



The following analysis, although independent of the main calculations in the paper, 
serves as an interesting warm-up example where we can study metric perturbations 
due to fluxes and D7 brane in a controlled AdS background. This will prepare us for 
the analysis of the next section where we study metric perturbations in the non-trivial 
OKS-BH geometry due to fluxes, D7 brane and strings. 

To start-off consider an AdS^ Black Hole (AdS-BH) metric given in the following 
way: 



f{u)dt^ + dx\ + dx\ + dxl + 



du 



2 -, 



{u)dx^dx'' 



[k.l] 



with fiu) = 1 — Uh is black hole horizon, m = is the boundary, and fj, or u 
0, 1, 2, 3, 4. The above metric is asymptotically AdS |TB[ as 



u'g' 



\u=0 



Vij 



(A.2) 



where rjij is the Minkowski metric, with i or j = 0,1,2,3. Note that we have ex- 
pressed the metric in terms of u = ^ coordinate. We could as well express everything 
in terms of r coordinates (as we did in the main text). 

To this background, first let us now add a D7 brane by switching off the black- 
hole factor. The D brane Born Infeld (DBI) action for D7 brane in 5 dimensional 
AdS-BH background is given by []T7 



5, 



DBI 



d^xC{x) = gs J d^x^/g{x) cos^$(a;) ^Jl + g'"'{x)df,^{x)d^^{x) 



(A.3) 



where the metric g^^iy is the metric of the resulting background, i.e g^vi^) ^^^^ 
In the above, gs is coupling constant, and is a scalar field describing the D7 
brane. The Euler-Lagrangian equation for the $(a;) field leads to: 



□$ + 3tan $ 

2 1 



g"'^d^<^dp<^ 



= 



(A.4) 



where □ 



9^ [ V99'"'d. 



. As ( |A.4p is a second order non linear differential equation. 



the solution can be written in the form 17 



i=0 



u 



(A.5) 



- 96 - 



Now considering first order 0{gs) in the perturbation and terms only upto 0{u^), 
we find (see [|l^ for details), 



^(u,lt,t) = cu ^ (A. 6) 

6 

where c is a constant. Using the lagrangian ( |A.3|) and the solution ( |A.6|) , we can find 
the stress energy tensor for the D7 brane: 



1=1 



where in the last line we have expressed this in powers of Qs to make it more generic. 
In fact, as we will see below, the energy momentum tensor can be assumed to come 
from various sources that include the D7 brane as well as background fluxes. 

Once we introduce D7 brane as well as the black-hole factor / we need to modify 
the scenario. In particular there will be non-trivial axio-dilaton background that 
we have been ignoring. This will generally back-react on the geometry and make 
the system non-AdS. To preserve the AdS like configuration we will need more D7 
branes (and also orientifold planes so that we could lift the system to a F-theory 
configuration |Q) as well as other background fluxes. Let us denote the action of 
all the fluxes etc. in five dimensional space as S'fiuxes, the action of the 0-planes as 
5'pianes and the action of the string from the D7 branes to the black-hole as S'strings- 
This means that our total action is: 



Stotal — J d^xy/C R + Sbbi + •S'auxcs + '5'pi 

anes ~l~ 'S'stj-jngs 



(A.8) 



where R is the Ricci scalar for the metric G^^, = gjj^^, + Qgh^ui and S'dbi now include 
the total action of all the D7 branes. We then expect the Euler-Lagrange equation 
for G^y using the action ([A.8|) determines the Einstein tensor: 



<3A^) = E 9l0;:li^ = E 9ir^{x) (A.9) 



i=l 1=1 



where the energy-momentum tensor is determined by varying the actions S'dbi + 
S'fiuxes + S'pianes wrt the backgrouud metric G^i,. Now equating coefficients of gs in 
the above expansion, we have the following equation for h^,^: 

G^^JiKAx)] = Til\x) (A.IO) 

Notice that symmetric tensor, it has fifteen degrees of freedom. We use 

coordinate transformations to fix ten degrees of freedom and treat h^^, as a diagonal 
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matrix. Consider the coordinate transformation x° — )■ a;'" 
transformation 



+ Qse"". Under this 



hf,^ h 



(A.ll) 



which is a gauge transformation. Thus we can fix five components (say /iqi = ^02 = 
^03 = ^04 = = 0) with the above transformations. Just hke in electromagnetism, 
this does not completely fix the gauge. We can add to another set of functions e'° 
which obey □e''^ = and leaves Einstein tensor unchanged. Demanding that e'° also 
leaves the fixed components invariant, we can further fix another five components 
(say /ii3 = = = /i24 = ^34 = 0). Thus we have diagonal h^^. Observe that 
this is not the most generic choice, but will suffice for this toy example provided we 
can have the energy-momentum tensors to be diagonal. We will also assume that the 
fluxes and branes/planes are independent of l^,t and thus the metric perturbations 
h^i, they induce should only depend on u and we should have hu = /i22 = ^33- With 
all these, we have 



U] 















\ 



h{u) 



h{u) 







h{u) 



V 















(A.12) 



l{u)) 



We observe that ( |A.1CI|) is a second order non linear differential equation for h^j^^. 
This will have a general solution: 



m[u) 



^ 00 



irtjU-' + TJijU-' log u 



(A.13) 



where m{u) = k(u), h{u), l{u) and rrij = kj, hj, Ij, rhj = kj, hj, Ij are constants. 

As we mentioned before the above form of the metric ( A.12 ) will make sense if 
the sources can also be made diagonal. In the presence of fluxes, branes, planes and 
non-trivial profile of the strings this doesn't look very difficult to achieve. We will 
therefore assume that the energy-momentum tensors are also diagonal, and the right 
hand side of ( [A.10| ) can be written as a Lorentz series as: 



-(1) ^(1) 

V{-7) ^ V(-6) ^ 



U 



7 



U 



6 



^1^(3) 



(A. 14) 
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Once ( |A.14| ) is made diagonal, we can satisfy ( A.IO ) easily. The explicit expressions 
for the Einstein tensor are: 



u 



1 



12m 



14 



ui - u^)i2u\l - 2m6))/(m) - 8u\l^k{u)\ (A.15) 



B(u) 



u 



+ ( - Au'ul"^ + 28u'ul' - Uu^'ul + 20u''ut) ^ 



+ (20M^nf - Qu\f - 2Au^\l + I2u^''ui - 2u^''') 

+ 88^1^4 + UQu^^i^ - I2u^^u]^ - {l'92u\]'^{ul 

2\ <fk{u) 



du 



u 



B-^{u)]l{u 



2u'\l + 2u\f^'^^^''^ 



+ [40u\l^ - 24u''u]^ - {9Qu\f{ul - uy)/B{u)Y{u) 



) du 



C{u) 



+ 



- Au^ui + 64m^<(-4 - u^)C-\u))h{u) 



(A.16) 



(A.17) 



where we have defined B{u) and C{u) in the following way: 



4\2 



C{u) 



(A.18) 



with ii = 11,22,33 from here on due to our assumed diagonal nature of the 
perturbation. Once the perturbations become non-diagonal the analysis will become 
more involved as we saw in section 3.3. 

By matching coefficients of various powers of u in ( [A.10|) we can solve for h^^, 
order by order. We obtain equations for the constants mj^fhj by plugging in (|A.13| ) in 
( |A.10 ) and using the explicit expressions ( |A.15| ) onwards and the expansion ( |A.14|) . 
Doing so, we observe that it is sufficient to have mj,mj = for j < —6. 

The equations obtained by matching coefficients of i for —6 < j < —3 are: 



4t, 



(1) 

OO(j-l) 



2(-3 - 2j + f)hj + (-6 + 2j)lj - 8kj + 2{f +j- 2)hj + 2/~- 
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4(3 - 2j + f)h, - 6(-3 + + 2(-23 - 2j + f)k, 
- 4(-2 + J + f)h, - 61, + 2(-2 + J + f)k, 

- (-2 + 6j)hj + (6 + 2j)kj - 6hj + 2kj (A.19) 

where because of the form of Q^i, all the above equations are independent of Uh- 
This will not be the case for j > —3. Finally, matching coefficients of ^^^j^ again for 
— 6 < j < —3 gives us: 

2(-3 - j)hj + 2(-3 + j)lj - 8kj = 

-4(2-j>j-6(-3 + j)^~+2(-23-j>j = 

(-2 + 6j>i + (-6 - 2j> = (A.20) 



At 



At 



(1) 

(1) 

440-1) 



For a given —6 < j < —3, we have six equations and six unknowns in ([A.19|) and 



( |A.20p and thus we can solve exactly. For j > —3, as we mentioned earlier, equations 
governing rrij, rhj will in general depend^^ on rrij-p, rhj^p for < p < j + 6. Once we 
solve for mj^p,mj^p, we can obtain exact solutions for all mj,rhj with j > —3. This 
way we have exact solutions for all j. For —6 < j < —3 the solutions are: 

. (3+j) .,(1) I f - - 61j + 5 (,) 

' {2f - 3f - Q2j - 5) (2j3 - 3f - 62j - 5)(-3 + j) ''^^-^ 

+ ....^}:'t (A.21) 



^ (3+j) (1) _ (-5-21j+j^+/) (1) 

' 3(2j3 - 3j2 - 62j - 5) oo(J-i) ((2/ - 3j2 - 62j - 5)(-3 + j)) "^^"'^ 



l-l + 3j) (1) . , 22^ 

+ 3(2j3 _ 3^-2 _ 62j - 5) 4^(^-1) ^ ^ ^ 



4^ 



^ (j2 - 2j - 35) (1) (45 + 68j-30j^-4j3+j-j 

' 3(2j3 - 3j2 - 62j - 5) ((-3 + j)(2j3 - 3j2 - 62j - 5)) ^^(^"i) 



l-2j-15+j^) (1) 
+ 3(2j3 _ 3^-2 _ 62j - 5) ^^(^-1) ^ ^ ^ 



r ^ (2j^ + 9j3 + II f - 21j - 181) (1) 
' 2((2j3 - 3f - 62j - 5)(j2 + 30j + 1)) "o^-^"^) 

{107 f + 324j^ - 1932j3 - 3614j2 + 697j - 1790) (i) 
~ 2((j2 + 30j + l)(-3 + jT(2j3 - 3j2 - 62j - 5)) ^-(^"D 



, (6j^ + - 17 f + 7j + 77) (,) 

^ 2((2j3 - 3j2 - 62j - 5)(j2 + 30j + 1)) ^^(^-i) ^ ' ' 



^^Also on Uh- This can be easily checked by plugging in the mode expansion for mj,rhj in Q^^ 
jiven above. 
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J _ (2/ + 9/ + llj'-21j- 181) (1) 
' ~ 6((2j3 - 3j2 - 62j - 5)(j2 + 30j + 1)) V^'-i) 

(33j5 + 32j'^ - AAAf + 1286j2 + 515j - 590) (i) 
~ 2((f + 30j + l)(-3 + j7(2j3 - 3j2 - 62j - 5)f''^^-'^ 
, (6/ + 7f - 17 f + 7j + 77) 

+ 6((2i3 _ 3^-2 _ 62j - 5){f + 30j + l))'^^^-!) I ■ ^) 

~ ^ (-13j'^-99j^ + 177j + 2131) (i) 
' 6((2j3 - 3j2 - 62j - 5)(j2 + 30j + 1)) oo^^-^) 

(2/ + 50^6 - 343j5 - 1342j4 + 5682j=^ - 2384j2 - 4869i + 6980) (i) 

2((j2 + 30j + l)(-3 + jY{2f - 3f - 62j - 5)) ^-(^"i) 
_ (31j^ - 33j^ + 71j + 907) (1) 

6((2j3-3j2-62j-5)(j2 + 30j + l)) ^' ' 

This way we get the background in the hmit where the back reactions from fluxes, 
branes, planes and strings are small. In the presence of off-diagonal energy-momentum 
tensor, the analysis will have to change but the underlying physics will remain un- 
changed. 
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B. Operator equations for metric fluctuations 



As we discussed in section 3.3, in the limit where the background fluxes, including 
the effects of the D7 brane, are very small the equation ( p.96| ) can be presented as 
an operator equation of the following form: 

« ^fjM ^ T^'r^x) (B.l) 

where x is a generic five- dimensional coordinate and T^J|f™^(x) is the energy momen- 
tum of the string. In this final form, the operator A^^ is a second order differential 
operator derived from ( p. 951) , with /x, i/, a, /3 = 0, 1, 2, 3, 4. 

Recall now that where we have used coordinate transformations to fix the five 
components of l^^, namely = 0,/i = 0,1,2,3,4. An additional residual gauge 
transformation allow us to eliminate another five degrees of freedom and we end up 
with five physical degrees of freedom i.e. five independent metric perturbation using 
which all other components could be expressed. Alternatively one can use certain 
combinations of the fifteen components to write five independent degrees of freedom 
for the metric fiuctuations. 

The above result is easy to demonstrate for the AdS space, as has already been 
discussed in [|^. For non-AdS spaces this is not so easy to construct. Therefore in 
the following we will take the complete set of ten components and using them we 
will determine the triangle operator A"^^ in (|3.96|) . If the ten components that can 
be labelled set: 



In 



{ho, hi, h2, hs, hi, h2, hs, hi, hs, hs} (B-2) 

then the operator A^^ in ( p.96|) gives rise to 77 equations that we present below. 
The warp factor h appearing in these equations can be taken to he h = h{r, vr, vr) 
because we will analyse fiuctuations close to the string and therefore our choice of 
background will be: 



Loo{t,r,x,y,z 
Loi{t,r,x,y,z 
Lo2{t,r,x,y,z 
Lmit,r,x,y,z 
Lnit,r,x,y,z 
Li2{t,r,x,y,z 



-g{r) + Khoit,r,x,y,z) 
Khi{t,r,x,y,z) 
Kh2{t,r,x,y,z)) 
Kh3it,r,x,y,z) 

h{rY/'^ + Khi(t,r,x,y,z) 
9{r) 

hih2(t,r,x,y,z) 
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Li3{t,r,x,y,z) 
L22{t,r,x,y,z) 
L23{t,r,x,y,z) 



Kli3{t,r,x,y,z) 
1 + Kl22{t,r,x,y,z) 
nl2z{t,r,x,y,z) 



r u \ + i^h3{t,r,x,y,z) 

L^sit, r, X, y, z) = ^^-p (B.3) 

Therefore using all the considerations, the explicit forms for A^^ can be presented 
as: 

^33 l_ Ag^2^|^ ^ g^^I ^ _ 10^2^3^^ + Aghi—\ 

16/ii V dr"^ \(?x2 y dr dr dr J 

1 / 9,6 9^ 9,3 9/1 (9 5 do d \ 

^" = "I^ + * & - ) 

°° 16/?J V drdy dr dy dr dy J 

All - --^ (-iQgh^ ^ 



16h 



dxdy ) 



Shg^y^^dtdxJ 
A°? = ^ ( 8h 



°3 - - J- /'s/,^^^ 

Shg"^ \ \dz^ dy"^ dx^ J dr dr dr dr J 



8/?.(?2 y (9^2 \(9a;2 (9^2 y dr dr dr 



^""^ 8/i5f2 y Qf2 + 4^ (7 (^^^2 dz"^) '^^ dr dr ^ '^^^ dr dr J 
^ "8^ ("^^'^a^) 



A IS 1 f d ^ ^ dg d \ 
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1 fQg2^_^ _ ^gl^^^\ 

^hg"^ \ dr dx dr dx J 



1 / dhd_\ 




or or or or \or J or or 

+ 8 ^h^ (— + — ^—]+ 12g^h— - 8 '/i^^ 




Idhig"^ i^^^ ^^^9r) ^ dr dr ^'^9' h 

1 / 3 dh 9 , 5 (9(7 (9 , 5 (9^ \ 

IQh^g^ \ or or oy oyor ) 

^i^^^ \ 9r 5r dr dr \dr J dr dr 




16/l2gr^ 



' o, sdg d „ s, 5 "9^ o, ^2,^9^ ^ (9/i (9 \ 

+ - sft'sgii + 89^'- a? - * * j 

16/12 gr2 V ^^^?/ 5r5yy 

1 f ,r d^ \ 

— IQh^g—— 

16h^g^ V dtdyj 

1 n, sdh d . „ ■i, 5 dq d , „ q, 5 (9^ \ 

(-^°« * * + + (aF - 9f) j 

1 Z^-. ^ f dh\^ y^adhdq n^zd'^h ^,5 oS^oA 

1 I Ag^K^^^ ^ +6g^h^^ ^ 24:hg^ ( + ISg^h^^^^ 
l&h^g"^ \ dr dr dr dr \dr J dr dr 

d'^a d'^h ( (9^ 9^ 92 ' 




1 / , 3 9 (9/i (9 , 5 (9g (9 , 5 9 (9^ \ 

^ -20/i2/ + Sgh^^— + IGh^g^^—- 

16h^g^\ drdt ^ dr dt ^ drdt J 
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16/i2gr2 



ar ar or or \or J or or 



— -16h^g^—— 




3 ^ 

^ dtdx J 



dr J dr'^ dr dr j 



\ dr dy drdy J 

\ dr drdt J 

\ dr dt J 



1 




16hH \dy' 



^^3h^ - ^Qh—— + Wgh—— + 8h — 
or dr dr dr dr dr 




16/i3 



drdx dr dx dr dx J 
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2 _ ^^^\ 
drdt dr dt j 



i2 



dxdy) 



^ , dh d , dhdq , nd q 
+ lOgh—— - IQh—^ + %h^^ 
or or or or or 



dhd\ 




or ay or oy ) 



dr dr 



drdx J 

3„2 



IGh^g'^ \ ^ ^ drdy ) 



IGh^g'^ \ ^ drdt^ 

1 r4/,3^^+8/.v^' 



IGh^g^ \ dr dx drdx J 

1 / 3 d^ \ 



img 
1 



^3g,2 y ^ dr dx ^ dr dx ) 



+ 12g h—-8hg—^+16hg -- 
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l&h^g"^ \ dxdy J 

A22 _ __}__ ( QU5„2 \ 

~ Idh^g^ ^ drdy J 

^ -Teh (^^'^^) 

aoo_ ^_(Ah^^9d _ 3 \ 

Idh^g^V^drdy ^ drdy J 

IGh'^g'^ \ ^ dr dy ^ dr dy ) 

or or or or or or or I 

16h-^g V dtdx J 

~ Wh^g V dtdy J 

IQh^g ^9 ^ drdx dr dx ^ dr dx) 

^12 ^ (^_8//,2 _ 8gh^^— + lOA— — ^ 

16h^g \ drdy dr dy dr dy J 



The above therefore summarises all the fluctuation operators for the OKS-BH back- 
ground. As one can see, the situation here is much more involved than the AdS-BH 
case. One might however try to simplify the 77 equations by imposing some symme- 
try in the background, much like the one that we discussed in the previous appendix. 
In the next appendix we will consider such a simplification for the OKS-BH back- 
ground by taking diagonal perturbations. 
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C. An example with diagonal perturbations in OKS-BH back- 
ground 



As we discussed in the previous section, the full analysis with all ten components 
of the metric fluctuations is rather difficult. In section 3.3 we did present a partial 
analysis taking all the components into account. However in that section we couldn't 
provide precise numerical answers to the metric fluctuations because of the underlying 
complexity of the problem. In this section we will take a middle path where we will 
only consider diagonal perturbations for the metric fluctuations much like what we 
did in Appendix A for the AdS case. This means that with the choice of coordinates 
and the string profile given by ( |3.57| ) and ( |3.58| ) we can now formally obtain the source 



in ( JB.ID provided we take into account not only the energy momentum tensor of the 



string but also other sources (see below). However since the string moves in the x 
direction, the perturbation created in y and z directions are equal and therefore we 
can demand: 

^22 = ^33 (C.l) 

This would mean that we have only four independent components of l^^, and it is 
enough to consider four independent equations sourced by the energy momen- 



tum tensors that can come from various sources like the string, fiuxes, D7 branes as 
well as 07 planes: 

ry-itotal _ ry-istring I ry-ifluxes I ry-iplanes I ry-ibranes /pi r\\ 

-'- HV = ^jiv ' fiu ' -'- fiu ' -'- tJ.f {^■^) 

At this stage we will assume that the total energy momentum tensor from all the 
above sources guarantee a diagonal perturbation in the system. We will also assume 
that the sources are all expressed in terms of the variable u where in terms of u, the 
UV is at M = whereas IR is at m = cxd. Therefore in the regime close to u = oo we 
expect certain aspects of QCD to be revealed from the gravity dual ( |3.101| ). 



We will now solve ( p.l|) order by order in ggNj, ggM /N and g^NjM. At zeroth 
order in gsNf, g'^NfM, ggM'^/N, the warp factor becomes h{u) = L^u^ and the metric 
( p.l01| ) reduces to that of AdS^. Hence at zeroth order in gsNj, ggM'^/N, g'^MNf our 
analysis will be similar to the AdS/CFT calculations As ( |B.1| ) is a second 

order non linear partial differential equation, we can solve it by Fourier decomposing 
X, y, z, t dependence of l^i, and writing it as a Taylor series in u in the following way: 

l,u = + 9sNf{a +bg,M)l^^l (C.3) 

where the subscript [0], [1] refer to the zeroth and the first order in {gsNf)J^ The 
zeroth order can then be succinctly presented as a Fourier series in the following 



^^Henceforth it also means we are keeping terms upto 0{gsM'^ /N, g'^NfM). 
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way 



77. 



,u,x,y,z) 



k=0 



dqidq2dq^duj 



-i{ujt-qix-q2y-q3z) (fe)[0] i 



(C.4) 



where s|it'''°' are expansion coefficients of tlie solution /[S- Tlie constant coefficients 
(a, 6) in /ji,! can be worked out easily. 

Similarly, we can also write the source in Fourier space as: 



T-itotal T-i[0]total 1 T-i[l]total 

in/ inj ~r in/ 



' fll/ 



(C.5) 



where as before, [0, 1] refer to the zeroth and first orders in (ggNf) respectively. The 
zeroth order can then be written as: 



dqidq2dq^duj 
(2^)^ 



(C.6) 



where t^^^l are expansion coefficients of source j^j^^^"^^^ at zeroth order in QgNf. These 
coefficients are obtained by using explicit expressions for T^°*^'(a;'^) given above in 
( |C.2|) . In terms of matrices we then expect the metric fluctuations la, i = 0, 1,2,3 
to take the following form: 



/ All All All \ / /oo \ 



^02 



^03 



A22 
^03 



00 All 



a33 
^02 



A33 
^03 



hi 



'■22 



Af^ / V /33 / 



1 



/ ^01 \ 

^02 



V J 



(C.7) 



where we have only switched on the 01, 02 and 12 components of the sources; and 
A^(^ are now deflned wrt the variable u instead of r, the radial coordinate. 

Manipulating the above matrix equation and using the explicit expressions for 
tf!l, we can extract a relation for all 522"*^°^ zeroth order in ggNj^ggM ( p.l| ) as: 



E 



3ku I dg 



3u 



• dg Q 



g du 2 



(C. 



^''Note that in all the subsequent mode expansions we will be ignoring the g{uc) dependences. 
Therefore for us Uc is close to the actual boundary so that r « t. A more careful analysis has been 
presented in section 3.3 wherein we took all the subtleties into account. 
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where A'^ can be given in terms of a series in in the following way: 



^0 = S{uj - vqi)9{u- Uo)^C^ 



j=4 



(C.9) 



once we 



where (j are w-independent constants that could be determined from 
have the explicit expressions for all the terms in ( |C2p . We can now use the above 
equations ( |C.8| ) and ( |C7| ) to write a relation between s^^"*'"' and 3^22^^^ ^^e following 
way: 



E 



.(fc)[o] 



,(fc)[o] 



'11 



+ 



"22 



u 



(C.IO) 



where again can also be expressed in series like (|U.9| ). The above two set of equa- 
tions have infinite number of variables. They can be solved if we know a generating 
function. Such a function could be determined for Sgg^'^^ in terms of and 3^22'^^^ 



as: 



.(fc)[o] 



'33 



C°(w,M,^) -gs[\ 



Jk)[0] 



'22 



(c.ii; 



Thus if we know s[^^'°^ and ^22^'°^ we can determine Sgg^'*^'. In the following we will 
determine these coefficients using Green's function. However before we go into it, let 
us write the last equation relating Sqq''^°' to the other coefficients: 



E 



-2qiu dg 



+ Mik 1 4?'°^ + I 2^7giT^ - I2g'qi 1 8^^^' - Sgq.ks 



dg 



g du 



du 



,(fc)[o] 
'33 



u 



k-l 



(C.12) 



where T)*^ are determined from the source (|C.2| ) like ( |C.9D above. Note that we 
could also write (|C.12| ) in terms of and s"22^'^^ using the generating function 

To solve the set of equations ( |C8| ) to ( |(J.12| ) we will be using Green's functions. 
Since all the equations are given in terms of series in u^, we first write the delta 
function as^^ 



u 



(C.13) 



so that we can equate coefficients on both sides of the equation. Using this it is 
straightforward to show that the the Green's function for (|C.8|) is given by: 



^°2(m,w,^) = ^c°(w,^; 



u 



(C.14) 



i=-4 



^^Such a way of expressing the delta function can be motivated from the standard completeness 
relation in quantum mechanics. Here of course the coefficient bi are some specified integers. 
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where c° are given, in terms of the bi coefficients appearing in the delta function 
( P.13|) , in the following way: 



2ibn 



2ibi + 86260 



^-2 



2ib3 + 2Aib2 + SQOibi + 6480z6o 
2^64 + ISibs + 216262 + 3240i6i + 58320i6o 



2^62 + 30z6i + 540260 



(C.15) 



Observe that the lower limit of the sum is from i = —4 because of the u ^ suppression 
in the LHS of 

It is now easy to write down the solution for the metric perturbation /f^' using 
the Green's function as: 



4 M, a;, z) 



fc=0 



duod^q 
duod^q 



-i{u)t-q\x-q2y-q-iz) 



duA\u)gUu\uj,t) 

ii'-t-q,x-q„,z)^mo]^^^ - uo)n^(C.16) 



where explicit expressions for 522'"°^ are given below. Observe that with 532''^"' known, 



(fc)[0] 



we obtain s^^'"^^ using (|C.10D . Then knowing 822''^'^', s^'^''^'^' we obtain 333''^'^' using ( |C.11|) 



(fc)[o] Jfc)[o] 



and finally Sqq'''^^ using ( |C.12|) . The explicit expressions for ^22''^'^' are: 



{km 



J0)[0] 
*22 



6{lj -vqi] 



- C4C°4«o + iCA + C5C°4)^ - (C4C° 2 + Csc" 3) "° 



3 



+ (C4C° 1 + C5C° 2)^ - (C4C° + C,C\ + Csc\)^ 



„{1)[0] 
''22 


= 6{io- 


^?l)C4C°4 


J2)[0] 
■^22 




-^gi)(C4C°3 + C5C°4) 


J3)[0] 
''22 




-t;gi)(C4C°2 + C5C°3) 


J4)[0] 
^22 


= \h. 


-^gi)(C4C°i+C5C°2) 


Am 
^22 




-^gi)(C4c[J + C5C°i+C8C°4 



(C.17) 



Observe also that all ^22^'°^ are proportional to 6{uj — qiv) which will eventually 
produce the Mach cone. Finally, using (|C.16|) in (|C.10|) , (|C.11|) , and ( |C.12|) , we can 
obtain rest of the metric perturbations at zeroth order in Qs (although we do not 
present the explicit expressions for the rest of the zeroth order perturbations here). 
Now we solve for the metric perturbation at linear order in gsNf,gsM where 



it is easiest to switch to coordinate u 



r-c(l-0' 



SO that the entire manifold is now 
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described by < C < 1 with Vc arbitrarily large and we can get a meaningful Taylor 
series expansion of the logarithms and other functions appearing in the equation. 
As we did for the zeroth order cases, we can decompose the first order in QgNf, QgM 
fluctuations via the following Fourier series: 



fc=0 



dqidq2dq^duj 



(C.18) 



where s^l?'^' are the corresponding Fourier modes. These Fourier modes will eventu- 
ally appear in the final equations for t^l. In fact we would also need the zeroth order 
perturbations /j?]. in the equations. Therefore we decompose as: 

dud^q 



,x,y,z) 



where / 



{km 



are now the Fourier modes for ll^l. These modes have one-to-one corre- 



spondence with given earlier and can be related by coordinate transformations. 
Using these equation ( P.7|) reads: 



7(fc)[0] , 

'■22 "T 



+ 



iujl 



(fc)[0] 



22 



(fc)[l]r4 
22 ^ 



1 dg 



h dg 3 dh ~dg ghr^ 
AgTc'AglC^ d( ^ ^ 



(C.20) 



rtil - CY 



3 ^ ]_dg^ ^ g_ 

Ar^gdC r,(l-C)^ r^dC 2 



'22 I22 



{km 



(C.21) 



where as before as in ( P.8|) , we could separate the equations relating sh 
from the rest of the other Fourier modes. We have also defined h as: 

^ = r\l- CY ["^ ""'^^ -0-B log2 r,(l - C) 

with A and B are defined in (|3.101|) . The other variables appearing above have 
already been defined earlier. Using these, and using the appropriate Green's function 
we can easily determine these Fourier modes. For example the explicit expressions 
for 522""^^ can be determined by first writing the delta function as: 

5{l/r) 



i=0 



(C.22) 



with hj defined in the following way: 



h = E 



(C.23) 



i=0 
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Then solve for the Greens function for (|C.20 ) to obtain 

oo 



(C.24) 



1=1 



where the coefficients appearing above are defined as: 
cl = -4?nAlog(r,) - 51og- - 



-Ig^^'^ [Alog(r,) - i?log^(re)] - ^ 

-4c} - 4)'°] [4A\og{r,) - 4B\og\r,) -A + 2Slog(r,)] - 2lg^^'^ [Alogir,) - B\og\r,)] 
2cK r rj _ 3 _ 4r| ^1 24y°Vj [Alog(r,) - glog^(r,)] ^ 



12 [Alog(rc) - 51og^ 



+ 



L^^(r,) [Alog(r,) - Elog^ 
2 



3A + 651og(rJ ^ 4 [Alog( 



r^g[rc) 
51og^(rc)l rl 



(C.25) 



Once everything is laid up, we can write the source in ( |C.2CI|) as a power series in 
C i.e. = SjC"' with dj derivable from (|C.2|) . This would finally give us the 

..... _ 



required Fourier coefficients to first order in QsNj as^^: 



Mil] 


= 6{u 


- vqi 


J2)[l] 
^22 


= 6{u 


- vqi 


J3)[l] 
"'22 


= 6{u 


- vqi 



c}aoCo (4^0 + qai)Co 



(C.26) 



with the following Fourier decomposition: 



dud^q 



-i(ujt-qix-q2y-q3z) 



= J2 I ^^e-*('^*-«i^-«^J^-^^^)4?W(a;,^)^^(C-Co)C'(C.27) 



k=0 



i2ny 



Once we know these modes, we can use them to write the relation for s\i in the 
following way: 



E 



, ^22 



'11 



9 



(C.28) 



79 



We are only solving upto 0(C^ 
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with given in the appendix. Observe that the above equation has exactly the same 



form as (|C.10|) except that these modes are written for hnear order perturbations. 
It is then no surprise that the generating function for Sgg^^^' takes exactly the same 
form as in ( |U.11|) : 



'33 



'22 



(C.29) 



Finally once we know the Fourier modes s-^ with i = 1,2,3 we can use them to 
write the equation for Sqo^'^'. The equation turns out to be rather involved with all 
zeroth and first order coefficients appearing together. Nevertheless one can present 
the following form for the equation: 



E 



+4q,hlg^^''he-'- + 4q,L 



-2 



■ dC Teg 



00 



dg L^qi 
dCrlil-O'g 



'00 ^ 



fc-1 



do h 



(km^k 



11 



+ 2g 



dg 



dcr^A^-cr 



,dh 1 



{k)[0].k 



-12g' 



r5(l-C)5 



{km,.k-i 



33 



4 '='33 



(C.30) 



where V^,C^ is determined from the source (|C.2|) . Plugging in the other Fourier 



modes, we have been able to solve for all the Sqq''^^' modes using the corresponding 
Green's function (we don't present the results here). 

To conclude therefore, using the equations ( |C.28| ), ( |C29| ) and ( |(J.3CI| ) we find 



Jk)[l] Jk)[l] 



'11 



and 533'*'^' etc. which in turn give us lQQ{t,(,x,y, z),lii{t,(,x,y, z) and 



l\^^{t, (, x,y, z) etc. using equations like ( |(J.27] ) 
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D. Coefficients in (|3.116|) and ( ^.12U|) 

A sample of the coefficients appearing in ( p.l2(]| ) are given below for the simplified 
OKS-BH geometry with approximate diagonal perturbations: 

„~2<h 



402 _ 420 _ ^ (^f\h'^n^\ 



,,,,, (16(7=^/12) 

g-2</. 

-■-32^(--V|.24.f«') 
'^J" - 32^ h'V') 

- 32^ - 

- 32^^ h'V') 

-"-3£^K«'--«""'f) 

-3^^ --V|) 
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32^/113/4 \^ dr^ ^ dr 



where h is the warp factor measured on the shce ( |3.14| ), and is the dilaton. To get 



the exphcit expressions for Aq,Cq in ( ^.1161 ) we need to replace h with L'^/r'^ and 



e"^"^ with l/g1 in the above expressions. We have written the nonzero A^l^Ci and 
therefore the terms not appearing above are all zeroes. Note also that A^j^ = B^^ , k = 
0, 1 and since we don't explicitly know the sources and the full background, we don't 
know the explicit expressions for the terms D]^ , El, etc. 
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E. Detailed Viscosity Analysis 



As in the previous section, here we work out the coefficients in the quadratic action 
in (gTTg) upto 0(g,N,,g.M^IN): 



9s V 27r IT J 



^.2 



271 



TT 



r' I + 24/ - + rt f + 150/ f ^ | log. 
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(logr) 2 - 8 
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2tx 
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31og(r) - 2 (g'{r) + Mil 



logr + 12 logr + 2 
3 



27r 



+ 



27rr 



r/i^ 



J{r) = - [c?o + c/ilog(r) + rf2(log(r))2] + _ [cq + eilog(r) + rf2(log(r))2] + 



where J7'(r) and M.{r) appear in (|3.173|) . The coefficients appearing in the equations 
above are defined as: 



729glM^a'^ 75g',M^Nf 122^ 
do--AgsNf+ ^g^4^2^r2 
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ei 
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2\ 2 
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„2 »,,f2 
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Finally the perturbation coefficients ai,bj appearing in the mode expansions for Q 
and n in (|3l8^) and (gl84D respectively are given to 0{gsNf, gsM'^ /N) as: 
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F. Detailed Entropy Analysis 

The coefficients in ( ^.216| ), imposing an UV cut-off at r = Tc, are: 

Cl' = r^e-''^ ( gQ^^ 2 ^ 8^2^ 2 -1^ _ 4^2 2^ _ 48^2^^\ 

^ f-240(7'/.r + U^g'^r' - 24r'hg'g^'^ + 24rV^] 

32^/^hg^L^r^hy^ V dr^ " dr ^ dr J 

Cl' = , , (-2AQg'hr + ufgh^ - 2ArV92'^ + ^Sr^hg^^) 

32^/^hg^L^r^h^l/A) \ dr dr dr J 

Cf = rAVL...u (-^^^9'hr + 12 fg'r^ - 2ArV^ - 2Ar^hg%''^f) 
32^/^hg^L^r^h^/^ \ dr dr dr J 



^6g-20 
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Fl = —L^^^^ i-2Agh^h) 

The nonzero coefficients in ( p.219| ) are (taking the approximation gi = g2 = ^ — -r + 
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The other coefficients &„(j)[i] in p.219| ) are defined as: 
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where A = f i + M 



N 



4. /' ^ = ^ISr^' a = 1 + a = 1 + We also have 



J0)[1] 
'00 



.{4)[1] _ 



'00 
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■ (Alog r + i^log^r) 
(Alog r + i?log^r) 



—Alog r — i?log r 



with every other ^' = 0. Note that (a) all coefficients are suppressed as 0{gsNf, ggM'^/N) 
as expected, and (b) an appropriate UV cap will make all Smn^^ independent of r, as 
we discussed in the main text. 
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